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I Abstract 

PQ . Aims. The largest Local Group spiral galaxy, M 31, has been completely imaged for the first time obtaining a luminosity lower limit 

HH ■ ~10^^ erg s~^ in the 0.2-4.5 keV band. Our XMM-Newton EPIC survey combines archival observations along the major axis, from 

' June 2000 to July 2004, with observations taken between June 2006 and February 2008 that cover the remainder of the D25 ellipse. 

^ , The main goal of the paper is a study of the X-ray source population of M 3 1 . 

O |. Methods. An X-ray catalogue of 1 897 sources was created; 914 were detected for the first time. Source classification and identification 

I ■ were based on X-ray hardness ratios, spatial extent of the sources, and by cross correlating with catalogues in the X-ray, optical, 

O ' infrared and radio wavelengths. We also analysed the long-term variability of the X-ray sources; this variability allows us to distinguish 

^ I between X-ray binaries and active galactic nuclei (AGN). Furthermore, supernova remnant classifications of previous studies that did 

^ . not use long-term variability as a classification criterion, could be validated. Inclusion of previous Chandra and ROSAT observations 

, ^ , ■ in the long-term variability study allowed us to detect additional transient or at least highly variable sources, which are good candidate 
X-ray binaries. 

I , Results. Fourteen of the 30 supersoft source (SSS) candidates represent supersoft emission of optical novae. Many of the 25 supernova 

^ . remnants (SNRs) and 31 SNR candidates lie within the lOkpc dust ring and other star forming regions in M 31; this connection 

ly-^ ■ between SNRs and star forming regions implies that most of the remnants originate in type II supernovae. The brightest sources in 

' X-rays in M 3 1 belong to the class of X-ray binaries (XRBs). Ten low mass XRBs (LMXBs) and 26 LMXB candidates were identified 

I based on their temporal variability. In addition 36 LMXBs and 17 LMXB candidates were identified due to correlations with globular 

. clusters and globular cluster candidates. From optical and X-ray colour-colour diagrams, possible high mass XRB (HMXB) candidates 

• ■ were selected. Two of these candidates have an X-ray spectrum as is expected for an HMXB containing a neutron star primary. 

' Conclusions. While our survey has greatly improved our understanding of the X-ray source populations in M 31, at this point 65% of 

\ the sources can still only be classified as "hard" sources, /. e. it is not possible to decide whether these sources are X-ray binaries or 

, Crab-like supernova remnants in M 31, or X-ray sources in the background. Deeper observations in X-ray and at other wavelengths 

■ would help classify these sources. 
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1. Introduction lite (Giacconi et al. 1974), the imaging X-ray optics flown on 
^ .,1.1 .11 1 A 1 , the Einstein X-ray observatory permitted the resolution of in- 
Our nearest neighbounng large spiral galaxy the Andromeda ^^^j^^^j ^.^ay sources in M 31 for the first time. In the en- 
galaxy, also known asM 31 orNGC 224, is an ideal target for an ^.^^^^.^ j observations of M 31, Trinchieri 
X-ray source population study of a galaxy similar to the Milky ^ p^j^j^^^^^ ^1991^ j^^^^^^^^^ ^P91^ ^^^^^ individual X-ray 
Way. Its proximity (distance 780 kpc, Holland 1998; Stanek & ^^^^^^^ ^^ ^^^^ ^j^^^ ^ ^ 4^^q36 ^-i^ ^f which 16 sources 
Gamavich 1998) and the moderate Galactic foreground absorp- ^^^^^^ variability (van Speybroeck et al. 1979; Collura et al. 
tion(A?H=7xl02°cm"^ Stark etal. 1992) allow a detailed study ^^^qx fj 
of source populations and individual sources. 

After early detections of M 31 with X-ray detectors mounted juiy 1990^ the bulge region of M 31 was observed with 

on rockets (e.g. Bowyer et al. 1974) and the Uhuru satel- t^e ROSAT High Resolution Imager (HRI) for ~ 48 ks. Primini 

~ Z ■ . . , TT 0.1 -1 et al. (1993, hereafter PFJ93) reported 86 sources brighter than 

htlter.sSLeS " ™' ^ l.^xlO^ erg s- m tMs observation. Of the ROSAT HRI 

* Based on observations obtained with XMM-Newton, an ESAsci- ^^^^^^^ ^^^^^ed within 7.5 of the nucleus, 18 sources were 

ence mission with instruments and contributions directly funded by ^^nd to vary when compared to previous Einstein observa- 

ESA Member States and NASA. lions and about three of the sources may be "transients". Two 

Tables 5 and 8 are only available in electronic form at the deep PSPC (Position Sensitive Proportional Counter) surveys 

CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via of M 31 were performed with ROSAT, the first in July 1991 

http://cdsweb.u-strasbg.fr/cgi-bin/qcat7J/A-hA/ (Supper et al. 1997, hereafter SHP97), the second in July/August 
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1992 (Supper et al. 2001, hereafter SHL2001). In total 560 X- 
ray sources were detected in the field of M 31; of these, 491 
sources were not detected in previous Einstein observations. In 
addition, a comparison with the results of the Einstein survey 
revealed long term variability in 18 sources, including 7 possi- 
ble transients. Comparing the two 7?05'Ar surveys, 34 long term 
variable sources and 8 transient candidates were detected. The 
derived luminosities of the detected M 3 1 sources ranged from 
5x10^^ erg s"^ to 5x10^^ erg s"^. Another important result ob- 
tained with ROSAT was the establishment of supersoft sources 
(SSSs) as a new class of M 31 X-ray sources (cf. Kahabka 1999) 
and the identification of the first SSS with an optical nova in 
M 31 (Nedialkov et al. 2002). 

Garcia et al. (2000) reported on first observations of the nu- 
clear region of M 31 with Chandra. They found that the nu- 
clear source has an unusual X-ray spectrum compared to the 
other point sources in M 31. Kong et al. (2002b) report on eight 
Chandra ACIS-I observations taken between 1999 and 2001, 
which cover the central ~ 17' x 17' region of M 31. They de- 
tected 204 sources, of which -50% are variable on timescales 
of months and 13 sources were classified as transients. Kaaret 
(2002) detected 142 point sources (Lx = 2x10^^ to 2x10^^ 
erg s"^ in the 0.1-lOkeV band) in a 47 ks Chandra/HRC ob- 
servation of the central region of M 31. A comparison with a 
ROSAT observation taken 1 1 yr earlier, showed that 46+26% of 
the sources with Lx > 5x10^^ erg s"^ are variable. Three dif- 
ferent M 31 disc fields, consisting of diff'erent stellar population 
mixtures, were observed by Chandra. Di Stefano et al. (2002) 
investigated bright X-ray binaries (XRBs) in these fields, while 
Di Stefano et al. (2004) examined the populations of supersoft 
sources (SSSs) and quasisoft sources (QSSs), including obser- 
vations of the central field. Using Chandra HRC observations, 
Williams et al. (2004) measured the mean fluxes and long-term 
time variability of 166 sources detected in these data. Voss & 
Gilfanov (2007) used Chandra data to examine the low mass 
X-ray binaries (LMXBs) in the bulge of M 31. Good candi- 
dates for LMXBs are the so-called transient sources. Studies 
of transient sources in M 31 are presented in numerous papers, 
e. g. WilHams et al. (2006b), Trudolyubov et al. (2006, hereafter 
TPC06), WilHams et al. (2005b), WilHams et al. (2006a, here- 
after WGM06), and Voss et al. (2008). 

Using XMM-Newton and Chandra data, Trudolyubov & 
Priedhorsky (2004) detected 43 X-ray sources coincident with 
globular cluster candidates from various optical surveys. They 
studied their spectral properties, time variability and log N-log S 
relations. 

Osborne et al. (2001) used XMM-Newton Performance 
Verification observations to study the variability of X-ray 
sources in the central region of M 31. They found 116 sources 
brighter than a limiting luminosity of 6x10^^ erg s"^ and exam- 
ined the ~ 60 brightest sources for periodic and non-periodic 
variability. At least 15% of these sources appear to be variable 
on a time scale of several months. Barnard et al. (2003) used 
XMM-Newton to study the X-ray binary RX J0042.6+41 15 and 
suggested it as a Z- source. Orio (2006) studied the population 
of SSSs and QSSs with XMM-Newton. Recently, Trudolyubov 
& Priedhorsky (2008) reported the discovery of 217s pulsations 
in the bright persistent SSS XMMU J004252.5+411540. Shaw 
Greening et al. (2009, hereafter SBK2009) presented the results 
of a complete spectral survey of the 335 X-ray point sources they 
detected in five XMM-Newton observations located along the 
major axis of M 31. They obtained background subtracted spec- 
tra and lightcurves for each of the 335 X-ray sources. Sources 
with more than 50 source counts were individually spectrally fit- 



ted. In addition, they selected 18 HMXB candidates, based on a 
power law photon index of0.8<r<1.2. 

Pietsch et al. (2005b, hereafter PFH2005) prepared a cata- 
logue of M 3 1 point-like X-ray sources analysing all observa- 
tions available at that time in the XMM-Newton archive which 
overlap at least in part with the optical D25 extent of the galaxy. 
In total, they detected 856 sources. The central part of the galaxy 
was covered four times with a separation of the observations 
of about half a year starting in June 2000. PFH2005 only gave 
source properties derived from an analysis of the combined ob- 
servations of the central region. Source identification and clas- 
sification were based on hardness ratios, and correlations with 
sources in other wavelength regimes. In follow-up work, (i) 
Pietsch & Haberl (2005) searched for X-ray burst sources in 
globular cluster (GIG) sources and candidates and identified two 
X-ray bursters and a few more candidates, while (ii) Pietsch et al. 
(2005a, hereafter PFF2005) searched for correlations with opti- 
cal novae. They identified 7 SSSs and 1 symbiotic star from the 
catalogue of PFH2005 with optical novae, and identified anaddi- 
tional XMM-Newton source with an optical nova. This work was 
continued and extended on archival Chandra HRC-I and ACIS-I 
observations by Pietsch et al. (2007, hereafter PHS2007). 

Stiele et al. (2008, hereafter SPH2008) presented a time vari- 
ability analysis of all of the M 31 central sources. They detected 
39 sources not reported at all in PFH2005. 21 sources were de- 
tected in the July 2004 monitoring observations of the low mass 
X-ray binary RX J0042.6-h4115 (PI Barnard), which became 
available in the meantime. Six sources, which were classified as 
"hard" sources by PFH2005, show distinct time variability and 
hence are classified as XRB candidates in SPH2008. The SNR 
classifications of three other sources from PFH2005 had to be 
rejected due to the distinct time variability found by SPH2008. 
Henze et al. (2009a) reported on the first two SSSs ever dis- 
covered in the M 3 1 globular cluster system, and Henze et al. 
(2009b) discussed the very short supersoft X-ray state of the 
classical nova M31N 2007-1 la. A comparative study of super- 
soft sources detected with ROSAT, Chandra and XMM-Newton, 
examining their long-term variability, was presented by Stiele 
et al. (2010). 

An investigation of the log N-log S relation of sources de- 
tected in the 2.0-10.0 keV range will be presented in a forth- 
coming paper (Stiele et al. 2011 in prep.). In this work the con- 
tribution of background objects and the spatial dependence of 
the log N-log S relations for sources of M 3 1 is studied. 

In this paper we report on the large XMM-Newton survey 
of M 31, which covers the entire D25 ellipse of M 31, for the 
first time, down to a limiting luminosity of ~10^^ erg s"^ in 
the 0.2-4.5 keV band. In Sect. 2 information about the obser- 
vations used is provided. The analysis of the data is presented 
in Sect. 3. Section 4 presents the combined colour image of 
all observations used. The source catalogue of the deep XMM- 
Newton survey of M 31 is described in Sect. 5. The results of 
the temporal variability analysis are discussed in Sect. 6. Cross- 
correlations with other M 31 X-ray catalogues are discussed 
in Sect. 7, while Sect. 8 discusses cross-correlations with cata- 
logues at other wavelengths. Our results related to foreground 
stars and background sources in the field of M 3 1 are presented 
in Sect. 9. Individual source classes belonging to M31 are dis- 
cussed in Sect. 10. We draw our conclusions in Sect. 11. 

2. Observations 

Figure 1 shows the layout of the individual XMM-Newton ob- 
servations over the field of M 31. The observations of the "Deep 
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Table 1. A selection of important X-ray surveys of M 3 1 . 



Paper 


S+ 


#ofSrc* 


Lt 

^2 -1 

erg cm s 


field 


comments 


Trinchieri & Fabbiano (1991, TF91) 


E 


108 


6.4x10^6-1.3x103^ 
(0.2-4 keV) 


entire set of Einstein 
imaging observations 


16 sources showed variability 


Primini etal. (1993, PFJ93) 


R (HRI) 


86 


> 1.8x10^6 
(0.2-4 keV) 


bulge region 


18 sources variable; ~3 transients 


Supper etal. (1997, 2001) 


R (PSPC) 


560 


5x10^5-5x1038 


whole galaxy 


two deep surveys 


(SPH97, SHL2001) 






(0. 1-2.4 keV) 




491 sources not detected with Einstein 

1 1 sources variable, 7 transients compared to Einstein 

34 sources variable, 8 transients between ROSAT surveys 


Osborne et al. (2001) 


X 


116 


> 6x1035 
(0.3-12 keV) 


centre 


examined the ~ 60 brightest sources for variability 


Kong et al. (2002b) 


C (ACIS-I) 


204 


> 2xl035 
(0.3-7 keV) 


central ~ 17' x 17' 


observations between 1999 and 2001 

~50% of the sources are variable, 13 transients 


Kaaret (2002) 


C (HRC) 


142 


2x1035-2x1038 
(0.1-lOkeV) 


centre 


one 47 ks observation; 46±26% of the sources 
with Lx > 5xl036 erg s"^ are variable 


Di Stefano et al. (2002) 


C (ACIS-I/S) 


28 


5x1035-3x1038 
(0.3-7 keV) 


3 disc fields 


bright X-ray binaries 


Di Stefano et al. (2004) 


C (ACIS-S S3) 


33 




3 disc fields + centre 


supersoft sources and quasisoft sources 


Williams et al. (2004) 


C (HRC) 


166 


1.4x1036-5x1038 
(0.1-lOkeV) 


major axis + centre 


^25% showed significant variability 


Trudolyubov & Priedhorsky (2004) 


C,X 


43 


~ 1035— 1039 

(0.3-10 keV) 


major axis + centre 


globular cluster study 


Pietsch et al. (2005b, PFH2005) 


X 


856 


4.4x1034-2.8x1038 
(0.2-4.5 keV) 


major axis + centre 


source catalogue 


Pietsch et al. (2005a, PFF2005) 


C, R,X 


21 


- 1035-^ 1038 

(0.2-1 keV) 


centre 


correlations with optical novae 


Orio (2006) 


C,X 


42 


6x1035— 1039 
(0.2-2 keV) 
(0.3-10 keV) 


major axis + centre 


supersoft sources and quasisoft sources 


Pietsch et al. (2007, PHS2007) 


C,X 


46 


~ 1035— 1038 

(0.2-1 keV) 


centre 


correlations with optical novae 


Voss & Gilfanov (2007) 


c 


263 


5xl033-1.5xl038 
(0.5-8 keV) 


bulge region 


low mass X-ray binary study 


Stiele et al. (2008, SPH2008) 


X 


39 
300 


7x1034-6x1037 
4.4xl034-2.8xl038 
(0.2-4.5 keV) 


centre 


re-analysis of archival and new 2004 observations 
time variability analysis; 149 sources with a significance 
for variability >3; 6 new X-ray binary candidates, 
3 supernova remnant classifications were rejected 


Shaw Greening et al. (2009, SBK2009) 


X 


335 


~ 1034— 1039 
(0.3-10 keV) 


5 fields along 
major axis 


background subtracted spectra and lightcurves for 
each source; 18 HMXB candidates, selected from their 
power law photon index 


Stiele et al. (2010) 


X 


40 




whole galaxy 


supersoft sources; comparing ROSAT, Chandra and 
XMM-Newton catalogues 



Notes: 

+ : X-ray satellite(s) on which the study is based: E for Einstein, R for ROSAT, C for Chandra, and X for XMM-Newton (EPIC) 
* : Number of sources 

^ : observed luminosity range in the indicated energy band, assuming a distance of 780 kpc to M 31 



XMM-Newton Survey of M 31" (PI Pietsch) mainly point at 
the outer parts of M 31, while the area along the major axis is 
covered by archival XMM-Newton observations (Pis Watson, 
Mason, Di Stefano). To treat all data in the same way, we 
re-analysed all archival XMM-Newton observations of M 31, 
which were used in Pietsch et al. (2005b). In addition we in- 
cluded an XMM-Newton target of opportunity (ToO) observa- 
tion of source CX0M31 J004059.2-h41 1551 and the four obser- 
vations of source RX J0042.6-h41 15 (PI Barnard). 

All observations of the "Deep XMM-Newton Survey of 
M 31" and the ToO observation were taken between June 2006 
and February 2008. All other observations were available via 
the XMM-Newton Data Archive^ and were taken between June 
2000 and July 2004. 

The journal of observations is given in Table 2. It includes 
the M 31 field name (Column 1), the identification number (2) 
and date (3) of the observation and the pointing direction (4, 5), 
while col. 6 contains the systematic off'set (see Sect. 3.4). For 
each EPIC camera the filter used and the exposure time after 
screening for high background is given (see Sect. 3.1). 



^ http://xinm.esac.esa.int/xsa/ 



3. Data analysis 

In this section, the basic concepts of the X-ray data reduction 
and source detection processes are described. 

3.1. Screening for high background 

The first step was to exclude times of increased background, 
due to soft proton flares. Most of these times are located at the 
start or end of an orbit window. We selected good time inter- 
vals (GTIs) - intervals where the intensity was lower than a cer- 
tain threshold - using 7-15 keV light curves constructed from 
source-free regions of each observation. The GTIs with PN and 
MOS data were determined from the higher statistic PN light 
curves. Outside the PN time coverage, GTIs were determined 
from the combined MOS light curves. For each observation, the 
limiting thresholds for the count rate were adjusted individually; 
this way we avoided cutting out short periods (up to a few hun- 
dred seconds) of marginally increased background. Short peri- 
ods of low background, which were embedded within longer pe- 
riods of high background, were omitted. For most observations, 
the PN count rate thresholds were 2-8 cts ks"^ arcmin"^. 
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Table 2. XMM-Newton log of the Deep Survey and archival M 31 observation overlapping with the optical D25 ellipse. 



M31 field 


Obs. id. 


Obs. dates 


Pointing direction 


Offset * 


EPIC PN 


EPIC MOSl 


EPIC M0S2 










RA/Dec (J2000) 




Filter 


^ exp 


Filter^ 


^ exp 


Filter^ 


^ exp 


(1) 




(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


Centre 1 


(cl) 


0112570401 


2000-06-25 


0:42:36.2 


41:16:58 


-1.9, +0.1 


medium 


zJ.4o(zJ.4oj 


medium 


2^.04(2^.04) 


medium 


zy. 04(2^.04) 


L,enire z 


(CZj 


0112570601 


2000-12-28 


0:42:49.8 


41:14:37 


-2.1, +0.2 


medium 


J.OZ.{ J.OZj 


medium 


A Al( A Al\ 
0.4Zi^ 0.4Zj 


medium 


A Al( A Al\ 
0.4Zi^ 0.4Zj 


Centre 3 


(c3) 


0109270101 


2001-06-29 


0:42:36.3 


41:16:54 


-3.2,-1.7 


medium 


Zi./i(Zi./ij 


medium 


11. Q'id'X Q<^ 


medium 


11 Qf^dl QA^ 
ZJ.oO(ZJ.oOj 


Ml 


(nij 


0109270701 


2002-01-05 


0:44:08.2 


41:34:56 


-0.3, +0.7 


medium 


/I fi ^AQ ^ 1 \ 

46. J 1(^46. J 1 ) 


medium 


J J.Ool^J J.Oo ) 


medium 


J J.O / (^J J.O / ) 


L,enire 4 


(C4j 


0112570101 


2002-01-06/07 


0:42:50.4 


41:14:46 


-1.0, -0.8 


thin 


4 / .OJi^4 / .OJ ) 


thin 


'\1 ^IC^l 87^ 
jZ.o lyjZ.o 1 ) 


thin 


jZ.oOi^jZ.oOJ 


Q 1 


(si) 


0112570201 


2002-01-12/13 


0:41:32.7 


40:54:38 


-2.1,-1.7 


thin 


40. / 3(^40. / J) 


thin 


Ji.oJl^Ji.oJj 


thin 


Ji.o4(^ji.o4j 


QO 

oz 


(s2) 


0112570301 


2002-01-24/25 


0:40:06.0 


40:35:24 


-1.1,-0.3 


thin 


ZZ.Zj{ZZ.Zj ) 


thin 


lA IKIA H\ 
Zh.Zj{ZH.Zj ) 


thin 


lA lAClA 1A\ 
Z^.Z^yZ^.Z^ ) 


N2 


(nzj 


0109270301 


2002-01-26/27 


0:45:20.0 


41:56:09 


-0.3,-1.5 


medium 


ZZ. 1 j\ZZ. Id) 


medium 


K llfK 11\ 
ZD .ZZyZj .ZZ ) 


medium 


l'\ lSifl'\ 1^\ 

zj.zoyzj.zo ) 


iN J 


(njj 


0109270401 


2002-06-29/30 


0:46:38.0 


42:16:20 


-2.3,-1.7 


medium 




medium 


4J. jU(^4J. jUj 


medium 


Al (^KAl (^1\ 
4J.0j(^4J.0j ) 


rl4 


(h4) 


0151580401 


2003-02-06 


0:46:07.0 


41:20:58 


+0.3, +0.0 


medium 


iU. i4mj. 14 ) 


medium 


LZ. /0{LZ. /Oj 


medium 


LZ. /0{LZ. /Oj 


RX 1 


(bl)* 


0202230201 


2004-07-16 


0:42:38.6 


41:16:04 


-1.3,-1.2 


medium 


io.3z(io.3z) 


medium 


i9.Zi(iy.Zi) 


medium 


in oi/in oi\ 
i9.Zi(iy.Zi) 


RX 2 


(b2) 


0202230301 


2004-07-17 


0:42:38.6 


41:16:04 


-1.0, -0.9 


medium 


0.0(0.0) 


medium 


0.0(0.0) 


medium 


0.0(0.0) 


RX 3 


(b3)* 


0202230401 


2004-07-18 


0:42:38.6 


41:16:04 


-1.7,-1.5 


medium 


12.30(12.30) 


medium 


17.64(17.64) 


medium 


17.68(17.68) 


RX 4 


(b4)* 


0202230501 


2004-07-19 


0:42:38.6 


41:16:04 


-1.4,-1.8 


medium 


7.94(7.94) 


medium 


10.12(10.12) 


medium 


10.13(10.13) 


S3 


(s3) 


0402560101 


2006-06-28 


0:38:52.8 


40:15:00 


-3.1,-3.0 


thin 


4.99(4.99) 


medium 


6.96(6.96) 


medium 


6.97(6.97) 


SSI 


(ssl) 


0402560201 


2006-06-30 


0:43:28.8 


40:55:12 


-4.4,-3.7 


thin 


14.07(9.57) 


medium 


24.56(10.65) 


medium 


24.58(10.66) 


SNl 


(snl) 


0402560301 


2006-07-01 


0:40:43.2 


41:17:60 


-2.7,-1.5 


thin 


41.23(35.42) 


medium 


47.60(39.40) 


medium 


47.64(39.44) 


SS2 


(ss2) 


0402560401 


2006-07-08 


0:42:16.8 


40:37:12 


-1.2,-1.3 


thin 


21.64(9.92) 


medium 


25.59(11.04) 


medium 


25.64(11.05) 


SN2 


(sn2) 


0402560501 


2006-07-20 


0:39:40.8 


40:58:48 


-0.8,-0.7 


thin 


48.79(21.45) 


medium 


56.13(23.85) 


medium 


56.17(23.86) 


SN3 


(sn3) 


0402560701 


2006-07-23 


0:39:02.4 


40:37:48 


-0.9, -2.0 


thin 


23.80(15.43) 


medium 


28.02(17.16) 


medium 


28.04(17.17) 


SS3 


(ss3) 


0402560601 


2006-07-28 


0:40:45.6 


40:21:00 


-1.8,-1.7 


thin 


27.77(20.22) 


medium 


31.92(22.49) 


medium 


31.94(22.5) 


S2 


(s21) 


0402560801 


2006-12-25 


0:40:06.0 


40:35:24 


-1.6, -0.7 


thin 


39.12(39.12) 


medium 


45.19(45.19) 


medium 


45.21(45.21) 


NNl 


(nnl) 


0402560901 


2006-12-26 


0:41:52.8 


41:36:36 


-1.5,-1.5 


thin 


37.9(37.9) 


medium 


43.08(43.08) 


medium 


43.1(43.1) 


NSl 


(nsl) 


0402561001 


2006-12-30 


0:44:38.4 


41:12:00 


-1.0,-1.3 


thin 


45.11(45.11) 


medium 


50.9(50.9) 


medium 


50.93(50.93) 


NN2 


(nn2) 


0402561101 


2007-01-01 


0:43:09.6 


41:55:12 


-0.0,-1.2 


thin 


41.73(41.73) 


medium 


46.45(46.45) 


medium 


46.47(46.47) 


NS2 


(ns2) 


0402561201 


2007-01-02 


0:45:43.2 


41:31:48 


-2.3,-1.7 


thin 


34.96(34.96) 


medium 


40.55(40.55) 


medium 


40.58(40.58) 


NN3 


(nn3) 


0402561301 


2007-01-03 


0:44:45.6 


42:09:36 


-1.4, -0.7 


thin 


31.04(31.04) 


medium 


34.81(34.81) 


medium 


34.81(34.81) 


NS3 


(ns3) 


0402561401 


2007-01-04 


0:46:38.4 


41:53:60 


-2.1, +0.3 


thin 


39.41(39.41) 


medium 


45.50(45.50) 


medium 


45.52(45.52) 


N2 


(n21) 


U4UZjoi DUi 


ZUU/-Ui-Uj 


0:45:20.0 


41:56:09 


—2.6, —1.3 


thin 


37.18(37.18) 


medium 


41.98(41.98) 


medium 


42.03(42.03) 


SSI 


(ssll) 


0505760201 


2007-07-22 


0:43:28.8 


40:55:12 


-2.5,-2.6 


thin 


30.07(23.90) 


medium 


34.01(26.70) 


medium 


34.02(26.72) 


S3 


(s31) 


0505760101 


2007-07-24 


0:38:52.8 


40:15:00 


-1.8,-1.0 


thin 


21.86(15.74) 


medium 


24.74(17.65) 


medium 


24.74(17.65) 


CX0M31 


(snl If 


0410582001 


2007-07-25 


0:40:59.2 


41:15:51 


-1.2, -0.3 


thin 


11.27(11.27) 


medium 


14.01(14.01) 


medium 


14.02(14.02) 


SS3 


(ss31) 


0505760401 


2007-12-25 


0:40:45.6 


40:21:00 


-1.0, +0.1 


thin 


23.56(22.82) 


medium 


28.18(25.8) 


medium 


28.2(25.82) 


SS2 


(ss21) 


0505760301 


2007-12-28 


0:42:16.8 


40:37:12 


+ 1.3,-0.1 


thin 


35.28(35.28) 


medium 


40.00(40.00) 


medium 


40.01(40.01) 


SN3 


(sn31) 


0505760501 


2007-12-31 


0:39:02.4 


40:37:48 


-1.6,-1.3 


thin 


24.26(24.26) 


medium 


28.77(28.77) 


medium 


28.78(28.78) 


S3 


(s32) 


0511380101 


2008-01-02 


0:38:52.8 


40:15:00 


-1.7,-3.3 


thin 


38.31(38.31) 


medium 


44.92(44.92) 


medium 


44.95(44.95) 


SSI 


(ssl2) 


0511380201 


2008-01-05 


0:43:28.8 


40:55:12 


-0.9,-1.4 


thin 


8.85( 8.85) 


medium 


11.28(11.28) 


medium 


11.29(11.29) 


SN2 


(sn21) 


0511380301 


2008-01-06 


0:39:40.8 


40:58:48 


-0.2, -0.4 


thin 


24.79(24.79) 


medium 


29.28(29.28) 


medium 


29.29(29.29) 


SSI 


(ssl3) 


0511380601 


2008-02-09 


0:43:28.8 


40:55:12 


-0.8,-1.8 


thin 


13.35(13.35) 


medium 


15.07(15.07) 


medium 


15.08(15.08) 



Notes: 

* : Systematic offset in RA and Dec in arcsec determined from correlations with 2MASS, USNO-Bl, LGGS and Chandra catalogues 
: All observations in full frame imaging mode 

^ : Exposure time in units of ks after screening for high background used for detection, for colour image in brackets 

* : Combination of the three observations is called b (see text), RX denotes RX J0042.6+41 15 
" : CX0M31 denotes CX0M31 J004059.2+411 



As many of the observations were affected by strong back- 
ground flares, the net exposure which can be used for our anal- 
ysis was Strongly reduced. The GTIs of the various observa- 
tions ranged over 6-56 ks, apart from observation b2 (ObsID 
0202230301) which had to be rejected, because it showed high 
background throughout the observation. The exposures for all 
three EPIC instruments are given in Cols. 8,10 and 12 of Table 2. 
The observations obtained during the summer visibility window 
of M 3 1 were afl'ected more strongly by background radiation 
than those taken during the winter window. The most afl'ected 
observations of the deep survey were reobserved. 

After screening for times of enhanced particle background, 
the second step was to examine the influence of solar wind 
charge exchange. This was done by producing soft energy 
(< 2 keV) background light curves. These lightcurves varied only 
for 10 observations, for which additional screening was neces- 
sary. The screening of enhanced background due to solar wind 
charge exchange was applied to the observations only for the 
creation of colour images, in order to avoid that these observa- 



tions will appear in the mosaic image with a tinge of red. The 
screening was not used for source detection. 

The third and last step includes the study of the back- 
ground due to detector noise. The processing chains take into 
account all known bad or hot pixels and columns and flag the 
afl'ected pixels in the event lists. We selected data with (FLAG 
& 0xfa0000)=0, excluded rows and columns near edges, and 
searched by eye for additional warm or hot pixels and columns 
in each observation. To avoid background variability over the PN 
images, we omitted the energy range from 7.2-9.2 keV where 
Strong fluorescence lines cause higher background in the outer 
detector area (Freyberg et al. 2004). 

An additional background component can occur during the 
EPIC PN ofl'set map calculation. If this period is afl'ected by high 
particle background, the ofl'set calculation will lead to a slight 
underestimate of the ofl'set in some pixels which can then result 
in blocks of pixels 4x4) with enhanced low energy signal.^ 



^ See also http://xinm2.esac.esa.int/docs/documents/ 
CAL-TN-0050- 1-0 . ps . gz 
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Figure 1. A deep optical image of M 31 (private communication: 
V. Burwitz) overplotted with the XMM-Newton fields of the survey. The 
area covered by individual EPIC observations is approximated by cir- 
cles with 14 arcmin radius. Fields observed in the "Deep XMM-Newton 
Survey of M 31" are marked with thicker lines. For presentation pur- 
poses, the ToO observation and the observations of RX J0042.6-h4115 
are omitted. 

These blocks will be found by the SAS detection tools and ap- 
pear as sources with extremely soft spectrum (so called supersoft 
sources). To reduce the number of false detections in this source 
class, we decided to include the task epreject in epchain, 
which locates the pixels with a slight underestimate of the offset 
and corrects this underestimate. To ensure that epreject pro- 
duces reliable results, difference images of the event lists ob- 
tained with and without epreject, were created. Only events 
with energies above 200 eV were used. We checked whether 
epreject removed all pixels with an enhanced low energy sig- 
nal. Only in observation nsl the difference image still shows a 
block of pixels with enhanced signal. As this block is also vis- 
ible at higher energies (PHA> 30) it cannot be corrected with 
epreject. Additionally, we ascertained that almost all pixels 
not affected during the offset map calculation have a value con- 
sistent with zero in the difference images, with two exceptions 
discussed in Sect. 5. 

3.2. Images 

For each observation, the data were split into five energy bands: 
(0.2-0.5) keV, (0.5-1. 0)keV, (1. 0-2.0) keV, (2.0-4.5) keV, and 
(4.5-1 2) keV. For the FN data, we used only single-pixel 
events (PATTERN = 0) in the first energy band, while for the 
other bands, single-pixel and double-pixel events were se- 
lected (PATTERN < 4). In the MOS cameras, single-pixel to 
quadruple-pixel events (PATTERN < 12) were used. We created 



images, background images and exposure maps (with and with- 
out vignetting correction) for PN, MOS 1 and MOS 2 in each of 
the five energy bands and masked them for the acceptable de- 
tector area. The image bin size is T'. The same procedure was 
applied in our previous M 31 and M 33 studies (PFH2005 and 
Pietsch et al. 2004). 

To create background images, the SAS task eboxdetect was 
run in local mode, in which it determines the background from 
the surrounding pixels of a sliding box, with box sizes of 5 x 5, 
10 X 10 and 20 x 20 pixels (lO^'xlO'', 20''x20''and 40''x40'0. 
The detection threshold is set to likemin= 15, which is a good 
compromise between cutting out most of the sources and leav- 
ing suflftcient area to derive the appropriate background. For the 
background calculation, a two dimensional spline is fitted to a 
rebinned and exposure corrected image (task esplinemap). The 
number of bins used for rebinning is controlled by the parameter 
nsplinenodes, which is set to 16 for all but the observations of 
the central region, where it was set to 20 (maximum value). For 
PN, the background maps contain the contribution from the "out 
of time (OoT)" events. 

3.3. Source detection 

For each observation, source detection was performed simul- 
taneously on 5 energy bands for each EPIC camera, using the 
XMM-SAS detection tasks eboxdetect and emldetect, as 
such fitting provides the most statistically robust measurements 
of the source positions by including all of the data. This method 
was also used to generate the 2XMM catalog (cf Watson et al. 
2009). In the following we describe the detection procedure 
used. 

The source detection procedure consists of two consecutive 
detection steps. An initial source list is created with the task 
eboxdetect (cf. Sect. 3.2). To select source candidates down to 
a low statistical significance level, a low likelihood threshold of 
four was used at this stage. The background was estimated from 
the previously created background images (see Sect. 3.2). 

This list is the starting point for the XMM-SAS task 
emldetect (v. 4.60.1). The emldetect task performs a 
Maximum Likelihood fit of the distribution of source counts 
(based on Cash C-statistics approach; Cash 1979), using a point 
spread function model obtained from ray tracing calculations. 
If P is the probability that a Poissonian fluctuation in the back- 
ground is detected as a spurious source, the likelihood of the 
detection is then defined as X = - In (P).^ The fit is performed 
simultaneously in all energy bands for all three cameras by sum- 
ming the likelihood contribution of each band and each camera. 
Sources exceeding the detection likelihood threshold in the full 
band (combination of the 15 bands) are regarded as detections; 
the catalogue is thus full band selected. 

The detection threshold used is 7, as in PFH2005. Some 
other parameters diflTer from the values used in PFH2005, as 
in this work a parameter setting optimised for the detection of 
extended sources was used (G. Lamer; private communication). 
The parameters in question are the event cut-out (ecut = 30.0) 
and the source selection radius (scut = 0.9) for multi-source 
fitting, the maximum number of sources into which one input 
source can be split (nmulsou= 2), and the maximum number of 
sources that can be fitted simultaneously (nmaxfit = 2). Multi- 



^ This is a simplified description as emldetect transforms the de- 
rived likelihoods to equivalent likelihoods, corresponding to the case 
of two free parameters. This allows comparison between detection runs 
with different numbers of free parameters. 
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Table 3. Count rate to energy conversion factors. The ECFs used for 
observations obtained before revolution 534 are marked with "OLD". 



Detector 


Filter 


Bl 


B2 


B3 


B4 


B5 












B ) 




EPIC PN 


thin 


11.33 


8.44 


5.97 


1.94 


0.58 




medium 


10.05 


8.19 


5.79 


1.94 


0.58 


EPIC MOS 1 


thin 


2.25 


1.94 


2.06 


0.76 


0.14 




medium 


2.07 


1.90 


2.07 


0.75 


0.15 


EPIC MOS 2 


thin 


2.29 


1.98 


2.09 


0.78 


0.15 




medium 


2.06 


1.90 


2.04 


0.75 


0.15 


EPIC MOS 1 


thin 


2.59 


2.04 


2.12 


0.76 


0.15 


OLD 


medium 


2.33 


1.98 


2.09 


0.76 


0.15 


EPIC MOS 2 


thin 


2.58 


2.04 


2.13 


0.76 


0.15 


OLD 


medium 


2.38 


1.99 


2.09 


0.75 


0.16 



PSF fitting was performed in a two stage process for objects with 
a detection likelihood larger than ten. All of the sources were 
also fitted with a convolution of a yS-model cluster brightness pro- 
file (Cavaliere & Fusco-Femiano 1976) with the XMM-Newton 
point spread function, in order to detect any possible extension 
in the detected signal. Sources which have a core radius signif- 
icantly larger than the PSF are flagged as extended. The free 
parameters of the fit were the source location, the source extent 
and the source counts in each energy band of each telescope. 

To derive the X-ray flux of a source from its measured count 
rate, one uses the so-called energy conversion factors (ECF): 



Flux = 



Rate 
ECF 



(1) 



These factors were calculated using the detector response, and 
depended on the used filter, the energy band in question, and 
the spectrum of the source. As we wanted to apply the conver- 
sion factors to all sources found in the survey, we assumed a 
power law model with photon index F = 1.7 and the Galactic 
foreground absorption of A^h = 7 x 10^^ cm"^ (Stark et al. 1992, 
see also PFH2005) to be the universal source spectrum for the 
ECF calculation. 

The ECFs (see Table 3) were derived with XSPEC4(v 11.3.2) 
using response matrices (V.7.1) available from the XMM- 
Newton calibration homepage^. As all necessary corrections 
of the source parameters {e. g. vignetting corrections) were in- 
cluded in the image creation and source detection procedure^, 
the on axis ECF values were derived (cf. Watson et al. 2009). 
The fluxes determined with the ECFs given in Table 3 are ab- 
sorbed (/. e. observed) fluxes and hence correspond to the ob- 
served count rates, which are derived in the eml detect task. 

During the mission lifetime, the MOS energy distribution 
behaviour has changed. Near the nominal boresight positions, 
where most of the detected photons hit the detectors, there has 
been a decrease in the low energy response of the MOS cameras 
(Read et al. 2006). To take this efl'ect into account, difl'erent re- 
sponse matrices for observations obtained before and after the 
year 2005 were used (see Table 3). 

For most sources, band 5 just adds noise to the total count 
rate. If converted to flux, this noise often dominates the total flux 
due to the small ECF. To avoid this problem we calculated count 



http : //heasarc . gs£c . gov/docs/xanadu/xspec 
^ http : //xvmZ . esac . esa . int/external/xitim_sw_cal/ 
calib/epic_£iles . shtml 
^ especially in the emldetect task 



rates and fluxes for detected sources in the "XID" (0.2-4.5) keV 
band (bands 1 to 4 combined). While for most sources this is a 
good solution, for extremely hard or soft sources there may still 
be bands just adding noise. This, then, may lead to rate and flux 
errors that seem to falsely indicate a lower source significance. A 
similar eff'ect occurs in the combined rates and fluxes, if a source 
is detected primarily by one instrument (e. g. soft sources in PN). 

Sources are entered in the XMM LP-total catalogue from the 
observation in which the highest source detection likelihood is 
obtained (either combined or single observations). For variable 
sources this means that the source properties given in the XMM 
LP-total catalogue (see Sect. 5 and Table 5) are those observed 
during their brightest state. 

We rejected spurious detections in the vicinity of bright 
sources. In regions with a highly structured background, the SAS 
detection task emldetect registered some extended sources. We 
also rejected these "sources" as spurious detections. In an addi- 
tional step we checked whether an object had visible contours in 
at least one image out of the five energy bands. The point-like or 
extended nature, which was determined with emldetect, was 
taken into account. In this way, "sources" that are fluctuations in 
the background, but which were not fully modelled in the back- 
ground images, were detected. In addition, objects located on 
hot pixels, or bright pixels at the rim or in the corners of the in- 
dividual CCD chips (which were missed during the background 
screening) were recognised and excluded from the source cata- 
logue, especially if they were detected with a likelihood larger 
than six in one detector only. 

To allow for a statistical analysis, the source catalogue only 
contains sources detected by the SAS tasks eboxdetect and 
emldetect as described above, /. e. the few sources that were 
not detected by the analysis program, despite being visible on 
the X-ray images, have not been added by hand as it was done 
in previous studies (SPH2008; PFH2005). 

To classify the source spectra, we computed four hardness 
ratios. The hardness ratios and errors are defined as: 



HR/: 



i+l 



Bi+\ + Bt 



^^ and EHR/ = 2V(^^^^ 



(Bi^i + Bt)^ 



for / = 1 to 4, where Bi and EBi denote count rates and corre- 
sponding errors in energy band /. 

3.4. Astrometrical corrections 

To obtain astrometrically-corrected positions for the sources of 
the five central fields we used the SAS-task eposcorr with 
Chandra source lists (Kong et al. 2002b; Kaaret 2002; Williams 
et al. 2004). For the other fields we selected sources from the 
USNO-Bl (Monet et al. 2003), 2MASS (Skrutskie et al. 2006) 
and Local Group Galaxy Survey (LGGS; Massey et al. 2006) 
catalogues^. 

3.4.1 . Astrometry of optical/infrared catalogues 

In a first step, we examined the agreement between the positions 
given by the various optical catalogues.^ A close examination 
of the shifts obtained, showed significant diff'erences between 



For the remainder of the subsection we will call all three catalogues 
"optical catalogues" for easier readability, although the 2MASS cata- 
logue is an infrared catalogue. 

^ From the LGGS catalogue only sources brighter than 21 mag were 
used in order to be comparable to the brightness limit of the USNO-B 1 
catalogue. 
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the positions given in the individual catalogues. In summary, be- 
tween the USNO-B 1 and LGGS catalogues we found an offset 
of: -a/ 197 in R.A. and (7/067 in Dec^; and between the USNO- 
Bl and 2MASS catalogues we found an offset of: -(7/108 in 
R.A. and a/ 204 in Dec. We chose the USNO-B 1 catalogue as a 
reference, since it covers the entire field observed in the Deep 
XMM-Newton survey, and in addition it provides values for the 
proper motion of the optical sources. 

Since the optical catalogues, as well as the Deep XMM- 
Newton catalogue, are composed of individual observations of 
sub-fields of M 3 1 , we searched for systematic drifts in the posi- 
tional zero points from region to region. However no systematic 
offsets were found. 

Finally, we applied the corrections found to the sources in 
the LGGS and 2MASS catalogues, to bring all catalogues to the 
USNO-B 1 reference frame. 

The offsets found between the USNO-B 1 and 2MASS cata- 
logues can be explained by the independent determination of the 
astrometric solutions for these catalogues. Given that the posi- 
tions provided in the LGGS catalogue are corrected with respect 
to the USNO-B 1 catalogue (see Massey et al. 2006), the offset 
found in right ascension was totally unexpected and cannot be 
explained. 

3.4.2. Corrections of the X-ray observations 

From the positionally corrected catalogues, we selected sources 
which either correlate with globular clusters from the Revised 
Bologna Catalogue (V.3.4, January 2008; Galleti et al. 2004, 
2005, 2006, 2007) or with foreground stars, characterised 
by their optical to X-ray flux ratio (Maccacaro et al. 1988) 
and their hardness ratio (see source selection criteria given in 
Table 6 and Stiele et al. 2008). For sources selected from the 
USNO-B 1 catalogue, we used the proper motion corrected 
positions. We then used the SAS-task eposcorr to derive the 
offset of the X-ray aspect solution. Four observations did not 
have enough optical counterparts to apply this method. The lack 
of counterparts is due to the very short exposure times resulting 
after the screening for high background (obs. s3, ssl2, ssl3) 
and the location of the observation (obs. snll). In these cases, 
we used bright persistent X-ray sources, which we correlated 
with another observation of the same field. We checked for 
any residual systematic uncertainty in the source positions and 
found it to be well characterised by a conservative Icr value of 
0V5. This uncertainty is due to positional errors of the optical 
sources as well as inaccuracy in the process of the determination 
of the offset between optical and X-ray sources, and is called 
systematic positional error. The appropriate offset, given in 
Col. 6 of Table 2, was applied to the event file of each pointing, 
and images and exposure maps were then reproduced with the 
corrected astrometry. 

Fields that were observed at least twice are treated in a spe- 
cial way, which is described in the following section. 

3.5. Multiple observations of selected fields 

The fields that were observed more than once were the central 
field, the fields pointing on RX 10042.6-^4115^^, two fields lo- 
cated on the major axis of M 31 (S2, N2) and all fields of the 
"Large Survey" located in the southern part of the galaxy (SSI, 

^ the offset in declination is negligible 
The combination of observations bl, b3 and b4 is called b. 



SS2, SS3, S3, SN3, SN2, SNl). To reach higher detection sensi- 
tivity we merged the images, background images and exposure 
maps of observations which have the same pointing direction 
and were obtained with the same filter setting. Subsequently, 
source detection, as described in Sect. 3.3, was repeated on the 
merged data. For the S2 field, there are two observations with 
different filter settings. In this case, source detection was per- 
formed simultaneously on all 15 bands of both observations, /. e. 
on 30 bands simultaneously. The N2 field was treated in the same 
way. For the central field images, background images and expo- 
sure maps of observations cl, c2 and c3 were merged. These 
merged data were used together with the data of observation c4 
to search for sources simultaneously; in this way it was possible 
to take into account the different ECFs for the different filters. 
One field was observed twice with slightly different pointing di- 
rection in observations snl and snll; simultaneous source de- 
tection was used for these observations also. 

3.6. Variability calculation 

To examine the time variability of each source listed in the total 
source catalogue, we determined the XID flux at the source po- 
sition in each observation or at least an upper limit for the XID 
flux. We used the task eml detect with fixed source positions 
when calculating the total flux. To get fluxes and upper limits 
for all sources in the input list we set the detection likelihood 
threshold to 0. 

A starting list was created from the full source catalogue, 
which only contains the identification number and position of 
each source located in the field examined. To give correct re- 
sults, the task emldetect has to process the sources from the 
brightest one to the faintest one. We, therefore, had to first or- 
der the sources in each observation by the detection likelihood. 
For sources not visible in the observation in question we set the 
detection likelihood to 0. This list was used as input for a first 
emldetect run. In this way we achieved an output list in which 
a detection likelihood was allocated to every source. For a final 
examination of the sources in order of detection likelihood, a 
second emldetect run was necessary. 

We only accepted XID fluxes for detections > 3 cr; otherwise 
we used a 3 cr upper limit. To compare the XID fluxes between 
the different observations, we calculated the significance of the 
difference 

^ _ -^max ~~ -^min 
y CTmax + CTmin 

and the ratio of the XID fluxes V = i^max/^min, where Fmax 
and Fmin are the maximum and minimum (or upper limit) source 
XID flux, and (Xmax and (Tmin are the errors of the maximum and 
minimum flux, respectively. This calculation was not performed 
whenever Fmax was an upper limit. Finally, the largest XID flux 
of each source was derived, excluding upper limits. 

3.7. Spectral analysis 

To extract the X-ray spectrum of individual sources, we se- 
lected an extraction region and a corresponding background re- 
gion which was at least as large as the source region, was lo- 
cated on the same CCD at a similar off axis angle as the source, 
and did not contain any point sources or extended emission. For 
EPIC PN, we only accepted single-pixel events for the spec- 
tra of supersoft sources, while for all other spectra single and 
double-pixel events were used. For the EPIC-MOS detectors. 
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single-pixel through to quadruple-pixel events were always used. 
Additionally, we only kept events with FLAG = for all three de- 
tectors. For each extraction region, we produced the correspond- 
ing response matrix files and ancillary response files. 

For each source, the spectral fit was obtained by fitting all 
three EPIC spectra simultaneously, using the tool XSPEC. For 
the absorption, we used the TBabs model, with abundances from 
Wilms et al. (2000) and photoelectric absorption cross-sections 
from Balucinska-Church & McCammon (1992) with a new He 
cross-section based on Yan et al. (1998). 

3.8. Cross correlations 

Sources were regarded as correlating if their positions over- 
lapped within their 3cr (99.73%) positional errors, defined as 
(Watson et al. 2009): 

Apos < 3.44 X ^0-2 + ^syst + 3 X 0-ccat (4) 

where (Xstat is the statistical and (Xsyst the systematic error of the 
X-ray sources detected in the present study. The statistical error 
was derived by emldetect. The determination of the system- 
atic error is described in Sect. 3.4. We use a value of 0V5, for 
all sources. The positional error of the sources in the catalogue 
used for cross-correlation is given by (Xccat- The values of (Xccat 
(68% error) used for the diff'erent X-ray catalogues can be found 
in Table 4. Exceptions to Eq. 4 are sources that are listed in more 
than one catalogue or that are resolved into multiple sources with 
Chandra. The first case is restricted to catalogues with compa- 
rable spatial resolution and hence positional uncertainty. 

To identify the X-ray sources in the field of M 31 we 
searched for correlations with catalogues in other wavelength 
regimes. The XMM-Newton source catalogue was correlated 
with the following catalogues and public data bases: 

Globular Clusters: Bologna Catalogue (V.3.5, March 2008; 

Galleti et al. 2004, 2005, 2006, 2007, 2009, cr,,^t = 0.''2; 

REV 3.5), Caldwell et al. (2009, cTccat = 0.''2), Hodge et al. 

(2009, o-ccat = 0.''5), Krienke & Hodge (2008, cr,,^t = 0.''2), 

Krienke & Hodge (2007, cr^cat = 0.''2), Fan et al. (2005), 

Magnier(1993,crccat = I'O 
Novae: Nova list of the M 31 Nova Monitoring Project^ ^ (cTccat 

is given for each individual source), PHS2007, Pietsch 

(2010) 

Supernova Remnants: Dodorico et al. (1980, 19 srcs), 
Walterbos & Braun (1992, 967 srcs) and Braun & Walterbos 
(1993, 58 srcs), Magnier et al. (1995, 233 srcs); An X-ray 
source is considered as correlating with a SNR, if the X-ray 
source position (including 3cr error) lies within the extent 
given for the SNR. 

Radio Catalogues: Gelfand et al. (2005, (Xccat is given for each 
individual source), Gelfand et al. (2004, (Xccat is given for 
each individual source), Kimball & Ivezic (2008, (Xccat = 
3''), Braun (1990, (Xccat is given for each individual source), 
NVSS (NRAO/VLA Sky Survey^^. Condon et al. 1998, cTccat 
is given for each individual source) 

H II Regions, H a Catalogue: Walterbos & Braun (1992, (Xccat 
is given for each individual source), Massey et al. (2007, 

C^ccat = 0.''2) 



^ ^ http : //www . mpe . mpg . de/ ~ m3 lnovae/opt/m3 1/M3 l_table . 
html 

http : //www . cv . nrao . edu/nvss/NVSSlist . shtml 



Table 4. X-ray source catalogues used for cross-correlation and the used 
positional errors 



X-ray catalogue* crl^^^ X-ray catalogue* cr*^ 



PFH2005 




DKG2004 


or3 


SPH2008 




WNG2006 


or3 


SHP97 




VG2007 


or4 


SHL2001 




OBT2001 




PFJ93 




O2006 


\" 


TF91 




SBK2009 


3"'^ 


Ka2002 


073 


D2002 


075 


KGP2002 




TP2004 


\" 


WGK2004 




ONB2010 


\" 



Notes: 

* : ^ indicates that the catalogue provides cTccat values for each source 
individually 

^ : value taken from indicated paper 

* : TF91: Trinchieri & Fabbiano (1991), PFJ93: Primini et al. (1993), 
SHP97: Supper et al. (1997), SHL2001: Supper et al. (2001), OBT2001: 
Osborne et al. (2001), D2002: Di Stefano et al. (2002), KGP2002: 
Kong et al. (2002b), Ka2002: Kaaret (2002), WGK2004: Williams et al. 
(2004), DKG2004: Di Stefano et al. (2004), TP2004: Trudolyubov & 
Priedhorsky (2004), PFH2005: Pietsch et al. (2005b), O2006: Orio 

(2006) , WNG2006: Williams et al. (2006b), VG2007: Voss & Gilfanov 

(2007) , SPH2008: Stiele et al. (2008), SBK2009: Shaw Greening et al. 
(2009), ONB2010: Orio et al. (2010) 

Optical Catalogues: USNO-B 1 (Monet et al. 2003, cTccat is given 

for each individual source). Local Group Survey (LSG; 

Massey et al. 2006, cTccat = 0.''2) 
Infrared catalogues: 2MASS (Skrutskie et al. 2006, (Xccat is 

given for each individual source). Mould et al. (2008, (Xccat = 

0.''8, for Table 2: cTccat = 0.''5) 
Data bases: the SIMBAD catalogue (Centre de Donnees as- 

tronomiques de Strasbourg; hereafter SIMBAD) , the NASA 

Extragalactic Database (hereafter NED) 



4. Colour image 

Figure 2 shows the combined, exposure corrected EPIC PN, 
MOS 1 and MOS 2 RGB (red-green-blue) mosaic image of the 
Deep Survey and archival data. The colours represent the X-ray 
energies as follows: red: 0.2-l.OkeV, green: 1.0-2.0keV and 
blue: 2.0-1 2 keV. The optical extent of M 31 is indicated by the 
D25 ellipse and the boundary of the observed field is given by 
the green contour. The image is smoothed with a 2D-Gaussian 
of 20'' FWHM. In some observations, individual noisy MOS 1 
and MOS 2 CCDs are omitted. The images have not been cor- 
rected for the background of the detector or for vignetting. 
The colour of the sources reflects their class. Supersoft sources 
appear in red. Thermal SNRs and foreground stars are orange to 
yellow. "Hard" sources (background objects, mainly AGN, and 
X-ray binaries or Crab-like SNRs) are blue to white. 

Logarithmically scaled XMM-Newton EPIC low back- 
ground images made up of the combined images from the PN, 
MOS 1 and MOS 2 cameras in the (0.2-4.5) keV XID band for 
each M 31 observation can be found in the Appendix. The im- 
ages also show X-ray contours, and the sources from the XMM 
LP-total catalogue are marked with boxes. 

http : //simbad . u- strasbg . f r/simbad 
http : //nedwww . ipac . caltech . edu 
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Figure 2. Combined EPIC PN, MOS 1 and 
MOS 2 RGB image of the Deep M 31 Survey 
including archival data. The optical extent of 
M 31 is indicated by the D25 ellipse and the 
boundary of the observed field is given by the 
green contour. The central region, marked with 
the yellow square, is shown in higher resolution 
in the upper right corner. For more details see 
Sect. 4. 



5. Source catalogue (XMM LP-total) 

The source catalogue of the Deep XMM-Newton survey of 
M 31 (hereafter XMM LP-total catalogue) contains 1 897 X-ray 
sources. Of these sources 914 are detected for the first time in 
X-rays. 

The source parameters are summarised in Table 5, which 
gives the source number (Col. 1), detection field from which 
the source was entered into the catalogue (2), source position 
(3 to 9) with 3o- (99.73%) uncertainty radius (10), likelihood 
of existence (11), integrated PN, MOS 1 and MOS 2 count rate 
and error (12,13) and flux and error (14,15) in the (0.2-4.5) keV 
XID band, and hardness ratios and errors (16-23). Hardness ra- 
tios are calculated only for sources for which at least one of the 
two band count rates has a significance greater than 2cr. Errors 
are the properly combined statistical errors in each band and can 
extend beyond the range of allowed values of hardness ratios as 
defined previously (-1.0 to 1.0; Eq. 2). The "Val" parameter (Col 
24) indicates whether the source is within the field of view (true 
or false, "T" or "F") in the PN, MOS 1 and MOS 2 detectors 
respectively. 

Table 5 also gives the exposure time (25), source existence 
likelihood (26), the count rate and error (27, 28) and the flux and 



error (29, 30) in the (0.2-4.5) keV XID band, and hardness ra- 
tios and errors (31-38) for the EPIC PN. Columns 39 to 52 and 
53 to 66 give the same information corresponding to Cols. 25 to 
38, but for the EPIC MOS 1 and MOS 2 instruments. Hardness 
ratios for the individual instruments were again screened as de- 
scribed above. From the comparison between the hardness ra- 
tios derived from the integrated PN, MOS 1 and MOS 2 count 
rates (Cols. 16-23) and the hardness ratios from the individual 
instruments (Cols. 31-38, 45-52 and 59-66), it is clear that the 
combined count rates from all instruments yielded a significantly 
larger fraction of hardness ratios above the chosen significance 
threshold. 

Column 67 shows cross correlations with published M 3 1 X- 
ray catalogues (c/ Sect. 3.8). We discuss the results of the cross 
correlations in Sects. 9 and 10. 

In the remaining columns of Table 5, we give information 
extracted from the USNO-Bl, 2MASS and LGGS catalogues 
(c/ Sect. 3.8). The information from the USNO-Bl catalogue 
(name, number of objects within search area, distance, B2, R2 
and I magnitude of the brightest^ ^ object) is given in Cols. 68 to 
73. The 2MASS source name, number of objects within search 
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area, and the distance can be found in Cols. 74 to 76. Similar 
information from the LGGS catalogue is given in Cols. 77 to 82 
(name, number of objects within search area, distance, V magni- 
tude, V-R and B-V colours of the brightest^^ object). To improve 
the reliability of source classifications we used the USNO-Bl 
B2 and R2 magnitudes to calculate 



log I — I = log (/x) + ^ , + 5.37, 



opt ) 



2x2.5 



and the LGGS V magnitude to calculate 
logfA] = log(/.)+?^+5.37, 



opt ) 



2.5 



(5) 



(6) 



following Maccacaro et al. (1988, see Cols. 83-86). 

The X-ray sources in the XMM LP-total catalogue are identi- 
fied or classified based on properties in X-rays (HRs, variability, 
extent) and of the correlated objects in other wavelength regimes 
(Cols. 87 and 88 in Table 5). For classified sources the class 
name is given in angled brackets. Identification and classifica- 
tion criteria are summarised in Table 6, which provides, for each 
source class (Col. 1), the classification criteria (2), and the num- 
bers of identified (3) and classified (4) sources. The hardness 
ratio criteria are based on model spectra. Details on the defini- 
tion of these criteria can be found in Sect. 6 of PFH2005. As 
we have no clear hardness ratio criteria to discriminate between 
XRBs, Crab-like supernova remnants (SNRs) or AGN we intro- 
duced a <hard> class for those sources. If such a source shows 
strong variability (i.e. V>10) on the examined time scales it 
is likely to be an XRB. Compared with SPH2008 the HR2 se- 
lection criterion for SNRs was tightened (from HR2< -0.2 to 
HR2-hEHR2< -0.2) to exclude questionable SNR candidates 
from the class of SNRs. If we applied the former criterion to 
the survey data, -35 sources would be classified as SNRs in ad- 
dition to those listed in Table 6. Most of the 35 sources are lo- 
cated outside the D25 ellipse, and none of them correlates with 
an optically identified SNR, a radio source, or an Hll region. In 
addition, the errors in HR2 are of the same order as the HR2 
values. It is therefore very likely that these sources do belong to 
other classes, since the strip between -0.3 <HR2<0 is populated 
by foreground stars, XRBs, background objects, and candidates 
for these three classes. Outcomes of the identification and clas- 
sification processes are discussed in detail in Sects. 9 and 10. 

The last column (89) of Table 5 contains the XMM-Newton 
source name as registered to the lAU Registry. Source names 
consist of the acronym XMMM31 and the source position as 
follows: XMMM31 Jhhmmss.s-i-ddmmss, where the right ascen- 
sion is given in hours (hh), minutes (mm) and seconds (ss.s) 
truncated to decimal seconds and the declination is given in de- 
grees (dd), arc minutes (mm) and arc seconds (ss) truncated to 
arc seconds, for equinox 2000. In the following, we refer to indi- 
vidual sources by their source number (Col. 1 of Table 5), which 
is marked with a "N-" at the front of the number. 

Of the 1 897 sources, 1 247 can only be classified as <hard> 
sources, while 123 sources remain without classification. Two 
of them (N^ 482, N^ 768) are highly aff'ected by optical loading; 
both "X-ray sources" coincide spatially with very bright opti- 
cal foreground stars (USNO-Bl R2 magnitudes of 6.76 and 6.74 
respectively). The spectrum of source N- 482 is dominated by 
optical loading. This becomes evident from the hardness ratios 
which indicate an SSS. For N- 768 the hardness ratios would al- 
low a foreground star classification. The obtained count rates and 
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Figures. Distribution of the source fluxes in the 0.2^.5 keV (XID) 
band. The diagrams show the number of sources at each flux bin, plotted 
versus the flux, using logarithmic scales. The inlay shows the number 
of sources for XID fluxes smaller than 5xlO~^^ erg cm~^ s~\ on linear 
scales. The blue histogram gives the distribution of sources classified or 
identified as either SSSs, SNRs, XRBs or GlCs. 



fluxes of both sources are aff'ected by the usage of epreject, 
which neutralises the corrections applied for optical loading. 
Therefore residuals are visible in the difl'erence images created 
from event lists obtained with and without epreject. As we 
cannot exclude the possibility that some of the detected photons 
are true X-rays - especially for source N- 768 -, we decided 
to include them in the XMM LP-total catalogue, but without a 
classification. 

5A. Flux distribution 

The faintest source (N^ 526) has an XID band flux of 5.8x10"^^ 
erg cm"^ s"^ The source with the highest XID Flux (N- 966, 
XID band flux of 3.75x10"^^ erg cm"^ s"^) is located in the 
centre of M 31 and identified as a Z- source LMXB (Barnard 
et al. 2003). This source has a mean absorbed XID luminosity of 
2.74x10^^ ergs-i. 

Figure 3 shows the distribution of the XID (0.2-4.5 keV) 
source fluxes. Plotted are the number of sources in a certain flux 
bin. We see from the inlay that the number of sources starts to 
decrease in the bin from 2.4 to 2.6x10"^^ erg cm"^ s"^ This XID 
flux roughly determines the completeness limit of the survey and 
corresponds to an absorbed 0.2-4.5 keV limiting luminosity of 
-2x10^^ erg s"^ 

Previous X-ray studies (Williams et al. 2004, and references 
therein) noted a lack of bright sources (Lx ^10^^ erg s"^; 0.1- 
10 keV) in the northern half of the disc compared to the southern 
half. This finding is not supported in the present study. Excluding 
the pointings to the centre of M 3 1 , we found in the remaining 
observations 13 sources in each hemisphere that were brighter 
than Lxabs ^10^^ erg s"^^^ The reason our survey does not sup- 



in B magnitude 



The luminosity is based on XID Fluxes. Using the total 0.2-12 keV 
band the result does not change (23 in the northern half and 24 in the 
southern half). 
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Table 6. Summary of identifications and classifications. 



Source class 


Selection criteria 


identified 


classified 


fg Star 


log( ) < - 1 .0 and HR2-EHR2 < 0.3 and HR3-EHR3 < -0.4 or not defined 


40 


223 


AGN 


Radio source and not classification as SNR from HR2 or optical/radio 


11 


49 


Gal 


optical id with galaxy 


4 


19 


GCl 


X-ray extent and/or spectrum 


1 


5 


sss 


HRl <0.0, HR2-EHR2< -0.96 or HR2 not defined, HR3, HR4 not defined 




30 


SNR 


HRl > -0. 1 and HR2+EHR2 < -0.2 and not a fg Star, or id with optical/radio SNR 


25 


31 


GIC 


optical id 


36 


16 


XRB 


optical id or X-ray variability 


10 


26 


hard 


HR2-EHR2> -0.2 or only HR3 and/or HR4 defined, and no other classification 




1247 




. XMM-Newton Deep Survey image over plotted with sources 
that have an absorbed 0.2-4.5 keV luminosity larger than 10^^ erg s~^ 
Striking features are the two patches located north and south of the cen- 
tre. The central region (same as in Fig. 2) is shown with higher resolu- 
tion in the upper right corner. 




100000 200000 300000 



Figure 5. Exposure map of all fields of the XMM LP-total catalogue. 
For details see Sect. 5.2. 



port the old results is that we found several bright sources in the 
outer regions of the northern half of the disk, which have not 
been covered in Williams et al. (2004, and references therein). 
In the central field of M 31, a total of 41 sources brighter than 
Lx^lO^^ erg s'^ (0.2-4.5 keV) were found. 

Figure 4 shows the spatial distribution of the bright sources. 
Striking features are the two patches located north and south of 
the centre. The southern one seems to point roughly in the di- 
rection of M 32 (N- 995), while the northern one ends in the 
globular cluster B 1 16 (N- 947). However there is no association 
to any known spatial structure of M 31, like e. g. the spiral arms. 



5.2. Exposure map 

Figure 5 shows the exposure map used to create the colour im- 
age of all XMM-Newton Large Survey and archival observations 
(Fig. 2). The combined MOS exposure was weighted by a factor 
of 0.4, before being added to the PN exposure. However, this 
map does not quite represent the exposures used in source de- 
tection; overlapping regions were not combined during source 
detection. 

From Fig. 5 we see that the exposure for most of the sur- 
veyed area is rather homogeneous. Exceptions are the central 
area, overlapping regions and observation h4. 
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5.3. Hardness ratio diagrams 

We plot X-ray colour/colour diagrams based on the HRs (see 
Fig. 6). Sources are plotted as dots if the error in both contribut- 
ing HRs is below 0.2. Classified and identified sources are plot- 
ted as symbols in all cases. Symbols including a dot therefore 
mark the well-defined HRs of a class. 

From the HR1-HR2 diagram (upper panel in Fig. 6) we note 
that the class of SSSs is the only one that can be defined based on 
hardness ratios alone. In the part of the HR1-HR2 diagram that 
is populated by SNRs, most of the foreground stars and some 
background objects and XRBs are also found. 

Foreground star candidates can be selected from the HR2- 
HR3 diagram (middle panel in Fig. 6), where most of them are 
located in the lower left corner. The HR3-HR4 diagram (lower 
panel in Fig. 6) does not help to disentangle the diff'erent source 
classes. Thus, we need additional information from correlations 
with sources in other wavelengths or on the source variability or 
extent to be able to classify the sources. 

5.4. Extended sources 

The XMM LP-total catalogue contains 12 sources which are fit- 
ted as extended sources with a likelihood of extension larger than 
15. This value was chosen so as to minimise the number of spu- 
rious detections of extended sources (H. Brunner; private com- 
munication), as well as keeping all sources that can clearly be 
seen as extended sources in the X-ray images. A convolution of a 
yS-model cluster brightness profile (Cavaliere & Fusco-Femiano 
1976) with the XMM-Newton point spread function was used 
to determine the extent of the sources (cf. Sect. 3.3). This model 
describes the brightness profile of galaxy clusters, as 



f(x,y) 



1 



■ xof + (j - yof 



-3/2 



(7) 



where denotes the core radius; this is also the extent parameter 
given by emldetect. 

Table 7 gives the source number (Col. 1), likelihood of de- 
tection (2), the extent found (3) and its associated error (4) in 
arcsec, the likelihood of extension (5), and the classification of 
the source (6, see Sect. 9.2) for each of the 12 extended sources. 
Additional comments taken from Table 5 are provided in the last 
column. 

The extent parameter found for the sources ranges from 6!2 
to 23r03 (see Fig. 7). The brightest source (N^ 1 795), which has 
the highest likelihood of extension and the second largest extent, 
was identified from its X-ray properties as a galaxy cluster lo- 
cated behind M 31 (Kotov et al. 2006). The iron emission lines 
in the X-ray spectrum yield a cluster redshift of z = 0.29. For 
further discussion see Sect. 9.2. 



6. Variability between XMM-Newton observations 

To examine the long-term time variability of each source, we 
determined the XID flux at the source position in each observa- 
tion or at least an upper limit for the XID flux. The XID fluxes 
were used to derive the variability factor and the significance of 
variability (cf. Sect. 3.6). 

The sources are taken from the XMM LP-total catalogue 
(Table 5). Table 8 contains all information necessary to exam- 
ine time variability. Sources are only included in the table if they 
are observed at least twice. Column 1 gives the source number. 
Columns 2 and 3 contain the flux and the corresponding error 
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Figure 6. Hardness ratios of sources detected by XMM-Newton EPIC. 
Sources with HR errors smaller then 0.20 on both HR(/) and HR(/ + 1) 
are shown as dots. Foreground stars and candidates are marked as big 
and small stars, AGN and candidates as big and small crosses, back- 
ground galaxies and galaxy clusters as big "X" and their candidates as 
small "X", SSS candidates as triangles, SNRs and candidates as big and 
small octagons, GlCs and XRBs as big squares and their candidates as 
small squares. 



H. Stiele et al: The deep XMM-Newton Survey of M 31 



13 



Table 7. Extended sources in the XMM LP-total catalogue 



SRC 


DET_ML 


EXT+ 


EEXT+ 


EXT_ML 


XFLUX* 


XEFLUX* 


class 


comment^ 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


141 


65.08 


11.22 


1.29 


23.68 


1.45 


0.20 


<GC1> 


GLG127(Gal), 37 W 025 A (IR, RadioS; NED) 


199 


275.16 


17.33 


1.05 


174.73 


4.31 


0.29 


<hard> 




252 


222.05 


14.64 


1.12 


81.60 


4.40 


0.49 


<GC1> 


5 optical objects in error box 


304 


299.75 


15.10 


0.92 


133.62 


2.20 


0.18 


<GC1> 


B242 [CHM09]; RBC3.5: <G1C> 


442 


33.76 


11.60 


1.71 


15.44 


1.62 


0.28 


<hard> 




618 


271.08 


6.20 


0.73 


42.86 


3.15 


0.21 


<hard> 




718 


77.75 


7.18 


1.23 


21.47 


0.58 


0.07 


Gal 


B052 [CHM09], RBC3.5 


1 130 


168.31 


10.80 


0.97 


44.23 


3.27 


0.31 


<hard> 




1543 


70.49 


11.87 


1.37 


28.63 


1.51 


0.19 


<GC1> 


[MLA93] 1076 PN (SIM,NED) 


1795 


11416.36 


18.79 


0.29 


4 169.74 


98.87 


1.43 


GCl 


GLG253 (Gal), [B90] 473, z=0.3 [KTV2006] 


1859 


107.09 


13.73 


1.40 


43.89 


1.23 


0.19 


<hard> 




1912 


332.06 


23.03 


1.23 


213.90 


5.43 


0.37 


<GC1> 


cluster of galaxies candidate 



Notes: 

+ : Extent and error of extent in units of I" ; I" corresponds to 3.8 pc at the assumed distance of M 31 
* : XID Flux and flux error in units of 1x10"^"^ erg cm"^ s"^ 
^ : Taken from Table 5 




Figure 7. Distribution of extent parameter. 



in the (0.2-4.5) keV XID band. The hardness ratios and errors 
are given in columns 4 to 11. Column 12 gives the type of the 
source. All this information was taken from Table 5. 

The subsequent 140 columns provide information related to 
individual observations in which the position of the source was 
observed. Column 13 gives the name of one of these observa- 
tions, which we will call observation 1. The EPIC instruments 
contributing to the source detection in observation 1, are indi- 
cated by three characters in the "obsl_val" parameter (Col. 14, 
first character for PN, second MOS 1, third MOS 2), each one 
being either a "T" if the source is inside the FoV, or "F" if 
it lies outside the FoV. Then the count rate and error (15,16) 
and flux and error (17,18) in the (0.2-4.5) keV XID band, and 
hardness ratios and error (19-26) of observation 1 are given. 
Corresponding information is given for the remaining observa- 
tions which cover the position of the source: obs. 2 (cols. 27- 
40), obs. 3 (41-54), obs. 4 (55-68), obs. 5 (69-82), obs. 6 
(83-96), obs. 7 (97-110), obs. 8 (111-124), obs. 9 (125-138), 
obs. 10 (13 9- 15 2). Whether the columns corresponding to obs. 3 
- obs. 10 are filled in or not, depends on the number of observa- 
tions in which the source has been covered in the combined EPIC 
FoV. This number is indicated in column 153. The maximum 
significance of variation and the maximum flux ratio (fvar_max) 
are given in columns 154 and 155. As described in Sect. 3.6, only 
detections with a significance greater than 3cr were used, other- 
wise the 3(T upper limit was adopted. Column 156 indicates the 
number of observations that provide only an upper limit. The 
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Figures. Variability factor of sources from the XMM LP-total cat- 
alogue in the 0.2-4.5 keV band derived from average fluxes of the 
XMM-Newton EPIC observations plotted versus maximum detected 
flux (erg cm~^ s"^. Source classification is indicated: Foreground stars 
and candidates are marked as big and small stars, AGN and candidates 
as big and small crosses, background galaxies and galaxy clusters as 
big "X" and their candidates as small "X", SSS candidates as triangles, 
SNRs and candidates as big and small octagons, GlCs and XRBs as big 
squares and their candidates as small squares. Sources with a statistical 
significance for the variability below 3 are marked in green. 



maximum flux (fmax) and its error are given in columns 157 
and 158. 

In a few cases a maximum flux value could not be derived, 
because each observation only yielded an upper limit. There can 
be two reasons for this: The first reason is that faint sources de- 
tected in merged observations may not be detected in the indi- 
vidual observations at the 3cr limit. The second reason is that in 
cases where the significance of detection was not much above 
the 3cr limit, it can become smaller than the 3cr limit when the 
source position is fixed to the adopted final mean value from all 
observations. 
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Table 10. Sources with maximum flux larger than 8x10"^^ erg cm"^ a statistical significance of variability larger than 10 and a flux variability smaller than 5, ordered by flux. 



Source 


fvar 


svar 


fmax* 


efmax* 


class"^ 


Comment^ 


966 


1.63 


49.01 


46.73 


0.59 


XRB 


l(sv,z), 2, 10(v), 12(v), 13, 14, 20, 22(v), 25(LMXB), 27, 28(1.56) 


877 


3.13 


49.13 


16.06 


0.20 


<hard> 


l(sv), 2, 10(v), 12(v), 13, 14, 20(v), 22(v), 27, 28(3.05) 


745 


2.43 


26.89 


12.65 


0.18 


AGN 


13, 14 


1 157 


1.32 


11.10 


9.87 


0.25 


GIC 


l(sv), 2, 5, 10, 12, 13, 14, 20, 21, 22(v), 27, 28(1.37) 


1060 


2.13 


30.00 


9.04 


0.14 


<XRB> 


l(sv), 2, 10, 12, 13, 14, 20(v, NS-LMXB), 22(v), 27 


1 171 


4.14 


18.86 


9.02 


0.41 


QIC 


l(d,sv), 2(t, 53.4), 5, 10, 12, 13, 14, 16, 20, 22, 27, 28(2.47) 


1 116 


3.76 


51.98 


8.16 


0.10 


GIC 


l(sv), 2(t, 58.6), 3(t, 33), 5, 10, 12, 13, 14, 16, 20, 21, 22(v,t), 27 



Notes: 

* : maximum XID flux and error in units of 1x10"^^ erg cm~^ s~^ or maximum absorbed 0.2-4.5 keV luminosity and error in units of 7.3x10^^ erg s"^ 

: class according to Table 6 
^ : for comment column see Table 9 
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Figure 9. Variability factor of sources from the XMM LP-total catalogue in the 0.2-4.5 keV band (derived from the average fluxes of the XMM- 
Newton EPIC observations) plotted versus HRl in the left panel and HR2 in the right panel. For source classification see Fig. 8. Sources with a 
statistical significance of the variability below 3 are marked in green. 



Figure 8 shows the variability factor plotted versus maxi- 
mum detected XID flux. Apart from XRBs, or XRBs in GlCs, 
or candidates of these source classes, which were selected based 
on their variability, there are a few SSS candidates showing pro- 
nounced temporal variability. The sources classified or identi- 
fied as AGN, background galaxies or galaxy clusters all show 
Fvar < 4. Most of the foreground stars show Fyar < 4. 

Out of the 1407 examined sources, we found 317 sources 
with a variability significance > 3.0, /. e. 182 more than reported 
in SPH2008. For bright sources it is much easier to detect vari- 
ability than for faint sources, because the diflference between 
the maximum observed flux and the detection limit is larger. 
Therefore the significance of the variability declines with de- 
creasing flux. This is illustrated by the distribution of the sources 
marked in green in Fig. 8. 

Table 9 lists all sources with a variability factor larger than 
five. There are 69 such sources (34 in addition to SPH2008). The 
sources are sorted in descending order with respect to their vari- 
ability factors. Table 9 gives the source number (Col. 1), maxima 
of flux variability (2) and maxima of the significance parameter 
(3). The next columns (4, 5) indicate the maximum observed flux 
and its error. Column 6 contains the class of the source. Sources 
with Fvar > 10 that were not already classified as SSS or fore- 
ground stars, were classified as XRB. 



Time variability can also be helpful to verify a SNR can- 
didate classification. If there is significant variability, the SNR 
classification must be rejected, and if an optical counterpart 
is detected, the source has to be re-classified as foreground 
star candidate. Column 7 contains references to the individual 
sources in the literature. In some cases the reference provides 
information on the temporal behaviour and a more precise clas- 
sification (see brackets). The numbers given in connection with 
Voss & Gilfanov (2007) and Williams et al. (2006b) are the vari- 
ability factors obtained in these papers from Chandra data. From 
the 69 sources of Table 9, ten show a flux variability larger than 
100. With a flux variability factor > 690 source N- 523 is the 
most variable source in our sample. Source N- 57 has the largest 
significance of variability, with a value of ^ 97. The variability 
significance is below 10 for just 33 sources, 15 of which show 
significance values below 5. Thirty-five of the variable sources 
are classified as XRBs or XRB candidates, and eight of them 
are located in globular clusters. Nine of the variable sources are 
SSS candidates, while six variable sources are classified as fore- 
ground stars and foreground star candidates. 

Table 10 lists all "bright" sources with a maximum flux 
larger than 8x10"^^ erg cm"^ s"^ and a flux variability smaller 
than five (the description of the columns is the same as in 
Table 9). All seven sources listed in Table 10 (three in addition 
to SPH2008) have a significance of variability > 10. Apart from 
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two sources, they are XRBs (three in globular clusters) or XRB 
candidates. The most luminous source in our sample is source 

966 with an absorbed 0.2-4.5 keV luminosity of ^ 3.3x10^^ 
erg s"^ at maximum. 

Figure 9 shows the relationship between the variability factor 
and the hardness ratios HRl and HR2, respectively. The hard- 
ness ratios are taken from Table 5. The HRl plot shows that the 
sample of highly variable sources includes SSS and XRB can- 
didates, which occupy two distinct regions in this plot (see also 
Haberl & Pietsch 1999, for the LMC). The SSSs marked by tri- 
angles, appear on the left hand side, while the XRBs or XRB 
candidates have much harder spectra, and appear on the right. It 
seems that foreground stars, SSSs and XRBs can be separated, 
on the HR2 diagram, although there is some overlap between 
foreground stars and XRBs. 

Individual sources are discussed in the Sects. 9 and 10. 



7. Cross-correlations with other M 31 X-ray 
catalogues 

Cross-correlations were determined by applying Eq. 4 to the 
sources of the XMM LP-total catalogue and to sources reported 
in earlier X-ray catalogues. The list of X-ray catalogues used is 
given in Table 4. 

7.1. Previous XMM-Newton catalogues 

Previous source lists based on archival XMM-Newton obser- 
vations were presented in Osborne et al. (2001), PFH2005, 
Orio (2006), SPH2008, and SBK2009. Of these four studies, 
PFH2005 covers the largest area of M 31. Table 11 lists all 
sources from previous XMM-Newton studies that are not de- 
tected in the present investigation. 

In the ten observations covering the major axis, and a field 
in the halo of M 31, PFH2005 detected 856 X-ray sources with 
a detection likelihood threshold of 7 (cf. Sect. 3.3). Of these 856 
sources, only 753 sources are also present in the XMM LP-total 
catalogue, /. e. 103 sources of PFH2005 were not detected. This 
can be due to: the search strategy; the parameter settings used in 
the eml detect run; the determination of the extent of a source 
for the XMM LP-total catalogue; the more severe screening for 
GTIs for the XMM LP-total catalogue, which led to shorter fi- 
nal exposure times; the use of the epreject task and last but 
not least due to the SAS versions and calibration files applied. 
The search strategy of PFH2005 was optimised to detect sources 
located close to each other in crowded fields. This point espe- 
cially explains the non-detection of the bright PFH2005 sources 
[PFH2005] 281, 312, 316, 327, 332, 384 (£ > 50) in the present 
study, as four of them ([PFH2005] 312, 316, 327, 332) are lo- 
cated in the innermost central region of M 31 where source 
detection is complicated by the bright diff'use X-ray emission, 
while [PFH2005] 281 and 384 He in the immediate vicinity 
of two bright sources ([PFH2005] 280 and 381 at distances of 
7. 7 ''and 5.5'', respectively). The changes in the SAS versions 
and the GTIs, in particular, aff'ect sources with small detection 
likelihoods (£ < 10). 

The improvements in the SAS detection tools and calibration 
files should reduce the number of spurious detections, which in- 
crease with decreasing detection likelihood. However, this does 
not necessarily imply that all undetected sources with £ < 10 of 
PFH2005 are spurious detections. The changes in the SAS ver- 
sions, calibration files and GTIs do not only aff'ect the source 
detection tasks, but also can cause changes in the background 



images. These changes may increase the assumed background 
value at the position of a source, which would result in a lower 
detection likelihood. Going from mlmin=7 to mlmin=6, but 
leaving everything else unchanged, we detected an additional 
nine sources of PFH2005. One of the previously undetected 
sources ([PFH2005] 75) was classified as <SSS>, but correlates 
with blocks of pixels with enhanced low energy signal in the 
PN off'set map and was corrected by epreject. Another source 
classified as <SSS> ([PFH2005] 799) is only detected in the 
MOS 1 camera, but not in MOS 2. From an examination by eye, 
it seems that source [PFH2005] 799 is the detection of some 
noisy pixels at the rim of the MOS 1 CCD 6 and not a real X-ray 
source. 

SPH2008 extended the source catalogue of PFH2005 by 
re-analysing the data of the central region of M 31 and 
also including data of monitoring observations of LMXB RX 
J0042.6-h4115. Of the 39 new sources presented in SPH2008, 
24 are also listed in the XMM LP-total catalogue, /. e. 15 
sources of SPH2008 were not detected. Diff'erences between 
the two studies include the detection likelihood thresholds 
used for eboxdetect (SPH2008: likemin=5) and emldetect 
(SPH2008: mlmin=6), the lower limit for the likelihood of ex- 
tention (SPH2008: dmlextmin = 4; XMM LP-total: 15), the 
screening for GTIs, the use of the epreject task and the SAS 
versions and calibration files used. Concerning the GTIs, im- 
ages, background images and exposure maps SPH2008 followed 
the same procedures as in PFH2005. The arguments given above 
are therefore also valid here. From the 14 undetected sources, 
three sources were detected in SPH2008 with mlmin < 7. One 
source ([SPH2008] 882) was added by hand to the final source 
list, as SPH2008 could not find any reason why emldetect did 
not automatically find it. The two extended sources ([SPH2008] 
863, 869) detected with extent likelihoods between 4.7 and 5.1 
in SPH2008, are neither detected as extended nor as pointlike 
sources in the present study, where the extent likelihood has to 
be larger than 15. 

SBK2009 re-analysed the XMM-Newton observations lo- 
cated along the major axis of M 31, ignoring all observations 
pointing to the centre of the galaxy. They used a detection likeli- 
hood threshold often. Of the 335 sources detected by SBK2009, 
304 sources are also contained in the XMM LP-total catalogue, 
/. ^.31 sources are not detected. Of the 304 re-detected sources, 
two ([SBK2009] 298, 233) are found with a detection likefi- 
hood below ten. Of the 31 undetected sources, 27 were also 
not detected in PFH2005. The remaining four sources correlate 
with PFH2005 sources, which were not detected in the present 
study. SBK2009 state that they find 34 sources not present in the 
source catalogue of PFH2005. A possible reason for this may 
be that SBK2009 used diff'erent energy bands for source detec- 
tion. They also had five bands, but they combined bands 2 and 3 
from PFH2005 into one band in the range 0.5-2 keV, and on the 
other hand they split band 5 of PFH2005 into two bands from 
4.5-7 keV and from 7-12keV, respectively. This might also ex- 
plain why most of the additional found sources were classified 
as <hard>. 

Orio (2006) addressed the population of SSSs and QSSs 
based on the same archival observations as PFH2005. Orio 
(2006) detected 15 SSSs, 18 QSSs and 10 SNRs of which 
one ([O2006] Table 4, Src.3) is also Hsted as an SSS 
([O2006] Table 2, Src. 13). Of these sources two SSSs, four 
QSSs and two SNRs (among them is the source [O2006] Table 4, 
Src. 3) are not contained in the XMM LP-total catalogue. These 
seven sources are also not present in the PFH2005 catalogue. 
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Table 11. Sources from previous XMM-Newton studies that are not listed in the XMM LP-total catalogue. 



PFH2005 856 sources 
103 not detected 
6 not detected, LH>100: 

10 not detected, 20<LH<50: 

14 not detected, 15<LH<20: 

21 not detected, 10<LH<15: 
52 not detected, LH<10: 



327 (<SNR>,LH=2140.0), 384 (XRB,667.0), 332 (<SNR>,654.0), 
316 (<SNR>,259.0), 312 (<SNR>,241.0), 281(<hard>,160.0) 
75 (<SSS>), 423 (<fg Star>), 120 (<hard>), 505 (<hard>), 

220 (<SNR>), 304 (<fg Star>), 819 (<hard>), 799 (<SSS>), 413 (<SNR>), 830 (<hard>) 
427(<hard>), 734 (<hard>), 424 (<hard>), 518 (<SSS>), 

232 (<hard>), 339 (<hard>), 446 (<SSS>), 219 (<fg Star>), 567 (<hard>), 256 (<fg Star>), 
356 (<hard>), 248 (<hard>), 160 (<hard>), 399 () 
375 (<hard>), 17 (<hard>), 195 (<hard>), 417 (<SNR>), 

783 (<hard>), 803 (<hard>), 829 (<hard>), 135 (<hard>), 151 (<hard>), 131 (<hard>), 
426 (<hard>), 593 (<fg Star>), 526 (<hard>), 250 (<hard>), 62 (<hard>), 67 (<hard>), 
188 (<hard>), 186 (<AGN>), 510 (<hard>), 529 (<hard>), 754 (<hard>) 
599 (<hard>), 439 (<hard>), 809 (<hard>), 14 (<SNR>), 743 (<hard>), 
433 (<hard>), 5 (), 210 (<hard>), 97 (<hard>), 708 (<hard>), 476 (), 534 (<hard>), 501 (), 
170 (<hard>), 146 (SNR), 769 (), 838 (<hard>), 571 (<hard>), 816 (<hard>, 554 (), 627 (<hard>), 
464 (<fg Star>), 811 (<hard>), 655 (<hard>), 184 (<hard>), 447 (<hard>), 380 (<hard>), 
566 (<hard>), 137 (<fg Star>), 63 (), 48 (), 152 (<fg Star>), 291 (<hard>), 559 (<hard>), 
102 (<hard>), 740 (<hard>), 540 (<fg Star>), 240 (<hard>), 485 (), 668 (<hard>), 44 (), 

560 (<hard>), 836 (<hard>), 436 (<hard>), 484 (<fg Star>), 216 (<hard>), 362 (<hard>), 527 (<>), 179 (<hard>), 
834 (<hard>), 86 (<hard>), 455 () 

SPH2008 39 sources 
15 not detected 

3 not detected, 50<LH< 100: 874 (<SNR>,LH=85.5), 895 (<hard>,75.9), 882 ( ,56.4) 

6 not detected, 10<LH<50: 869 (), 885 (<SNR>), 863 (<hard>), 875 (<SSS>), 893 (<hard>), 866 (<hard>) 

6 not detected, LH<10: 870 (<SNR>), 891 (<hard>), 889 (<hard>), 872 (<SNR>), 867 (<hard>), 862 (<SNR>) 



SBK2009 335 sources 
3 1 not detected 



4 (<hard>), 18 (<hard>), 29 (<hard>), 32 (<hard>), 34 (<hard>), 45 (<SSS>), 67 (<hard>), 
102 (<hard>), 106 (<hard>), 117 (<hard>), 149 (<hard>), 152 (<hard>), 179 (<hard>), 
183 (<hard>), 184 (<hard>), 188 (<hard>), 191 (<hard>), 192 (<AGN>), 202 (<hard>), 
204 (<fg star>), 217 (<hard>), 249 (<hard>), 250 (<hard>), 260 (<hard>), 274 (<hard>), 
279 (<hard>), 285 (<hard>), 295 (<hard>), 306 (<hard>), 325 (<hard>), 333 (<hard>) 



The nine bright variable sources from Osborne et al. (2001) 
were all detected. 

7.2. Chandra catalogues 

The Chandra catalogues used for cross-correlations were pre- 
sented in Sect. 1 (see also Table 4). 

Details of the comparison between the XMM LP-total cat- 
alogue and the different Chandra catalogues can be found in 
Table 12. Here, we only give a few general remarks. A non- 
negligible number of Chandra sources not reported in the XMM 
LP-total catalogue have already been classified as transient or 
variable sources. Thus, it is not surprising that these sources 
were not detected in the XMM-Newton observations (parts of: 
Voss & Gilfanov 2007; WilHams et al. 2006b, DKG2004). One 
Chandra source (nl-66) lies outside the field of M 31 covered 
by the XMM-Newton observations. For the innermost central re- 
gion of M 3 1 , the point spread function of XMM-Newton causes 
source confusion and therefore only Chandra observations are 
able to resolve the individual sources, especially if they are faint 
compared to the diff'use emission or nearby bright sources (Kong 
et al. 2002b; Kaaret 2002; WilHams et al. 2004; Di Stefano et al. 
2004; WilHams et al. 2006b; Voss & Gilfanov 2007). This ex- 
plains why a certain number of these sources are not detected in 
XMM-Newton observations. 

Of the 28 bright X-ray sources located in globular clusters 
(Di Stefano et al. 2002), two were not found in the XMM- 
Newton data (see Table 12). They are also not included in the 
source catalogue of PFH2005 and SPH2008. Hence, both ob- 
jects are good candidates for being transient or at least highly 
variable sources (cf. Sect. 10.4.2). Another study of the globu- 
lar cluster population of M 3 1 is presented by Trudolyubov & 
Priedhorsky (2004). Their work is based on XMM-Newton and 



Chandra data and contains 43 X-ray sources. Of these sources 
three were not found in the present study. One of them ([TP2004] 
1) is located well outside the field of M 31 covered by the Deep 
XMM-Newton Survey^^ The second source ([TP2004] 21) cor- 
relates with r3-7 1 , which is discussed above (see Di Stefano et al. 
(2002) in Table 12). The transient nature of the third source 
([TP2004] 35), and the fact that it was not observed in any 
XMM-Newton observation taken before 2004 was already re- 
ported by Trudolyubov & Priedhorsky (2004). The source was 
first detected with XMM-Newton in the observation from 31 
December 2006. 

7.3. ROSAT catalogues 

Of the 86 sources detected with ROSAT HRI in the cen- 
tral -34' of M 31 (PFJ93), aH but eight sources ([PFJ93] 
1,2,31,33,40,48,63,85) are detected in the XMM-Newton ob- 
servations. Six of these eight sources ([PFJ93] 1,2,31,33,63,85) 
have already been discussed in PFH2005 and classified as tran- 
sients. Sources [PFJ93] 40 and 48 correlate with [PFH2005] 312 
and 332, respectively, which are discussed in Sect. 7.1. In ad- 
dition to these eight sources, PFH2005 did not detect source 
[PFJ93] 51. This source was detected in the XMM-Newton ob- 
servations centred on RX J0042.6-I-41 15 and was thus classified 
as a recurrent transient (see SPH2008). 

In each of the two ROSAT PSPC surveys of M 31, 396 in- 
dividual X-ray sources were detected (SHP97 and SHL2001). 
From the SHP97 catalogue 130 sources were not detected. Of 
these sources 48 are located outside the FoV of our XMM- 
Newton M 31 survey. From the SHL2001 catalogue, 93 sources 



The source was observed with XMM-Newton on 11 January 2001. 
Obs. id.: 0065770101 
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Table 12. Sources detected in previous Chandra studies that are not present in the XMM LP-total catalogue. 
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Kong et al. (2002b) 204 sources 
58 not detected 

5 transient: r3-46,r3-43,r2-28,rl-23,rl-19 

20 variable: r3-53,r3-77,r3406,r3-76,r2-52,r2-31,r2-23,rl-31,r2-20,rl-24,rl-28,rl-27,rl-33,rl-21,rl-20,rl-7,r2-15,rl-17^ 
33 unclassified: r3402,r3-92,r3-51,r3-75,r3-91,r3-89,r3401,r3-88,r2-44,r2-55,r2-54,r3-32,r2-53,rl-30,r3-99,rl-22,rl-26,rl-18,^^ 
r2-41,r2-40,r3-71,r2-50,r2-49,r2-38,r3-97,r2-46,r3-12,r3-66,r3-104,r3-82,r3-5,r3-4 



Kaaret (2002) 142 sources 
26 not detected 

3 transient: J004217.0+41 1508,1004243.8+41 1604,1004245.9+41 1619 

7 variable: J004232.7+41 131 1,1004242.0+41 1532,1004243. 1+41 1640,J004244.3+41 1605,1004245.2+41 161 1,J004245.5+41 1608, 

J004248.6+411624 

16 unclassified: J004207.3+410443,J004229. 1+412857,1004239.5+41 1614,1004239.6+41 1700,1004242.5+41 1659,1004242.7+41 1503, 
J004243. 1+41 1604,1004244.2+41 1614,1004245.0+41 1523,1004246. 1+41 1543,1004247.4+41 1507,1004249. 1+41 1742, 
J004251.2+41 1639,1004252.3+41 1734,1004252.5+41 1328,1004318.5+410950 



Williams et al. (2004) 166 sources 
28 not detected 

12 transient: sl-79,sl-80,sl-82,r3-46,r2-28,rl-23,rl-19,r2-69,rl-28,rl-35,rl-34,nl-85 

7 variable: r2-31,rl-31,rl-24,rl-20,rl-7,rl-17,rl-16 

9 unclassified: sl-81,r2-68,sl-85,rl-30,rl-22,rl-26,rl-18,nl-77,nl-84 



Voss & Gilfanov (2007) 261 sources 
104 not detected 

11 transient: 6,12,29,32,41,51,59,84,118,130,146 

15 variable: 3,5,8,9,18,22,24,27,44,63,92,96,99,149,169 

78 unclassified: 4,19,21,25,26,30,37,39,40,42,48,49,53,56,57,58,60,62,65,70,73,75,76,77,80,82,84,86,87,89,91,94,97,98,104,109,114, 
115,117,119,122,124,129,133,138,141,143,144,145,150,152,158,162,164,167,171,173,182,183,188,189,191,193,194, 
197,202,205,206,210,213,217,219,220,225,256,257,263 



Williams et al. (2006b) 45 sources 
25 not detected 

25 transient: nl-26,nl-85,nl-86,nl-88,nl-89,rl-19,rl-23,rl-28,rl-34,rl-35,r2-28,r2-61,r2-62,r2-66,r2-69,r2-72,r3-43,r3-46,sl-18, 
sl-27,sl-69,sl-79,sl-80,sl-82,s2-62 



DKG2004 43 sources 
15 not detected 

9 transient: s2-62,sl-27,sl-69,sl-18,nl-26,r2-62,rl-35,r2-61,r2-66 
5 unclassified: s2-27,s2-10,nl-29,nl-46,r2-54 
1 not in FoV: nl-66 



Di Stefano et al. (2002) 28 sources 
2 not detcted 

2 unclassified: 17 (= r2-15), 28 (= r3-71) 



Notes: 

Variability information (transient, variable) is taken from the papers. "Unclassified" denotes sources which are not indicated as transient or variable sources in the papers. 



are not detected, 60 of which lie outside the XMM-Newton FoV. 
For information on individual sources see Table 13. 

Forty-four (out of 302) sources from SHP97 and 27 (out of 
293) sources from SHL2001, have 7?0^Ar detection likelihoods 
larger than 15, but are not listed in the XMM LP-total catalogue. 
These sources have to be regarded as transient or at least highly 
variable. 

7.4. Einstein catalogue 

The list of Einstein X-ray sources in the field of M 31 re- 
ported by TF91 contains 108 sources, with 81 sources taken 
from the Einstein HRI data with an assumed positional error of 
3'' (reported by Crampton et al. 1984), and 27 sources based 
on Einstein IPC data with a 45'' positional error. Applying 
the above mentioned correlation procedure to the Einstein HRI 
sources, 64 of these sources are also detected in this work and 
listed in the XMM LP-total catalogue, /. ^.17 sources are not 
detected ([TF91] 29, 31, 35, 39, 40, 43, 46, 50, 53, 54, 65, 66, 
72, 75, 78, 93, 96). For the Einstein IPC sources only the 1 cr 
positional error was used to search for counterparts among the 
XMM-Newton sources. Of the 27 Einstein IPC sources six re- 
main without a counterpart in our catalogue ([TF91] 15, 99, 
100, 106, 107, 108), of which [TF91] 15 and 108 are located 
outside the field of M 31 covered by the XMM LP-total cata- 



logue. Sources [TF91] 50 and 54 correlate with [PFH2005] 312 
and 316, respectively. Both sources were already discussed in 
Sect. 7.1. Apart from [TF91] 106, which is suggested as a pos- 
sible faint transient by SHL2001, the remaining 18 sources are 
also not detected by PFH2005. They classified those sources as 
transient. 

8. Cross-correlations with catalogues at other 
wavelengths 

The XMM LP-total catalogue was correlated with the catalogues 
and public data bases given in Sect. 3.8. Two sources (from the 
XMM LP-total and from the reference catalogues) were be con- 
sidered as correlating, if their positions matched within the un- 
certainty (see Eq. 4). 

However, the correlation of an X-ray source with a source 
from the reference catalogue does not necessarily imply that the 
two sources are counterparts. To confirm this, additional infor- 
mation is needed, like corresponding temporal variability of both 
sources or corresponding spectral properties. We should also 
take into account the possibility that the counterpart of the ex- 
amined X-ray source is not even listed in the reference catalogue 
used (due to faintness for example). 

The whole correlation process will get even more challeng- 
ing if an X-ray source correlates with more than one source from 



1 8 H. Stiele et al. : The deep XMM-Newton Survey of M 3 1 

Table 13. Sources from the ROSAT FSFC catalogues that are not present in the XMM LP-total catalogue. 



SHP97 396 sources 




130 not detected 




48 outside FoV: 


1,2,3,4,5,7,8,14,31,41,72,91,98,104,120,125,159,202,209,271,276,285,286,290,300,312,314, 




^9n lAl 2f^2 "Kf^l ^71 "KIA ^87 ^88 ^8Q ^QO ^09 ^Q^ ^Qzl ^Q'^ ^QA 


1 transient! 


69 


zi noi Qeiecieu, i^rKiz. 


1 Q 94 97 1^ Zl^ S9 '^Q ^8 71 1 1 ZlQ 1 ^1 lf\A 97^ 97'^ ^07 ^9Q ^^8 ^77 

iy,ZH,z / , J j,^o, jz, jy,Dj,oo, / 1, i j j, i^-y, idi,zdh,z i j,a i j, ju / , jzy, j ju, j jo, j / / 


10 noi QeiecieQ, izsi-/ri<iD. 


19 1 ZlQ 89 11^ 114 198 1 9^0 9^9 98^ ^^ZL ^Azl ^79 ^7A 

iz, 1 j,H-y,oz,yj, 1 1 J, 1 iH-, izo, iyo,zju,zoz,zo J, J jo^, J /z, j /o 


not QeiecieQ, i-/ri^ij. 


1 ^^^'T 14—9^ f\\ "Kl(^r\ 9^ zL^r'1 8 9^ ZlSf'SI 9^ ^nr'9n 1 ^ fsf^CXfs 9^ ^7('zLS^^ 9^ 78(^90 8n^'1 ^ 81 (lf\ f\\ 
lOi^i-zlrl— ZO.O ),jL\j\j.L ^,^Ji^lo.Z 1 .Z ^,DUi^ZU. 1 ),\)\i\j\).L ),V / [^-JJV.z. ), 1 oyZXj.D ^,oUi^lO. J ^,o H^ZO.O ), 




00(^jJ. / j,yji^J^O.Uj,lUZi^lO.^J,lZOi^Zl /.JJ,14H^o4J.JJ,l^Ji^40.yj,14Di^D/ J. / j,lDOi^l 1 l\y\J.\J), 




171 l^ 1 89(^4^4 Zl^ 1 8^^^Q 8^ 1 QO^^I 1 1 1 Q1 C^A 1 CtlC^A 9ni('1 HI ^^ 91 zLf'Zinn 9^ 91 '^^'9S1 

1 / H^Jh-. J ), loZi^M-J^.H- ), loDi^jy.o ), lyUi^l Lj.yj), ly l\p^■ .D ), lyA\p^■.J ^,Zl^i^^UU.Z ),AID\AD l.Kj), 




939n04 4') 24'Sr26 (W 260r'S4 ft\ 263n8 1 26'Sr24 6"! 268r'S4 270^40 4') 277n 5 6"! 309^81 K\ 






SHL2001 396 sources 




93 not detected 




60 outside FoV: 


1,2,3,4,5,6,7,8,9,10,11,12,14,15,16,21,22,32,39,58,67,69,75,77,81,83,85,90,93,125,141,146, 




160,164,192,202,243,260,282,296,298,302,325,326,328,355,371,372,378,379,383,388,389,390, 




391,392,393,394,395,396 


4 not detected, LH<12: 


62,96,238,269 


2 not detected, 12<LH<15: 


231,361 


27 not detected, LH>15: 


51(LH=28.4),104(901.2),121(94.1),126(46.2),143(34.7),168(131.9),171(43.0),173(317.8),190(215.8), 




207(98.0),208(298.8),226(73.1),230(75.6),232(1165.6),240(218.4),246(39.9),248(219.6),256(60.0), 




267(22.2),271(52.8), 322(2703. 3),324(147.7),344(40.7),356(15.3),365(19.0),380(17.4),384(15. 8) 



the reference catalogue. In this case we need a method to de- 
cide which of the correlating sources is the most likely to cor- 
respond to the X-ray source in question. Therefore, the method 
used should indicate how likely the correlation is with each one 
of the sources from the reference catalogue. Based on these like- 
lihoods one can define criteria to accept a source from the refer- 
ence catalogue as being the most likely source to correspond to 
the X-ray source. 

The simplest method uses the spatial distance between the 
X-ray source and the reference sources to derive the likelihoods. 
In other words, the source from the reference catalogue that is 
located closest to the X-ray source is regarded as the most likely 
source corresponding to the X-ray source. 

An improved method is a "likelihood ratio" technique, were 
an additional source property (e. g. an optical magnitude in deep 
field studies) is used to strengthen the correlation selection pro- 
cess. This technique was applied successfully to deep fields to 
find optical counterparts of X-ray sources (e. g. Brusa et al. 
2007). A drawback of this method is that one a priori has to 
know the expected probability distribution of the optical magni- 
tudes of the sources belonging to the studied object. In our case, 
this means that we have to know the distribution function for all 
optical sources of M 31 that can have X-ray counterparts, with- 
out including foreground and background sources. Apart from 
the fact that such distribution functions are unknown, an addi- 
tional challenge would be the time dependence of the magnitude 
of the optical sources (e. g. of novae) and of the connection be- 
tween optical and X-ray sources (e. g. optical novae and SSSs). 
Therefore it is not possible to apply this "likelihood ratio" tech- 
nique to the sources in the XMM LP-total survey. The whole 
correlation selection process becomes even more challenging if 
more than one reference catalogue is used. 

To be able to take all available information into account, we 
decided not to automate the selection process, but to select the 
class and most likely correlations for each source by hand (as it 
was done e. g. in PFH2005). Therefore the source classification, 
and thus the correlation selection process, is based on the cross 
correlations between the diff'erent reference catalogues, on the 
X-ray properties (hardness ratios, extent and time variability), 
and on the criteria given in Table 6. For reasons of completeness 
we give for each X-ray source the number of correlations found 



in the USNO-Bl, 2MASS and LGGS catalogues in Table 5. The 
caveat of this method is that it cannot quantify the probability of 
the individual correlations. 

9. Foreground stars and background objects 

9.1. Foreground stars 

X-ray emission has been detected from many late-type - spectral 
types F, G, K, and M - stars, as well as from hot OB stars (see 
review by Schmitt 2000). Hence, X-ray observations of nearby 
galaxies also reveal a significant fraction of Galactic stars. With 
typical absorption-corrected luminosities of Lo.2-iokeV <10^^ erg 
s"^ single stars in other galaxies are too faint to be detected 
with present instruments. However, concentrations of stars can 
be detected, but not resolved. 

Foreground stars (fg Stars) are a class of X-ray sources 
which are homogeneously distributed over the field of M 31 
(Fig. 10). The good positional accuracy of XMM-Newton and 
the available catalogues USNO-Bl, 2MASS and LGGS allow 
us to efficiently select this type of source. The selection crite- 
ria are given in Table 6. The optical follow-up observations of 
Hatzidimitriou et al. (2006) and Bonfini et al. (2009) have con- 
firmed the foreground star nature of bright foreground star candi- 
dates selected in PFH2005, based on the same selection criteria 
as used in this paper. Somewhat different criteria were applied 
for very red foreground stars, with an LGGS colour V - R > 1 
or USNO-Bl colour B2 - R2 > 1. These are classified as fore- 
ground star candidates, if /x//opt < -0.65 and /x//opt,R < -1.0. 
A misclassification of symbiotic systems in M 3 1 as foreground 
objects by this criterion can be excluded, as symbiotic systems 
typically have X-ray luminosities below 10^^ erg s"^ which is 
more than a factor 100 below the detection limit of our survey. 

If the foreground star candidate lies within the field cov- 
ered by the LGGS we checked its presence in the LGGS im- 
ages (as the LGGS catalogue itself does not list bright stars, 
because of saturation problems). Otherwise DSS2 images were 
used. Correlations with bright optical sources from the USNO- 
Bl catalogue, with an fx/ f opt in the range expected for fore- 
ground stars, that were not visible in the optical images were 
rejected as spurious. We found 223 foreground star candidates. 
Fourty sources were identified as foreground stars, either be- 
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Table 14. Infrared colours and spectral types of foreground stars that 
show flares. 







J mag 


H mag 


K mag 


SpT* 


err+ 




473 


12.984 


12.681 


12.558 


KO 


0.4 




714 


14.310 


13.618 


13.458 


MO 


0.2 




780 


14.251 


13.595 


13.351 


M3 


0.1 




1551 


12.666 


12.009 


11.806 


M2 


0.1 


* 

• 


1585 


13.488 


12.899 


12.650 


M2 


0.1 




1676 


10.460 


9.878 


9.798 


Kl 


0.2 




1742 


13.722 


13.138 


12.896 


Ml 


0.2 

















Notes: 

* : spectral type 

+ : error (in subtypes) 



- • 1 



Figure 10. The spatial distribution of foreground stars and candidates, 
classified in the XMM LP-total catalogue. The image shows the homo- 
geneous distribution of the sources over the covered field (marked with 
green dots). 



cause they are listed in the globular cluster catalogues as spec- 
troscopically confirmed foreground stars or because they have 
a spectral type assigned to them in the literature (Bonfini et al. 
2009; Hatzidimitriou et al. 2006, SIMBAD). 

Two of the foreground star candidates close to the centre of 
M 31 (N^ 826, 1 110) have no entry in the USNO-Bl and 
LGGS catalogues, and one has no entry in the USNO-Bl R2 
and B2 columns (N- 976). However, they are clearly visible on 
LGGS images, they are 2MASS sources and they fulfil the X- 
ray hardness ratio selection criteria. Therefore, we also classify 
them as foreground stars. 

The following 19 sources were selected as very red fore- 
ground star candidates: 54, 118, 384, 391, 393, 
585, 646, 651, 711, 1 038, 1 119, 1 330, 
1 396, 1 429, 1 506, 1 605, 1 695, 1 713 and 
1 747. A further 10 sources (N^ 210, 269, 278, 
310, 484, 714, 978, 1 591, 1 908 and 1 930) 
fulfil the hardness ratio criteria, but violate the /x//opt criteria 
and are therefore marked as "foreground star candidates" in the 
comment column of Table 5. 

Six sources (N^ 473, 780, 1 551, 1 585, 1 676, 
N- 1742), classified as foreground star candidates, have X- 
ray light curves that in a binning of 1 000 s showed flares (see 
Fig. 11). These observations strengthen the foreground star clas- 
sification. A seventh source (N- 714) is classified as a foreground 
star candidate, since its hardness ratios and its fx/ f opt ratio in the 
quiescent state fulfil the selection criteria of foreground star can- 
didates. In addition, the source shows a flare throughout obser- 
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Figure 12. Distribution of the source fluxes in the 0.2-4.5 keV (XID) 
band. The diagram shows a histogram of the number of foreground stars 
and candidates per flux bin, in logarithmic scales. 



vation ss3. Hence, the fx/ f opt ratio for this observation, in which 
the source is brightest, is too high to be consistent with the range 
of values expected for foreground stars. 

Table 14 gives the J, H and K magnitudes taken from the 
2MASS catalogue for each of the six flaring foreground stars. 
Using the standard calibration of spectral types for dwarf stars 
based on their near infrared colours (from the fourth edition of 
Allen's astrophysical quantities, edt. A.Cox, p. 151) we derived 
the spectral classification for the objects, using both H-K and 
J-K. The spectral types (and "error") we give in Table 14 are 
derived from averaging the two classes (derived from the two 
colours). The spectral types are entirely consistent with those 
expected for flare stars (usually K and M types). 

Figure 12 shows the XID flux distribution for foreground 
stars and foreground star candidates, which ranges from 
6.9x10"^^ erg cm"^ s"^ to 2.0x10"^^ erg cm"^ s"^ Most of the 
foreground stars and candidates (257 sources) have fluxes below 
5x10"^"^ erg cm"^ s"^ 
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Figure 11. X-ray light curves of foreground stars and candidates that, with a binning of 1000 s, show flares. 
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9.1 .1 . Comparing XMM-Newton, Chandra and ROSAT 
catalogues 

In the combined ROSAT PSPC survey (SHP97, SHL2001) 55 
sources were classified as foreground stars. Of these, 14 sources 
remain without counterparts in the present XMM-Newton sur- 
vey. Five of these 14 sources are located outside the field ob- 
served with XMM-Newton. Forty-one ROSAT foreground star 
candidates have counterparts in the XMM LP-total catalogue. 
Of these counterparts, 16 were classified as foreground star can- 
didates and four were identified as foreground stars (spectral 
type from Bonfini et al. 2009; Hatzidimitriou et al. 2006, or 
SIMBAD). In addition 12 sources were listed as <hard>, two 
as AGN candidates and one as a globular cluster candidate in 
the XMM LP-total catalogue. The counterparts of three ROSAT 
sources remain without classification in the XMM LP-total cat- 
alogue. 

Another three ROSAT sources have more than one counter- 
part in the XMM-Newton data. Source [SHP97] 109 correlates 
with sources N^ 597, N^ 604, N^ 606, and N^ 645. The former 
three are classified as <hard>, while source N- 645 is classi- 
fied as a foreground star candidate. However source N- 645 has 
the largest distance from the position of [SHP97] 109 compared 
to the other three XMM-Newton counterparts. Furthermore, this 
source had a flux below the ROSAT detection threshold (about 
a factor 2.6) in the XMM-Newton observations and is about a 
factor 3-34 fainter than the three other possible XMM-Newton 
counterparts. Thus it is very unlikely that [SHP97] 109 repre- 
sents the X-ray emission of a foreground star. 

Source [SHL2001] 156 has two XMM-Newton coun- 
terparts and is discussed in Sect. 10.1.3. The third source 
([SHL2001] 374) correlates with sources N^ 1 922 and N^ 1 924. 
The two XMM-Newton sources are classified as <hard> and 
as a foreground star candidate, respectively. In the source cat- 
alogue of SHL2001 source [SHP97] 369 is Hsted as the coun- 
terpart of [SHL2001] 374. The source in the first ROSAT survey 
has a smaller positional error and only correlates with source N- 
1 924. Although this seems to indicate that source N- 1 924 is 
the counterpart of [SHL2001] 374, we cannot exclude the pos- 
sibility that [SHL2001] 374 is a blend of both XMM-Newton 
sources, as these two sources have similar luminosities in the 
XMM-Newton observations. 

Kong et al. (2002b) classified four sources as foreground 
stars. For two sources (N^ 960=r2-42 and N^ 976=r3-33) the 
classification is confirmed by our study. The third source (N- 
1000=r2-19) remained without classification in the XMM LP- 
total catalogue, as it is too soft to be classified as <hard> and the 
optical counterpart found in the LGGS catalogue does not fulfil 
the /x//opt criteria. The fourth source (r2-46) was not detected in 
the XMM-Newton observations. 

The foreground star classification of three sources (si -74, sl- 
45, nl-82) in WilHams et al. (2004) is confirmed by the XMM 
LP-total study (N^ 289, N^ 603, N^ 1 449). For source N^ 289 
the spectral type FO was determined (Hatzidimitriou et al. 2006). 

The source list of DKG2004 contains six sources (s2-46, s2- 
29, s2-37, sl-45, sl-20, r3-122) that are classified as foreground 
stars. All six sources are confirmed as foreground star candidates 
by our XMM-Newton study (cf. Table 5). For source N- 696 
(=sl-20) Hatzidimitriou et al. (2006) obtained the spectral type 
GO. 

Of the four sources listed as foreground stars in Voss & 
Gilfanov (2007) only one source (N^ 936= [VG2007] 168) 
was confirmed as a foreground star, based on the entry in the 
RBCV3.5 and CaldweU et al. (2009). The second source (N^ 



1 118=[VG2007] 180) is Hsted in the RBC V3.5 and CaldweU 
et al. (2009) as a globular cluster. The third source (N- 
829=[VG2007] 181) does not have a counterpart in the USNO- 
Bl, 2MASS or LGGS catalogues, nor does it fulfil the hardness 
ratio criteria for foreground stars. Hence, the source is classified 
as <hard>. The fourth source ([VG2007] 81) is not spatially re- 
solved from its neighbouring source [VG2007] 79 in our XMM- 
Newton observations (source N- 1 078). Hence source N- 1 078 
is classified as <hard>. 

9.2. Galaxies, galaxy clusters and AGN 

The majority of background sources belong to the class of ac- 
tive galactic nuclei (AGN). This was shown by the recent deep- 
est available surveys of the X-ray background (Mushotzky et al. 
2000; Hasinger et al. 2001; Brandt & Hasinger 2005). The class 
of AGN is divided into many sub- sets. The common factor in all 
the sub- sets is that their emission emanates from a small, spa- 
tially unresolved galactic core. The small size of the emitting 
region is implied by the X-ray flux variability observed in many 
AGN, which is on time scales as short as several minutes (to 
years). The observed X-ray luminosities range from 10^^ to 10"^^ 
erg s"\ sometimes even exceeding 10^^ erg s"^ Although AGN 
show many difl'erent properties, like the amount of radio emis- 
sion or the emission line strengths and widths, they are believed 
to be only difl'erent facets of one underlying basic phenomenon 
(cf. Urry & Padovani 1995): the accretion of galactic matter onto 
a supermassive black hole (-10^-10^ M©) in the centre of the 
galaxy. 

It is difficult and, to some extent, arbitrary to distinguish be- 
tween active and normal galaxies, since most galaxies are be- 
lieved to host a black hole at the position of their kinetic centre 
(Bender et al. 2005). In normal galaxies the accretion rate to the 
central supermassive BH is so low, that only weak activity can be 
detected - if at all. The overall thermal emission of the nuclear 
region is due to bremsstrahlung from hot gas. The total X-ray 
luminosity of a normal galaxy can reach some 10^^ erg s"\ at 
most. It consists of diff'use emission and emission of unresolved 
individual sources. 

Galaxy clusters (GCls) are by far the largest and most mas- 
sive virialised objects in the Universe. Their masses lie in the 
range of 10^^- lO^^M© and they have sizes of a few megaparsecs 
(Mpc). A mass-to-Hght ratio of M/L^ 200 Mq/L© indicates that 
galaxy clusters are clearly dominated by their dark matter con- 
tent. Furthermore, galaxy clusters allow us to study the baryonic 
matter component, as they define the only large volumes in the 
Universe from which the majority of baryons emit detectable ra- 
diation. This baryonic gas, the hot intracluster medium (ICM), is 
extremely thin, with electron densities of ^ 10^-10^ m"^, and 
fills the entire cluster volume. Owing to the plasma temperatures 
oik^T ^ 2-10 keV, the thermal ICM emission gives rise to X-ray 
luminosities of Lx - 10^^- 3x10^^ erg s"^ Therefore galaxy clus- 
ters are the most X-ray luminous objects in the Universe next to 
AGN. 

We identified four sources as background galaxies and 1 1 
as AGN, and classified 19 galaxy and 49 AGN candidates. 
The classification is based on SIMBAD and NED correlations 
and correlations with sources listed as background objects in 
the globular cluster catalogues (RBC V3.5 and Caldwell et al. 
2009). Sources are classified as AGN candidates, if they have 
a radio counterpart (NVSS; Braun 1990; Gelfand et al. 2004) 
with the additional condition of being neither a SNR nor a 
SNR candidate from X-ray hardness ratios, as well as not being 
listed as a "normal" background galaxy in Gelfand et al. (2004). 
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Figure 13. The spatial distribution of background sources and candi- 
dates, classified in the XMM LP-total catalogue. AGN are marked with 
blue dots, "normal" galaxies with red dots and galaxy clusters with 
green dots. 

Most AGN will be classified as <hard> ((HR2-EHR2)> -0.2, 
see Table 6) because of their intrinsic power law component. 
Additional absorption in the line of sight by the interstellar 
medium of M 31 will lead to an even higher HR2. Only the few 
AGN with a dominant component in the measured flux below 
1 keV may lead to a classification <SNR> or <fg Star> in our 
adapted scheme. 

One (N- 995) of the four identified galaxies is M 32. An 
overview of previous X-ray observations of this galaxy is given 
in PFH2005. They also discuss the fact that Chandra resolved 
the X-ray emission of M 32 into several distinct point sources 
(maximum separation of the three central Chandra sources 8V3). 
Although M 32 is located closer to the centre of the FoV in the 
observations of field SSI, than it was in the si observation used 
in PFH2005, XMM-Newton still detects only one source. The 
remaining three sources (N- 88, N- 403, N- 718) are identified 
as galaxies, because they are listed as background galaxies in 
both the RBC V3.5 and Caldwefl et al. (2009). For source N^ 
403 (B 007) NED gives a redshift of 0.139692 ± 0.000230 (Kim 
et al. 2007). 

Eleven X-ray sources are identified as AGN. The first one 
(N- 363) correlates with a BL Lac object located behind M 31 
(NED, see also PFH2005). The second source (N^ 745) cor- 
relates with a Seyfert 1 galaxy (5C 3.100), which has a red- 
shift of ^ 0.07 (SIMBAD). The third source (N^ 1 559) cor- 
relates with a quasar (Sharov 21) that showed a single strong 
optical flare, during which its UV flux has increased by a fac- 
tor of ~20 (Meusinger et al. 2010). The remaining sources were 



Table 15. Spectral fit parameters for extended sources 



Src ID 


TVh/IO^i cm-2 




redshift 


xVdof 


141 
252 
304 
1543 


1 1Q+1.63 
^•^^-0.88 

2.68!?;g 

9 74+6-91 


9 17+2.30 
'-0.68 
1 QC+0.64 
^•^^-0.29 

2.08:2-3J 


24+^-24 
22+^-^^ 

^•^^-0.07 
-0.05 
^•'^^-0.26 


78.5/53 
56.4/151 
50.9/57 
32.9/34 



spectroscopically confirmed (from our optical follow-up obser- 
vations) to be AGN (D. Hatzidimitriou, private communication; 
and Hatzidimitriou et al. (2010) in prep.). 

In Sect. 5.4 the 12 extended sources in the XMM LP-total 
catalogue were presented. Kotov et al. (2006) showed that the 
brightest of these sources (N- 1795) is a galaxy cluster located 
at a redshift of z = 0.29. For the remaining 11 sources. X-ray 
spectra were created and fitted with the MEKAL model in XSPEC. 
Unfortunately, for most of the examined sources the spectral pa- 
rameters (foreground absorption, temperature and redshift) are 
not very well constrained. Nevertheless four sources (N- 141, 
N^ 252, N^ 304, N^ 1543) with temperatures in the range of 
~ 1-2 keV and proposed redshifts between 0.1-0.6 were found 
(Table 15). Inspection of optical images (DSS2 images and if 
available LGGS images) revealed an agglomeration of optical 
sources at the positions of these four extended X-ray sources. 
Thus they are classified as galaxy cluster candidates. 

Although, B242 (the optical counterpart of source N- 304) 
is listed as a globular cluster candidate in the RBC3.5 catalogue, 
Caldwell et al. (2009) classified this source as a background ob- 
ject. Our findings from the X-rays favour the background ob- 
ject classification. Hence a globular cluster classification for this 
source seems to be excluded. 

Source N- 1912 was already classified as a galaxy clus- 
ter candidate in PFH2005. The spectrum confirms this classi- 
fication. The best fit parameters are A^h= 1.29^^4^x10^^ cm"^, 
T = 2.8^^;^ keV and redshift of 0.06!^-^3^ 

A plot of the spatial distribution of the classified / identified 
background sources is given in Fig. 13, which shows that these 
sources are rather homogeneously distributed over the observed 
field. However, in the fields located along the major axis of 
M 3 1 we mainly see AGN, which are bright enough to be visible 
through M 3 1 , while most of the galaxies and galaxy clusters are 
detected in the outer fields. 

9.2.1. Comparing XMM-Newton, Chandra and ROSAT 
catalogues 

Of the ten ROSAT FSFC survey sources classified as background 
galaxies one is located outside the field of the Deep XMM- 
Newton Survey. The remaining objects are confirmed to be back- 
ground sources and are classified or identified as galaxies or 
AGN. The only case which is worth discussing in more detail 
is the source pair [SHP97] 246 and [SHL2001] 252. From the 
XMM-Newton observations it is evident that this source pair 
is not one source, as indicated in the combined ROSAT PSPC 
source catalogue (SHL2001), but consists of three individual 
sources (N^ 1269, N^ 1279 and N^ 1280). [SHL2001] 252 
correlates spatially with all three XMM-Newton sources, while 
[SHP97] 246 correlates only with source N^ 1 269, which is 
identified as a foreground star of type K2 (SIMBAD). The two 
other XMM-Newton counterparts of [SHL2001] 252 are classi- 
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fied as a galaxy candidate and an AGN candidate, respectively. 
In summary, [SHL2001] 252 is most likely a blend of both back- 
ground sources and maybe even a blend of all three XMM- 
Newton sources, while [SHP97] 246 seems to be the X-ray coun- 
terpart of the foreground star mentioned above. 

Kong et al. (2002b) classified source r3-83 (N^ 1 132) as an 
extragalactic object, as it is listed in SIMBAD and NED as an 
emission line object. Following PFH2005, we classified source 
N^ 1 132 as <hard>. The BL Lac object (N^ 363) was also de- 
tected in Chandra observations (Williams et al. 2004). 



10. M 31 sources 

10.1. Supersoft sources 

Supersoft source (SSS) classification is assigned to sources 
showing extremely soft spectra with equivalent blackbody tem- 
peratures of ~15-80eV. The associated bolometric luminosities 
are in the range of 10^^-10^^ erg s"^ (Kahabka & van den Heuvel 
1997). 

Because of the phenomenological definition, this class is 
likely to include objects of several types. The favoured model 
for these sources is that they are close binary systems with a 
white dwarf (WD) primary, burning hydrogen on the surface (cf. 
Kahabka & van den Heuvel 1997). Close binary SSSs include 
post-outburst, recurrent, and classical novae, the hottest symbi- 
otic stars, and other LMXBs containing a WD (cataclysmic vari- 
ables, CVs). Symbiotic systems, which contain a WD in a wide 
binary system, may also be observed as SSSs (Kahabka & van 
den Heuvel 1997). Because matter that is burned can be retained 
by the WD, some SSS binaries may be progenitors of type-la 
supernovae (cf. van den Heuvel et al. 1992). 

The XMM LP-total catalogue contains 30 SSS candidates 
that were selected on the basis of their hardness ratios (see Fig. 6 
and Table 6). 

10.1.1. Spatial and flux distribution 

Figure 14 shows the spatial distribution of the SSSs. Clearly 
visible is a concentration of sources in the central field. There 
are two explanations for that central enhancement. The first is 
that the central region was observed more often than the remain- 
ing fields and therefore there is a higher chance of catching a 
transient SSS in outburst. The second reason is that the major 
class of SSSs in the centre of M 31 are optical novae (PFF2005, 
PHS2007). Optical novae are part of the old stellar population 
which is much denser in the centre of M 3 1 . 

Figure 15 gives the distribution of 0.2-l.OkeV source fluxes 
for all SSSs (black) and for those correlating with optical no- 
vae (blue). The unabsorbed fluxes were determined assuming 
a 50 eV blackbody model (PFF2005). The two brightest SSSs 
(Fx >10"^^ erg cm"^ s"^) consist of a persistent source with 
217 s pulsations (N^ 1061; Trudolyubov & Priedhorsky 2008) 
and the nova M31N 2001-1 la (N^ 1 416; Smirnova et al. 2006). 
A large fraction of SSSs are rather faint, with fluxes below 
5x10"^"^ ergcm"^ s"^. Four sources have absorption-corrected 
luminosities below 10^^ erg s"^ (0.2-l.OkeV), which was indi- 
cated as the limiting luminosity for SSSs. That does not neces- 
sarily imply that these sources are not SSSs, since it is possible 
that the blackbody fit chosen does not represent well the proper- 
ties of these sources. A higher absorption or a lower temperature 
would lead to increased unabsorbed luminosities. We also have 
to take into account that we might have observed the source dur- 





ISarcmin 



Figure 14. The spatial distribution of SSSs classified in the XMM LP- 
total catalogue. The positions of the SSSs are marked with red and green 
dots. Sources that correlate with optical novae are given in green. An 
enhancement of sources in the central field is clearly visible. 

ing a phase of rising or decaying luminosity, /. e. not at maximum 
luminosity. 

10.1 .2. Correlations with optical novae 

By cross-correlating with the nova catalogue indicated in 
Sect. 3.8, 14 of the 30 SSSs can be classified as X-ray coun- 
terparts of optical novae. Of these 14 novae, eight (N- 748, N- 
993, 1 006, 1 046, 1 051, 1 076, 1 100, and 
1 236) are already discussed in PFF2005 and PHS2007. Nova 
M31N 2001-1 la was first detected as a supersoft X-ray source. 
Motivated by that SSS detection, Smirnova et al. (2006) found an 
optical nova at the position of the SSS in archival optical plates 
which had been overlooked in previous nova searches. Nova 
M31N 2007-06b has been discussed in Henze et al. (2009a). The 
remaining four novae are discussed individually in more detail 
below. 

As was shown in the XMM-Ncwton/ Chandra M 31 nova 
monitoring project^^, it is absolutely necessary to have a homo- 
geneous and dense sample of deep optical and X-ray observa- 
tions in order to study optical novae and their connections to su- 
persoft X-ray sources. In the optical, the outer regions of M 31 
are regularly observed down to a limiting magnitude of ~ 17 mag 
(Texas Supernova Search (TSS); Quimby 2006), while in X-rays 
only "snapshots" are available. Hence, the correlations of optical 

http : //www . mpe . mpg . de/~ m3 lnovae/opt/m3 1/M3 l_table . 
html 
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Figure 15. Distribution of the source fluxes in the 0.2-1.0 keV band. 
The diagram shows the number of SSSs per flux bin plotted versus the 
flux in logarithmic scale. The blue histogram gives the distribution of 
SSSs correlating with optical novae. 



80 



> 



20 







M31N 1997-lOc..... - 






^^^^ .a,.\0-''-'" 


y 




....\-\tf •••• ^ 






..... • ■■ 

yk - 

3-^0 

MO*-- — -^^"^"-"^^-^-^ 


■::----^^^y 

^ ■...wo*"a-' 

O 1 

^ 1 


1 





0.8 



2.4 



3.2 



Figure 16. Column density-temperature confidence contours inferred 
from the fit to the XMM-Newton EPIC PN spectrum of M31N1997- 
10c. The formal best fit parameters are indicated by the star. Also drawn 
are lines of constant bolometric luminosity (in erg s~^). The vertical 
dashed line indicates the Galactic foreground absorption in the direc- 
tion ofMSl. 



Table 16. 3cr upper limits for the absorption-corrected luminosities for 
Nova M31N 1997-lOc 



novae with detected SSSs have to be regarded as lucky coinci- 
dences. That also means that the identified nova counterparts are 
detected at a random stage of their SSS evolution which does 
not allow us to constrain the exact start or end point of the SSS 
phase, nor the maximum luminosity of the SSS. We also can- 
not exclude the possibility that some of the SSSs observed in the 
outer parts of M 3 1 correspond to the supersoft phase of opti- 
cal novae for which the optical outburst was missed. In the outer 
regions of M 31, the samples of optical novae and X-ray SSSs 
are certainly incomplete, due to the rather high luminosity limit 
in the optical monitoring, and the lack of complete monitoring 
in X-rays, respectively. So one should be cautious in deriving 
properties of the disc nova population of M 3 1 from the avail- 
able data. 

Nova M31 N 1 997-1 Oc was detected on 2 October 1997 at a B- 
band magnitude of 16.6 (ShA 58; Sharov & Alksnis 1998). An 
upper limit of 19 mag on 29 September 1997 was reported by 
the same authors. They classified this source as a very fast nova. 
In the XMM-Newton observation cl (25 June 2000), an SSS (N^ 
871), located within ~1 V9 of the optical nova, was detected. The 
source was fitted with an absorbed blackbody model. The formal 
best fit parameters of the XMM-Newton EPIC PN spectrum are: 
absorption A^h ~ 3.45x10^^ cm'^ and feT ^ 41 eV. The unab- 
sorbed luminosity in the 0.2-1 keV band is ^ 5.9x10^^ erg s"^. 
Confidence contours for absorption column density and black- 
body temperature are shown in Fig. 16. In the subsequent XMM- 
Newton observation of that region taken about half a year later 
(c2; 27 December 2000) the source is not detected. Although the 
source position is covered in observations c3 (29 June 2001), 
c4 (6/7 January 2002) and b (16-19 July 2004) the source was 
not re-detected. Using the count rates derived for the variabil- 
ity study (see Sect. 6) and assuming the same spectrum for the 
source as in observation cl, upper limits of the source luminosity 
can be derived, which are given in Table 16. 



observation Lx/10^^ erg s~^ (0.2-1 .0 keV) 

c2 10.8+ 

c3 1.9 

c4 1.0 



Notes: 

+ : The count rate detected in observation c2 gives a luminosity of 
2.4+2.8X 10^^ erg s~\ which results in the upper limit given in the Table. 
The fact that this upper limit is higher than the luminosity detected in 
observation cl is, at least in part, attributed to the very short eff'ective 
observing time of less than 6 000 s. 



Nova M31 N 2005-01 b was discovered on 19 January 2005 at a 
white light magnitude of 16.3 by R. Quimby.^^ An SSS (N^ 764) 
that correlates with the optical nova (distance: 4V3; 3(T error: 
5V5) was found in observation ss2 taken on 8 July 2006, which 
is 535 days after the discovery of the optical nova. Due to the 
severe background screening applied to observation ss2, there is 
not enough statistics to obtain a spectrum of the X-ray source. 
To get an estimate of the spectral properties of that source we 
created a spectrum in the 0.2-0.8 keV range of the unscreened 
data. Although the spectrum was background corrected, we can- 
not totally exclude a contribution from background flares. The 
spectrum is best fitted by an absorbed blackbody model with 
an absorption of A^h ~ 1.03x10^^ cm"^ and a blackbody tem- 
perature of ~ 45 eV. The unabsorbed 0.2-1 keV luminosity 
is Lx -1.0x10^^ erg s"^ In another XMM-Newton observation 
taken 1 073 days after the optical outburst (ss21; 28 December 
2007) the X-ray source is no longer visible. The 3(T upper limit 
of the unabsorbed source luminosity is ~ 3.3x10^^ erg cm"^ s"^ 
in the 0.2-4.5 keV band, assuming the spectral model used for 
source detection. 
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Figure 17. XMM-Newton EPIC spectrum of nova M31N 2005-Olc. 
The absorbed black body fit to the data is shown in the upper panel. 
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Figure 18. Column density (A^h) - temperature (^b T) confidence con- 
tours inferred from the blackbody fit to the XMM-Newton EPIC spec- 
trum of M31N 2005-Olc (see Fig. 17). The formal best fit parameters 
are indicated by the star. Also drawn are lines of constant bolometric lu- 
minosity and the vertical dashed line indicates the Galactic foreground 
absorption (see Fig 16). 



Nova M31 N 2005-01 c was discovered on 29 January 2005 at 
a white light magnitude of 16.1 by R. Quimby.^^ In the XMM- 
Newton observation from 02 January 2007 (ns2, 703 days after 
optical outburst) an SSS was detected (N- 1 675) at a position 
consistent with that of the optical nova (distance: 0V9). The X- 
ray spectrum (Fig. 17) can be well fitted by an absorbed black- 
body model with the following best fit parameters: absorption 
A^H = 1.58+^-^^x1021 cm-2 andy^B?^ = 40±6eV. The unabsorbed 
0.2-1 keV luminosity is Lx ~ 1.2x10^^ erg s"^ Confidence con- 
tours for absorption column density and blackbody temperature 
are shown in Fig. 18. 

Nova M31 N 2005-09b was discovered in optical images taken 
on 01 and 02 September 2005 at white light magnitudes of ~ 1 8.0 
and -16.5 respectively. From 31 August 2005, an upper limit of 
-18.7 mag was reported (Quimby et al. 2005). The nova was 
spectroscopically confirmed (Pietsch et al. 2006) and classified 



as a possible Fell or hybrid nova^^. An X-ray counterpart (N- 
92) was detected in the XMM-Newton observation s3 (299 days 
after the optical outburst). Its position is consistent with that of 
the optical nova (distance: 0V57). As observation s3 was heav- 
ily aflfected by background flares, we only could estimate the 
spectral parameters from the unscreened data (see also para- 
graph about Nova M3 IN 2005-Olb). A blackbody fit of the 0.2- 
0.8 keV gives A^h - 2.7x10^^ cm'^ kr ^ 35 eV, and an unab- 
sorbed 0.2-1 keV luminosity of Lx -5.4x10^^ erg s"^. The X- 
ray source was no longer visible in observation s3 1 , which was 
taken 391 days after observation s3. 

10.1.3. Comparing XMM-Newton, Chandra and ROSAT 
catalogues 

The results and a detailed discussion of a study of the long- 
term variability of the SSS population of M 31 are presented 
in Stiele et al. (2010). In summary our comparative study of SSS 
candidates in M 31 detected with ROSAT, Chandra and XMM- 
Newton demonstrated that strict selection criteria have to be ap- 
plied to securely select SSSs. It also underlined the high vari- 
ability of the sources in this class and the connection between 
SSSs and optical novae. 



10.2. Supernova remnants 

After an supernova explosion the interaction between the ejected 
material and the ISM forms a supernova remnant (SNR). The 
SNR X-ray luminosities typically vary between 10^^ and 10^^ 
erg s-i (0.2-10 keV). 

SNRs can be divided into two categories, (i) sources where 
the thermal components dominate the X-ray spectrum below 
2 keV, and (ii) the so-called "plerions" or Crab-like SNRs with 
power law spectra. The former are located in areas of the X-ray 
colour/colour diagrams that overlap only with foreground star 
locii. If we assume that we have identified all foreground star 
candidates from the optical correlation and inspection of the op- 
tical images, the remaining sources can be classified as SNR can- 
didates using the criteria given in Table 6. Similar criteria were 
used to select supernova remnant candidates in XMM-Newton 
observations of M 33 (Pietsch et al. 2004; Misanovic et al. 2006). 
Ghavamian et al. (2005) and Long et al. (2010) confirmed the 
supernova remnant nature of many of these candidates based on 
optical and radio follow-up observations. They also used a hard- 
ness ratio criterion to select supernova remnant candidates from 
Chandra data. 

An X-ray source is classified as a SNR candidate if it either 
fulfils the hardness ratio criterion given in Table 6 (these are 25 
such sources), or if it correlates with a known optical or radio 
SNR candidate (six sources). The sources assigned the classifi- 
cation of a SNR candidate based on the latter criterion alone, are 
marked in the comment column of Table 5 with the flag 'only 
correlation'. As these six SNR candidates would be classified as 
<hard> on the basis of their hardness ratios, they are good can- 
didates for being "plerions". SNRs are taken as identified when 
they coincide with SNR candidates from the optical or radio and 
fulfil the hardness ratio criterion. For a discussion of detection 
of SNRs in diflTerent wavelength bands see Long et al. (2010). 
All together, we identified 25 SNRs and 31 SNR candidates in 
the XMM LP-total catalogue. 

This number is in the range expected from an extrapola- 
tion of the X-ray detected SNRs in the Milky Way as shown 
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Figure 19. Distribution of SNR fluxes in the 0.2-4.5 keV (XID) band. 
The diagrams show the number of identified and classified SNRs at each 
flux bin, plotted versus the flux. 



below. Assuming that our own Galaxy contains about 1 440 X- 
ray sources which are brighter than ~lxlO^^ erg s"^ (Revnivtsev 
et al. 2006), and that it contains -110 SNRs detected in X- 
rays (Green 2009), we would expect to detect ~50 SNRs in the 
XMM LP-total catalogue (0.4 x (1 897sources - 263fgStars)). 
This number is in good agreement with the number of identified 
and classified SNRs. 

The XID fluxes for SNRs range between 5.9x10-^^ 
ergcm"^ s"^ for source N- 1234 and 1.5x10"^^ erg cm"^ s"^ 
for source N- 419. These fluxes correspond to luminosities of 
4.3x10^^ ergs"^ to 1.1x10^^ erg s"^ A diagram of the flux 
distribution of the detected SNRs and candidates is shown in 
Fig. 19. 

Among the 25 identified SNRs, there are 20 SNRs from the 
PFH2005 catalogue. Source [PFH2005] 146, which correlates 
with the radio source [B90] 11 and the SNR candidate BA146, 
was not found in the present study. Source [SPH2008] 858, 
which coincides with a source reported as a ring-like extended 
object in Chandra observations that was also detected in the 
optical and radio wavelength regimes and identified as a SNR 
(Kong et al. 2003), was re-detected (N^ 1 050). Of the 31 SNR 
candidates ten have been reported by PFH2005. In the follow- 
ing, we first discuss in more detail the remaining four identified 
SNRs, that appear in the new catalogue but were not included in 
PFH2005: 

XMMM31 J003923.5+404419 (N^ 182) was classified as a 
SNR candidate from its [SlI]:Ha ratio. It appears as an 'irreg- 
ular ring with southerly projection' (Dodorico et al. 1980, and 
Fig. 20) and correlates with a radio source (Pooley 1969). X-ray 
radiation of that source was first detected in the present study. 




Figure 20. Ua, R, SII and OlII images, taken from the LOG 
Survey. Over-plotted is a circle at the position of source 
XMMM31 J003923.5+404419 with a radius of 5^5 (3o- posi- 
tional error of the X-ray source). The ring-like SNR is clearly visible in 
the Ra and S II bands. 




Figure 21. Ha, R, SII and OIII images, taken from the LGG 
Survey. Over-plotted is a circle at the position of source 
XMMM31 J004413.5+411954 with a radius of 376 (3cr positional 
error of the X-ray source). The SNR 'appears as a bright knot'. 



Braun & Walterbos (1993). The source has counterparts in the 
radio (Braun 1990) and X-ray (SHP97) range. It was reported as 
a SNR by SHP97. 

XMMM31 J00451 0.5+41 3251 and 

XMMM31 J00451 2.3+420029 (N^ 1587 and N^ 1593, 
respectively) are new X-ray detections and correlate with the 
radio sources: #354 and #365 in the list of Braun (1990). Source 
N- 1 587 also correlates with source 37W209 from the catalogue 
of Walterbos et al. (1985). No optical counterparts were reported 
in the literature. 

In the following, we discuss two SNR candidates in more 
detail: 



XMMM31 J00441 3.5+41 1954 (N^ 1410) was classified as a 
SNR candidate from its [SlI]:Ha ratio (Braun & Walterbos 
1993; Magnier et al. 1995). From Fig. 21 we can see that the 
source 'appears as a bright knot\ as was already reported by 



XMMM31 J004434.8+412512 (N^ 1481) lies in the periph- 
ery of a super-shell with [SII]:Hq^ >0.5 (Braun & Walterbos 
1993, src 490). Located next to this source is a SNR candi- 
date reported in Magnier et al. (1995, src 3-086), which has a 



H. Stiele et al: The deep XMM-Newton Survey of M 31 



27 



• . 284 

^ # ^ alpha^ 


• 

1481 f 

- • o- ■ 

• 

R band ^ ^ 


Sll band 

• 


*OIII band 
• 

- A ' 1481 t 4 



Figure 22. Ha, R, SII and OIII images, taken from the LGG 
Survey. Over-plotted is a blue circle at the position of source 
XMMM31 J004434.8+412512 with a radius of 5'! 9 (3(T positional er- 
ror of the X-ray source). Source [SPH97] 284 is indicated by a black 
circle with a radius of 21'' (3cr positional error), source 3-086 by the 
magenta circle with a radius of W\ the position of the radio counter- 
part is marked by the yellow circle. 



radio counterpart from the NVSS catalogue. N- 1 481 also cor- 
relates with ROSAT source [SPH97] 284, which was identified 
as a SNR in SPH97 due to its spatial correlation with source 3- 
086. Figure 22 shows the XMM-Newton error circle over-plotted 
on LGGS images. From the XMM-Newton source position it 
looks more likely that the X-rays are emitted from the HII re- 
gion rather than from the SNR candidate visible in the optical 
and radio wavelengths. Nevertheless the XMM-Newton source 
detected is point-like and its hardness ratios lie in the range ex- 
pected for SNRs. If the X-ray emission originated from the H ii- 
region, it should have been detected as spatially extended emis- 
sion. Thus, N- 1 481 is classified as SNR candidate. A puzzling 
fact, however, is the pronounced variability between ROSAT and 
XMM-Newton observations of Fyar = 9.82 with a significance 
of 5'var ~ 4 (see Table 17), which is not consistent with the long 
term behaviour of SNRs. There is still the possibility that the de- 
tected X-ray emission does not belong to either the H ii-region 
or a SNR at all. 



XMMM31 J004239.8+40431 8 (N^ 969) was already observed 
with ROSAT (SHP97, SHL2001) and Chandra (WilHams et al. 
2004, si -84). No optical counterpart is visible on the LGGS im- 
ages. The X-ray spectrum, which is shown in Fig. 23, is well 
fitted by an absorbed non-equilibrium ionisation model with the 
following best fit values: an absorption of A^h = 1-76^q^qX10^^ 
cm"^, a temperature of RbT = 219^^9eV, and an ionisation 
timescale of r = 1.75+^-^^ x lOScm'^ The unabsorbed 0.2- 
5keV luminosity is Lx ~ 6.5x10^^ erg s"^. The soft spectrum 
with the temperature of -200 eV is in good agreement with spec- 
tra of old SNRs e. g. in the SMC (Filipovic et al. 2008). Although 
the unabsorbed luminosity is rather high for an old SNR, it is still 
in the range found for other SNRs (cf. Kong et al. 2002a; Gaetz 
et al. 2007). Hence, XMMM31 J004239.9+404318 is classified 
as a SNR candidate. 



XMMM31 J004239.8+404318 XMM-Newton, EPIC 
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Figure 23. 0.2-3.0 keV EPIC spectrum of source 969. The best fit 
absorbed non-equilibrium ionisation model is indicated by the solid 
lines. 



10.2.1 . Comparing SNRs and candidates in XMM-Newton, 
Chandra and ROSAT catalogues 

The second ROSAT PSPC catalogue (SHL2001) contains 16 
sources classified as SNRs. The counterparts of 12 of these 
sources are also classified as SNRs or SNR candidates in the 
XMM LP-total catalogue. 

Table 17 fists the XMM-Newton, ROSAT, and Chandra 
fluxes of all SNRs and SNR candidates from the XMM LP-total 
catalogue that have counterparts classified as SNRs in ROSAT or 
Chandra source lists. In addition, the maximum flux variability 
and the maximum significance of the variability (following the 
variability calculation of Sect. 3.6) are given. Three SNRs that 
have ROSAT counterparts show variability changing in flux by 
more than a factor of five. The most variable source (N- 1 066) 
is discussed below, the second source was discussed in Sect. 10.2 
(XMMM31 J004434.8-h412512, N^ 1 481), and the third source 
(N- 1 234) is embedded in the diff'use emission of the central 
area of M 3 1 . In this environment the larger PSF of ROSAT re- 
sults in an overestimate of the source flux, since the contribution 
of the difl'use emission could not be totally separated from the 
emission of the point source. 

The remaining four ROSAT sources classified as SNRs and 
their XMM-Newton counterparts are discussed in the following 
paragraph. 

SHP97 report that [SHP97] 203 and [SHP97] 211 
(=[SHL2001] 206) correlate with the same SNR 
([DDB80] 1.13), have the same spectral properties and 
have luminosities within the range of SNRs. A correlation with 
the ROSAT HRI catalogue (PFJ93) reveals that the true X-ray 
counterpart of [DDB80] 1.13 is located between the two ROSAT 
PSPC sources. Furthermore, PFJ93 report that this SNR is 
located 'within 19'' of a. brighter X-my source' which matches 
positionally with [SHP97] 211. These findings are confirmed by 
XMM-Newton and Chandra observations. The X-ray counter- 
part of [DDB80] 1.13 is source N^ 1 066 in the XMM LP-total 
catalogue (or [PFH2005] 354 or r3-69 in Kong et al. 2002b). 
The second source, which correlates with [SHP97] 211, is the 
XMM-Newton source N- 1 077, which has a "hard" spectrum 
and is ~ 6.7 times brighter than N^ 1 066. Hence, [SHP97] 211 
is a blend of the two XMM-Newton sources N- 1 066 and N- 
1 077. This also explains the pronounced variability between 
[SHL2001] 206 and N^ 1066 given in Table 17. Comparing 
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Table 17. Flux comparison of SNRs and SNR candidates from the XMM LP-total catalogue with counterparts classified as SNRs in ROSATmd Chandra catalogues 



XLPt SHP97 SHL2001 KGP2002+ WGK2004+ fvar svar reason why indicated vaiability is not reliable 

XID Flux with error in 10"^^ erg cm"^ s"^ 
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5.28 
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10.13 
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59.12+1.10 
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54.11 ±3.33 
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12.34 
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23.88 + 1.08 


53.50 ± 8.47 


79.39 ± 9.90 






3.32 


5.58 


1328 


9.25 + 0.74 




26.99 






1.00 


0.00 
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4.96 + 0.68 
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17.77 
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0.00 
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9.82 


4.07 


1535 


14.73 ± 1.31 


53.94 ±9.14 


34.41 ± 7.20 






3.66 


4.25 


1599 


16.08 ± 0.92 


54.39 ±10.03 


33.51 ± 6.97 






3.38 


3.80 


1637 


12.72 ± 1.33 




27.21 






1.00 


0.00 



ROSAT source is a blend of two XMM-Newton sources 
embedded in diffuse emission in central area of M 31 



see Sect. 10.2 



Notes: 

+ : KGP2002: Kong et al. (2002b), WGK2004: WilHams et al. (2004) 

ROSAT and Chandra count rates are converted to 0.2-4.5 keV fluxes, using WebPIMMS and assuming a foreground absorption of A^h = 6.6x10^^ cm"^ and a photon index of F = 1.7: 
ECFsHP97 = 2.229x10"^^ erg cm"2cts"\ ECFshl2001 = 2.249xl0"^^ergcm"2 cts"\ and ECFkgp2002 = 8.325xlO~^^ ergcm"^ cts"^ For WGK2004 the luminosity given in Table 2 of 
WGK2004 was converted to XID flux using FxioLergcm-^ s'^] =6.654x10"^^ xLwGK2004[10^^erg s"^]. 



Table 18. Flux comparison of SNRs and SNR candidates from the XMM LP-total catalogue which have counterparts in ROSAT, and/or Chandra catalogues that are not classified as SNRs 



N^ 


XLPt 


Chandra 


PFJ93 


SHP97 


SHL2001 


fvar 


svar 


remark* 






XID Flux with error in 10" 


erg cm"^ s~^ 










(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


294 


18.50 ±0.85 






53.27 ± 6.69 


46.78 ±7.87 


2.88 


5.16 




472 


3.15 ±0.69 








26.09 ± 6.07 


8.28 


3.76 


468 brt 


969 


53.51 ± 1.35 


84.51 ± 15.97+ 




34.55 ±6.91 


89.06 ± 11.92 


2.58 


3.96 




1079 


4.19 ±0.59 






20.06 ± 6.24 




4.79 


2.53 


brt 


1291 


14.55 ± 0.75 


16.04* 


>24.0 


35.22 ± 8.47 


40.93 ± 7.87 


2.81 


3.33 




1741 


4.12 ±0.65 


4.17^ 








1.01 




brt 


1793 


3.70 ± 0.52 






26.08 ± 6.46 




7.06 


3.46 


1 799 brt 



Notes: 

* : Source number (from XMM LP-total catalogue) of another (brighter) XMM LP-total source which correlate with the same ROSAT source as the XMM LP-total source given in Col. 1 ; 
brt: XMM LP-total flux is below the /?05Ar detection threshold (5.3x10"^^ erg cm"^ ^-i) 

ROSAT and Chandra count rates are converted to 0.2-4.5 keV fluxes, using WebPIMMS and assuming a foreground absorption of A^h = 6.6x10^^ cm"^ and a photon index of F = 1.7: 
ECFsHP97 = 2.229X lO-^^ erg cm-^ cts-^ , ECFshl20oi = 2.249x10-^4 erg cm-^ cts-^ , ECFhri = 6.001x10-^4 erg cm'^ cts'^ , * : ECFdkg2004 = 5.56x10-^2 erg cm'^ cts'^ . + : For WGK2004 
the luminosity given in Table 2 of WGK2004 was converted to XID flux using FxiD[erg cm"^ s-^] = 6.654x10"^^ xLwGK2004[10^^erg s"M- * : For VG2007 the luminosity given in Table 2 
of VG2007 was converted to XID flux using FxioLergcm^^ ^-ij =9.433x10-^^ xLvG2007[10^Srg s"^]. 



the Chandra detections of the SNR counterpart with the 
XMM-Newton flux gives a variabiHty factor of Fyar ~ 1.12. 

The distance between [SPH97] 203 and [DDB80] 1.13 is 
> 20''. [SPH97] 203 was reported only in the first ROSAT 
PSPC catalogue. It was not detected in the observations of the 
second ROSAT PSPC catalogue or in any XMM-Newton or 
Chandra observation of that region. Thus it seems very likely 
that [SPH97] 203 was either a transient source or a false detec- 
tion. In both cases [SPH97] 203 cannot be a SNR. As the field 
of [DDB80] 1.13 was observed many times with Chandra, and 
as Chandra has detected weak SNRs in the central part of M 3 1 
(Kong et al. 2003, and e. g. [DDB80] 1.13), Chandra should have 
detected X-ray emission corresponding to the ROSAT source 
[SPH97] 203, if it really belonged to a SNR. 

The remaining two ROSAT SNRs correlate with XMM- 
Newton sources, which were not classified as SNRs or SNR 
candidates. Source [SHP97] 258 correlates with source N- 1 337 
and has a 3(T positional error of 30''. From the improved spatial 
resolution of XMM-Newton the total positional error reduces to 
2 '.'3. Hence, we can see that the X-ray source belongs to a fore- 



ground star candidate (cf. Table 5) and not to the very nearby 
SNR. Source [SHL2001] 129 correlates with sources N^ 743 
and N- 761, which are classified as a GIC and a GIC candi- 
date, respectively. The SNR candidate listed as the counterpart 
of [SHL2001] 129 is located between these two XMM-Newton 
sources. In addition PFH2005 gives a third source which lies 
within the error circle of [SHL2001] 129 and which is classified 
as an AGN candidate. Thus it is very likely that [SHL2001] 129 
is a blend of these three XMM-Newton sources and that the cor- 
relation with the SNR candidate has to be considered as a chance 
coincidence. 

From the sources listed as SNRs in the diff'erent Chandra 
studies many are re-detected. Nevertheless two SNRs from 
Chandra were not detected in the XMM-Newton observations. 
Source nl-85 has been reported as spatially correlated with an 
optical SNR by Williams et al. (2004), but has also been classi- 
fied as a repeating transient source in the same paper. An XMM- 
Newton counterpart to nl-85 was detected neither in the study of 
PFH2005 nor in the XMM LP-total catalogue. The transient na- 
ture of this source is at odds with the SNR classification. Source 
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CX0M31 J004247.8+411556 (Kong et al. 2003), which corre- 
lates with the radio source [B90] 95, is located in the vicinity 
of two bright sources and close to the centre of M 31. Due to 
XMM-Newton 's larger point spread function this source cannot 
be resolved by XMM-Newton in this environment. The larger 
PSF of XMM-Newton is also the reason why source N- 1 050 
has a significant variability in Table 17, since this source is lo- 
cated within the central diff'use emission of M 3 1 . 

Finally, we wanted to determine whether any of the XMM 
LP-total SNRs and SNR candidates were previously observed, 
but not classified as SNRs. In total there are seven such sources. 

One of them (N- 1741) is classified as a SNR candidate 
based on its XMM-Newton hardness ratios, and correlates with 
the Chandra source nl-48 (DKG2004). The fluxes obtained 
with XMM-Newton and Chandra are in good agreement (see 
Table 18), but below the 7?05Ar detection threshold (5.3 xlO'^^ 
erg cm"^ s"^). 

For a further four sources, the corresponding sources were 
only detected previously with ROSAT. One of them (N- 
1 793=[SHP97] 347) also correlates with a radio source (source 
472 of Braun 1990) and is therefore identified as a SNR. The 
rather high flux variability between the ROSAT and XMM- 
Newton observations (see Table 18) can be attributed to source 
N^ 1 799, which is located within 19^9 of N^ 1 793. This sug- 
gests that [SHP97] 347 is a combination of both XMM-Newton 
sources, but as [SHP97] 347 was not detected in SHL2001, we 
cannot exclude a transient source or false detection as an expla- 
nation for the ROSAT source. Source N^ 472 (=[SHL2001] 84), 
source N^ 294 (=[SHP97] 53=[SHL2001] 56), and source N^ 
1 079 (=[SHP97] 212) are SNR candidates based on their hard- 
ness ratios. The pronounced flux variability of source N- 472 is 
due to source N^ 468, which is located within 18V5 of N^ 472 
and is ~8.6 times brighter than N- 472. The observed flux for 
source N- 1 079 was below the ROSAT threshold. Furthermore, 
the ROSAT source [SHP97] 212 was classified as a SNR, but did 
not appear in the SHL2001 catalogue. Hence ROSAT may have 
detected an unrelated transient instead. 

Sources corresponding to the remaining two XMM- 
Newton sources were detected with ROSAT and Chandra. 
Source N^ 969 was detected in both ROSAT PSPC surveys 
([SHP97] 185=[SHL2001] 186) and correlates with Chandra 
source sl-84 (WilHams et al. 2006b). We classify it as a SNR 
candidate due to its hardness ratios and X-ray spectrum (see 
XMMM31 J004239.8-h404318). Counterparts for source N^ 
1 291 were reported in the literature as [PFJ93] 84, [SHP97] 251, 
[SHL2001] 255, [VG2007] 261 and source 4 in Table 5 of Orio 
(2006). Based on the XMM-Newton hardness ratios and the cor- 
relation with radio source [B90] 166 (Braun 1990), we identified 
the source as a SNR. For sources N^ 294, N^ 969, and N^ 1 291 
the variability between diff'erent observations may not be real 
because of systematic cross-calibration uncertainties. Therefore, 
we keep the <SNR> and SNR classifications for these sources. 

10.2.2. The spatial distribution 

To examine the spatial distribution of the XMM-Newton SNRs 
and SNR candidates, we determined projected distances from 
the centre of M 31. The distribution of SNRs and SNR candi- 
dates (normalised per deg^) is shown in Fig. 24. It shows an en- 
hancement of sources around ~3 kpc, which corresponds to the 
SNR population in the 'inner spiral arms' of M 31. In addition, a 
second enhancement of sources around -10 kpc is detected; this 
corresponds to the well known dust ring or star formation ring 
in the disc of M 31 (Block et al. 2006). Only a few sources are 
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Figure 24. Projected radial distribution of SNRs and SNR candidates 
from the XMM LP-total catalogue. An enhancement in the source dis- 
tribution corresponding to the 10 kpc dust ring of M 31 is visible. 

located beyond this ring. Figure 25 shows the spatial distribution 
of the SNRs and SNR candidates from the XMM LP-total cata- 
logue plotted over the IRAS 60yum image (Wheelock et al. 1994). 
We see that most of the SNRs and SNR candidates are located on 
features that are visible in the IRAS image. This again demon- 
strates that SNRs and SNR candidates are coincident with the 
dust ring at ~ 10 kpc. In addition, the locations of star forming 
regions obtained from GALEX data (Kang et al. 2009, and pri- 
vate communication) are indicated in Fig. 25. We see that many 
of the SNRs and SNR candidates are located within or next to 
star forming regions in M 3 1 . 

10.3. X-ray binaries 

X-ray binaries consist of a compact object plus a companion star. 
The compact object can either be a white dwarf (these systems 
are a subclass of CVs), a neutron star (NS), or a black hole (BH). 
A common feature of all these systems is that a large amount of 
the emitted X-rays is produced due to the conversion of gravita- 
tional energy from the accreted matter into radiation by a mass- 
exchange from the companion star onto the compact object. 

X-ray binaries containing an NS or a BH are divided into two 
main classes, depending on the mass of the companion star: 

- Low mass X-ray binaries (LMXBs) contain companion stars 
of low mass (M ^ 1 M©) and late type (type A or later), 
and have a typical lifetime of ~10^"^ yr (Fabbiano 2006). 
LMXBs can be located in globular clusters. Mass transfer 
from the companion star into an accretion disc around the 
compact object occurs via Roche-lobe overflow. 

- High mass X-ray binaries (HMXBs) contain a massive O 
or B star companion (Mstar ~ lOM©, Verbunt & van den 
Heuvel 1995) and are short-lived with lifetimes of ~10^"^ yr 
(Fabbiano 2006). One has to distinguish between two main 
groups of HMXBs: super-giant and the Be/X-ray binaries. 
In these systems wind-driven accretion onto the compact ob- 
ject powers the X-ray emission. Mass-accretion via Roche- 
lobe overflow is less frequent in HMXBs, but is still known 
to occur in several bright systems (e. g. LMC X-4, SMC X- 
1, CenX-3). HMXBs are expected to be located in areas 
of relatively recent star formation, between 25-60 Myr ago 
(Antoniou et al. 2010). 
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Figure 26. Distribution of the source fluxes of XRBs and GIC sources 
in the 0.2-4.5 keV (XID) band. The diagram shows the number of iden- 
tified and classified XRBs and GlCs at each flux bin, plotted versus the 
flux. In addition, the individual distribution of (field) XRBs (in red) as 
well as GlCs (in green) are given. 



Figure 25. An IRAS GOjim image (Wheelock et al. 1994), which clearly 
shows the dust ring located at ~ 10 kpc, over-plotted with the location of 
SNRs and candidates (red dots) from the XMM LP-total catalogue. The 
coincidence between the SNRs and candidates and the structures of the 
image is visible. In addition the locations of star forming regions, which 
were obtained from GALEX data (Kang et al. 2009), are indicated by 
blue dots. Furthermore the two ellipses (green) at 3 and 10 kpc from the 
centre correspond to the enhancemnets of sources from Fig. 24. 



We should expect about 45 LMXBs in M 31, following a 
similar estimation as the one presented in Sect. 10.2. Here the 
number of LMXBs in the Galaxy was estimated from Grimm 
et al. (2002). In the XMM LP-total catalogue 88 sources are 
identified/classified as XRBs. This is not surprising as we may 
expect M 31 to have a higher fraction of XRBs than the Galaxy 
since it is an earlier type galaxy composed of a higher fraction 
of old stars. 

XRBs are the main contribution to the population of "hard" 
X-ray sources in M 3 1 . Despite some more or less reliable can- 
didates, not a single, definitely detected HMXB is known in 
M 31. The results of a new search for HMXB candidates are 
presented in Sect. 10.3.2. The LMXBs can be separated into two 
sub-classes: the field LMXBs (discussed in this section) and 
those located in globular clusters. Sources belonging to the latter 
sub-class are discussed in Sect. 10.4. 

The sources presented here are classified as XRBs, because 
they have HRs indicating a <hard> source and are either tran- 
sient or show a variability factor larger than ten (see Sect. 6). 

In total 10 sources are identified and 26 are classi- 
fied as XRBs by us, according to the classification crite- 
ria given in Table 6. Apart from source N- 57 (XMMM31 
J003833.2+402133, see below), the identified XRBs had been 
reported as X-ray binaries in the literature (see comment col- 
umn of Table 5). Figure 26 (red histogram) shows the flux dis- 



tribution of XRBs. We see that this class contains only rather 
bright sources. This is not surprising as the classification crite- 
rion for XRBs is based on their variability, which is more eas- 
ily detected for brighter sources (cf. Sect. 6). The XID fluxes 
range from 1.4x10-^'^ erg cm'^ s'^ (N^ 378) to 3.75x10-^^ 
erg cm"^ s"^ (N- 966), which correspond to luminosities from 
1.0x10^^ erg s"^ to 2.7x10^^ erg s"^. 

It is clear from Fig. 27, which shows the spatial distribution 
of the XRBs, that nearly all sources classified or identified as 
XRBs (yellow dots) are located in fields that were observed more 
than once (centre and southern part of the disc). This is partly a 
selection eflfect, caused by the fact that these particular fields 
were observed several times, thus allowing the determination of 
source variability. For sources located outside these fields, es- 
pecially the northern part of the disc, the transient nature must 
have been reported in the literature to mark them as an XRBs. 
The source density of LMXBs, which follows the overall stel- 
lar density, is higher in the centre than in the disc of M 3 1 . One 
would not expect HMXBs in the central region which is domi- 
nated by the bulge (old stellar population). From Fig. 28, which 
shows the spatial distribution of the XRBs over-plotted on an 
IRAS 60yum image (Wheelock et al. 1994), we see that only a 
few sources, classified or identified as XRBs, are located in the 
vicinity of star forming regions. 

References for the sources, selected from their temporal vari- 
ability, are given in Table 9. TPC06 report on four bright X- 
ray transients, which they detected in the observations of July 
2004 and suggested to be XRB candidates. We also found these 
sources and classified source N- 705 and identified sources N- 
985, 1 153, 1 177 as XRBs. One of the identified XRBs 
(N- 1 177) shows a very soft spectrum. Williams et al. (2005b) 
observed source N- 1 153 with Chandra and HST. From the lo- 
cation and X-ray spectrum they suggest it to be an LMXB. They 
propose that the optical counterpart of the X-ray source is a star 
within the X-ray error box , which shows an optical brightness 
change (in B) by ^1 mag. Source N- 985 was first detected in 
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Figure 27. The spatial distribution of XRBs and candidates from the 
XMM LP-total catalogue. The positions of the XRBs and candidates 
are marked with yellow dots; the two XRB candidates classified from 
their variability compared with ROSAT observations are marked with 
blue dots. An increase in the number density of sources in the central 
field is clearly visible. In addition the two new HMXB candidates pre- 
sented in Sect. 10.3.2 (red dots), and the three HMXB candidates of 
SBK2009 that satisfy our U-B/B-V selectrion criterion (green dots, see 
Sect. 10.3.2) are shown. XRBs which correlate with globular clusters 
are shown in Fig. 31. 



January 1979 by TF91 with the Einstein observatory. WGM06 
rediscovered it in Chandra observations from 2004. Their coor- 
dinated HST ACS imaging does not reveal any variable optical 
counterpart. From the X-ray spectrum and the lack of a bright 
star, WGM06 suggest that this source is an LMXB with a black 
hole primary. 

In the following subsections we discuss three transient XRBs 
in more detail. 

XMMM31 J003833.2+4021 33 (N^ 57) was first detected in the 
XMM-Newton observation from 02 January 2008 (s32) at an 
unabsorbed 0.2-lOkeV luminosity of ~ 2x10^^ erg s"^. From 
two observations, taken about 0.5 yr (s31) and 1.5 yr (s3) earlier, 
we derived upper limits for the fluxes, which were more than a 
factor of 100 below the values obtained in January 2008. 

The combined EPIC spectrum from observation s32 
(Fig. 29(a)) is best fitted with an absorbed disc blackbody plus 
power-law model, with A^h = 1.68^Q4gXlO^^ cm"^, temperature 
at the inner edge of the disc ^in = 0.462+0.01 3 keV and power- 
law index of 2.55^^ ^^. The contribution of the disc blackbody 
luminosity to the total luminosity is ~ 59 %. Formally accept- 
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Figure 28. The spatial distribution of XRBs and candidates from the 
XMM LP-total catalogue. Shown are the same sources as in Fig 27, but 
over-plotted on an IRAS 60yL/m image (Wheelock et al. 1994), which 
shows the dusty star forming region in M 31. In addition the locations 
of star forming regions, which were obtained from GALEX data (Kang 
et al. 2009), are indicated by cyan dots. 



able fits are also obtained from an absorbed disc blackbody and 
an absorbed bremsstrahlung model (see Table 19). 

We did not find any significant feature in a fast Fourier trans- 
formation (FFT) periodicity search. The combined EPIC light 
curve during observation s32 was consistent with a constant 
value. 

To identify possible optical counterparts we examined the 
LGGS images and the images taken with the XMM-Newton op- 
tical monitor during the X-ray observation (UVWl and UVW2 
filters). The absence of optical/UV counterparts and of variabil- 
ity on short timescales, as well as the spectral properties suggest 
that this source is a black hole LMXB in the steep power-law 
state (McClintock & Remiflard 2006). 

GX0M31 J004059.2+41 1551 : Galache et al. (2007) reported 
on the detection of a previously unseen X-ray source in a 5 ks 
Chandra ACIS-S observation from 05 July 2007. In an XMM- 
Newton ToO observation (snll, Stiele et al. 2007) taken about 
20 days after the Chandra detection, the source (N- 523) was 
still bright. The position agrees with that found by Chandra. We 
detected the source at an unabsorbed 0.2- lOkeV luminosity of 
-1.1x10^^ erg s-^ 

The combined EPIC spectrum (Fig. 29(b)) can be well 
fitted with an absorbed disc blackbody model with A^h = 
(2.00 ± 0.16)xl0^^ cm"^ and with a temperature at the inner 
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Table 20. The cumulative number of sources from the XMM LP-total 
catalogue that were classified as <hard>, and in addition do not corre- 
late with a source listed in the /?05Ar catalogues. Four diff'erent limiting 
fluxes are indicated. 



XMMM31 J003833.2+402133 



XMM-Newton, EPIC 



XID flux limit # of sources 



erg cm ^ s ^ 



5.5 E-15 
lE-14 
5E-14 
lE-13 



541 
242 
7 
1 



edge of the disc of fe^in = 0.538 ± 0.017 keV (Table 19). The 
spectral parameters and luminosity did not change significantly 
compared to the Chandra values of Galache et al. (2007). 

We did not find any significant feature in an FFT periodicity 
search. The combined EPIC light curve was consistent with a 
constant value. 

The examination of LGGS images and of images taken with 
the XMM-Newton optical monitor (UVWl and UVW2 filters) 
during the X-ray observation did not reveal any possible opti- 
cal/UV counterparts. 

The lack of bright optical counterparts and the X-ray param- 
eters (X-ray spectrum, lack of periodicity, transient nature, lumi- 
nosity) are consistent with this source being a black hole X-ray 
transient, as already mentioned in Galache et al. (2007). 

XMMU J004144.7+411110 (N^ 705) was detected by 
Trudolyubov et al. (2006) in XMM-Newton observations bl-b4 
(July 2004) at an unabsorbed luminosity of 3.1-4.4x10"^^ erg s"^ 
in the 0.3-7 keV band, using a DISKBB model. We detected the 
source in observation snll (25 July 2007) with an unabsorbed 
0.2-lOkeV luminosity of ~1.8xlO^^ erg s"^ using also a 
DISKBB model. 

In observation snll, the source was bright enough to al- 
low spectral analysis. The spectra can be well fitted with an 
absorbed power-law, disc blackbody or bremsstrahlung model 
(Table 19). The obtained spectral shapes (absorption and temper- 
ature as well as photon index) are in agreement with the values 
of Trudolyubov et al. (2006). 

An FFT periodicity search did not reveal any significant pe- 
riodicities in the 0.3 s to 2 000 s range. 

No optical counterparts were evident in the images taken 
with the XMM-Newton optical monitor UVWl and UVW2 dur- 
ing the snll observation, nor in the LGGS images. The lack of 
a bright optical counterpart and the X-ray parameters support 
that this source is a black hole X-ray transient, as classified by 
Trudolyubov et al. (2006). 
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Figure 29. EPIC spectra of the transient sources (a) 
XMMM31 J003833.2+402133, (b) CX0M31 J004059.2+411551 
and (c) XMMU J004144.7+411110. The histograms show the best-fit 
model: PL+DISCBB in (a), DISCBB in (b) and (c). 



10.3.1 . Sources from the XMM-LP total catalogue that were 
not detected by ROSAT 

To search for additional XRB candidates, we selected all sources 
from the XMM LP-total catalogue, that were classified as 
<hard> and which did not correlate with a source listed in the 
/?05Ar catalogues (PFJ93, SHP97 and SHL2001). The flux dis- 
tribution of the selected sources is shown in Fig. 30, and Table 20 
gives the number of sources brighter than the indicated flux limit. 

Possible, new XRB candidates are sources that have an XID 
flux that lies at least a factor of ten above the ROSAT detection 



threshold (5.3x10"^^ erg cm"^ s"^). These sources fulfil the vari- 
ability criterion used to classify XRBs (cf. Sect. 6). The XMM 
LP-total catalogue lists five sources without ROSAT counter- 
parts that have XID fluxes above 5.3x10"^"^ erg cm"^ s"^ . These 
are: 239, 365, 910, 1 164, and 1 553. Between 
the ROSAT and XMM-Newton observations more than ten years 
have elapsed. On this time scale AGN can also show strong vari- 
ability. To estimate the number of AGN among the five sources 
listed above, we investigated how many sources of the identi- 
fied and classified background objects from the XMM LP-total 
catalogue with an XID flux larger than 5.3x10"^"^ erg cm"^ s"^ 
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Table 19. Spectral parameters of the transient sources. 



M 31 field Model A^h k^T Vcos i* Photon ^1 Instrument 

(xlO^i cm-2) (keV) (km) Index (d.o.f) 

XMMM31 J003833.2+402133 (N^ 57) 



s32 PL+DISCBB ^-^^^Ti, 0.462 ±0.013 106^?^ 2.55+?g 173.89(145) 2.04 PN+M1+M2 

s32 DISCBB 1.06 ±0.06 0.511 ±0.009 95 ± 4 270.01(147) 1.46 PN+M1+M2 

s32 BREMSS 1.91 ±0.07 1.082+0029 208.65(147) 2.12 PN+M1+M2 











CX0M31 J004059.2+411551 (N 


^523) 








snll 
snll 


DISCBB 
BREMSS 


2.00 ±0.16 
3.13 ±0.19 


0.538 ±0.017 75 ±6 

1 007+0060 




97.70(79) 
93.17(79) 


1.12 
1.72 


PN+M1+M2 
PN+M1+M2 










XMMU J004 144.7+41 1 1 10 (N^ 


705) 









snll DISCBB 2.32+Jo^ 0.586!OJoo 26!^^ 29.74(23) 0.18 PN+M1+M2 

snll BREMSS 3.72+}-JJ 1.216+[J-3g 29.48(23) 0.29 PN+M1+M2 

snll PL 6.17+}]2 3.23:^[J-^6 31.57(23) 1.12 PN+M1+M2 



Notes: 

* : efi'ective inner disc radius, where / is the inclination angle of the disc 

^ : unabsorbed luminosity in the 0.2- 10.0 keV energy range in units of 10^^ erg s~^ 
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Figure 30. Distribution of the source fluxes in the 0.2-4.5keV (XID) 
band. The diagram shows the number of sources from the XMM LP- 
total catalogue that were classified as <hard>, and in addition do not 
correlate with a source listed in the ROSAT catalogues at each flux bin 
plotted versus the flux, using logarithmic scales. 



were not detected by ROSAT. The result is that ROSAT detected 
all background sources with an XID flux larger than 5.3x10"^^ 
erg cm"^ s"^ that are listed in the XMM LP-total catalogue. 
Thus, the probability that any of the five sources listed above 
is a background object is very small, in particular if the source 
is located within the D25 ellipse of M 31. Therefore, the two 
sources located within the D25 ellipse are listed in the XMM LP- 
total catalogue as XRB candidates, while the remaining three 
sources, which are located outside the D25 ellipse, are classified 
as <hard>. All five sources are marked in the comment column 
of Table 5 with 'XRB cand. from 7?0^Ar corr.' . 



10.3.2. Detection of high mass X-ray binaries 

As already mentioned, until now not a single secure HMXB in 
M 31 has been confirmed. The reason for this is that the de- 
tection of HMXBs in M 31 is diflScult. Colbert et al. (2004) 
showed that the hardness ratio method is very ineflftcient in se- 
lecting HMXBs in spiral galaxies. The selection process is com- 
plicated by the fact, that the spectral properties of BH HMXBs, 
which have power-law spectra with indices of ~1 2 are sim- 
ilar to LMXBs and AGN. Therefore the region in the HR dia- 
grams where BH HMXB are located is contaminated by other 
hard sources (LMXBs, AGN, and Crab like SNRs). For the NS 
HMXBs, which have power-law indices of ~1, and thus should 
be easier to select, the uncertainties in the hardness ratios lead at 
best to an overlap - in the worst case to a fusion - with the area 
occupied by other hard sources (Colbert et al. 2004). 

Based on the spectral analysis of individual sources of 
M 31, SBK2009 identified 18 HMXB candidates with power- 
law indices between 0.8 and 1.2. One of these sources 
([SBK2009] 123) correlates with a globular cluster, and hence 
it is rather an LMXB in a very hard state rather than an HMXB 
(c/ Trudolyubov & Priedhorsky 2004). Four of their sources 
([SBK2009] 34, 106, 149, and 295) do not have counterparts in 
the XMM LP-total catalogue. 

Eger (2008) developed a selection algorithm for HMXBs in 
the SMC, which also uses properties of the optical companion. 
X-ray sources were selected as HMXB candidates if they had 
HR2+EHR2>0.1 as well as an optical counterpart within 2V5 
of the X-ray source, with -0.5 <B-V< 0.5 mag, -1.5 <U-B< 
-0.2 mag and V<17mag. 

We tried to transfer this SMC selection algorithm to M 3 1 
sources. In doing so, we encountered two problems: The first 
problem is that the region of the U-B/B-V diagram is also popu- 
lated by globular clusters (LMXB candidates) in M 3 1 . The sec- 
ond problem is that due to the much larger distance to M 31, 
the range of detected V magnitudes of HMXBs in the SMC 
of ~13<V<17mag translates to a ~19<V<23mag criterion for 
M 3 1 . Thus the V magnitude of optical counterparts of possible 
HMXB candidates lies in the same range as the optical counter- 
parts of AGN. Therefore the V mag criterion, which provided 
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most of the discriminatory power in the case of the SMC, fails 
totally in the case of M 3 1 . 

A few of the sources selected from the optical colour-colour 
diagram and HR diagrams are bright enough to allow the cre- 
ation of X-ray spectra. That way two additional (/. e. not given 
in SBK2009) HMXB candidates were found. 

In addition, we determined the reddening free Q parameter: 

2 = (U - B) - 0.72(B - V) (8) 

(for definition see e. g. Cox 2001) which allowed us to keep 
only the intrinsically bluest stars, using Q < -0.4 (0-type stars 
typically have Q< -0.9, while -0.4 corresponds to a B5 dwarf or 
giant or an AO supergiant, Massey et al. 2007). U-B and B-V 
were taken from the LGGS catalogue. 

XMMM31 J004557.0+414830 (N^ 1716) has an USNO-Bl 
(R2= 18.72 mag), a 2MASS and an LGGS (V=20.02mag; Q= 
-0.44) counterpart. The EPIC spectrum is best fitted (x^^^ = 
0.93) by an absorbed power-law with A^h= 7.4+^-^x10^^ cm"^ 
and photon index F = 1.2 ± 0.4. The absorption corrected X-ray 
luminosity in the 0.2-10 keV band is ~7.1xl0^^ erg s"^ 

XMMM31 J004506.4+420615 (N^ 1579) has an USNO-Bl 
(B2=20.87mag), a 2MASS and an LGGS (V=20.77mag; 
Q= -1.04) counterpart. The EPIC PN spectrum is best fitted 
(X^^^ = 1.6) by an absorbed power-law with A^h= 0.48^^ ^x10^^ 
cm"^ and photon index F = l-0^0 5' absorption corrected 
X-ray luminosity in the 0.2-10 keV band is -8.6x10^^ erg s"^ . 

To strengthen these classifications spectroscopic optical 
follow-up observations of the optical counterparts are needed. 
An FFT periodicity search did not reveal any significant period- 
icities for either of the two sources and the light curves do not 
show eclipses. 

From the sources reported as HMXB candidates in 
SBK2009, three sources ([SBK2009] 21, 236, and 256) are lo- 
cated in the region of the U-B/B-V diagram, that we used. 
Another three sources ([SBK2009] 123, 172, and 226) are lo- 
cated outside that region. The remaining sources of SBK2009 
have either no counterparts with a U-B colour entry in the LGGS 
catalogue ([SBK2009] 99, 234, 294, and 302) or have no optical 
counterpart from the LGGS catalogue at all ([SBK2009] 9, 160, 
197, and 305). The reddening free Q parameter for the SBK2009 
sources that have counterparts in the LGGS catalogue are given 
in Table 21. 

10.4. Globular cluster sources 

A significant number of the luminous X-ray sources in the 
Galaxy and in M 31 are found in globular clusters. X-ray 
sources corresponding to globular clusters are identified by 
cross-correlating with globular cluster catalogues (see Sect. 3.8). 
Therefore changes between the XMM LP-total catalogue and the 
catalogue of PFH2005 in the classification of sources related to 
globular clusters are based on the availability of and modifica- 
tions in recent globular cluster catalogues. 

In total 52 sources of the XMM LP-total catalogue corre- 
late with (possible) globular clusters. Of these sources 36 are 
identified as GlCs because their optical counterparts are listed as 
globular clusters in the catalogues given in Sect. 3.8, while the 
remaining 16 sources are only listed as globular cluster candi- 
dates. 



Table 21. Reddening free Q parameter for HMXB candidates of 
SBK2009. 







[SBK2009] 


LGGS counterpart 


Q 


312 


21 


J004001. 50+403248.0 


-0.34 


1668 


236 


J004538.23+421236.0 


-0.77 


1724 


256 


J004558.98+420426.5 


-0.81 


1109 


123 


J00430L51+413017.5 


-hi. 77 


1436 


172 


J004420.98+413546.7 


-0.65 






J00442L01+413544.3* 


-0.29 


1630 


226 


J004526.68+41563L5 


-0.92 






J004526.58+415633.1* 


-0.72 



Notes: 

* : counterparts listed in SBK2009 



The range of source XID fluxes goes from 3.1x10 
erg cm-2 s'^ (N^ 924) to 2.7x10-^2 ^^^-2 ^-i ^^o ^ 957^^ 

in luminosity from 2.3x10^^ erg s"^ to 2.0x10^^ erg s"^ (Fig. 26; 
green histogram). Compared to the fluxes found for the XRBs 
discussed in Sect. 10.3, 14 sources that correlate with GlCs have 
fluxes below the lowest flux found for field XRBs. The reason 
for this finding is that the classification of field XRBs is based 
on the variable or transient nature of the sources, which can only 
be to detected for brighter sources (cf. Sect. 6) and not just by 
positional coincidence that is also possible for faint sources. 

Figure 31 shows the spatial distribution of the GIC sources. 
X-ray sources correlating with GlCs follow the distribution of 
the optical GlCs, which are also concentrated towards the central 
region of M 3 1 . 

The three brightest globular cluster sources, 
which are located in the northen disc of M 31, 
are 1057 (XMMM31 J004252.0-h413109), 

694 (XMMM31 J004 143. l-h4 13420), and 1692 
(XMMM31 J004545.8-h413941). They are aU brighter than 
8.4x10^^ erg s"^ Source N- 694 was classified as a black 
hole candidate, due to its variability observed at such high 
luminosities. A detailed discussion of the three sources is given 
in Barnard et al. (2008). 

XMMM31 J004303.2-h412121 (N^ 1 118) was identified as 
a foreground star in PFH2005, based on the classification in the 
"Revised Bologna Catalogue" (Galleti et al. 2004). Galleti et al. 
(2004) took the classification from Dubath & Grillmair (1997), 
which is based on the velocity dispersion of that source. Recent 
'HST images unambiguously reveal that this [B147] is a well re- 
solved star cluster, as recently pointed out also by Barmby et al. 
(lOOiy (Galleti et al. 2007). That is why source 1 1 18 is now 
identified as an XRB located in globular cluster B 147. 

10.4.1 . Integrated optical properties of the globular clusters 
in which the X ray sources are located 

For each X-ray source which correlates with a globular cluster or 
globular cluster candidate in the optical, we investigated its in- 
tegrated V-I colour and derived age estimates. Table 22 lists the 
name of the optical counterpart, its classification according to 
RBC V.4 (Galleti et al. 2009), the age classification of Caldwefl 
et al. (2009), the V magnitude and V-I colour given in RBC V.4, 
the dereddened V-I colour, and the age estimate derived by our- 
selves. 

The integrated V-I colours of the clusters can be found in 
RBC V.4 and can be used to provide estimates of the ages of the 
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Figure 31. The spatial distribution of X-ray sources correlating with 
GlCs and candidates from the XMM LP-total catalogue (yellow dots). 
An enhancement of sources towards the central region of M 3 1 is clearly 
visible. 



clusters, in conjunction with reddening values. We have adopted 
a reddening of E(B-V)= 0. 10 ± 0.03, which is the average of the 
reddenings of all M 31 clusters in Rich et al. (2005). Using these 
values, we have derived (V-I)o for our clusters. In most cases 
(V-I)o > 1.0 suggesting clusters older than ^2 Gyr according 
to Sarajedini et al. (2007). The histogram in Fig. 32 shows the 
distribution of (V-I)o for our clusters, with the approximate age- 
ranges marked. 

In general there is good agreement between the Caldwell 
et al. (2009) and our age estimates. This result indicates that the 
great majority of the objects are indeed old globular clusters. 

Figure 33 shows the distribution of (V-I)o for all confirmed 
and candidate globular clusters, listed in the RBC V.4, which are 
located in the XMM LP-total field, and which have V as well 
as I magnitudes given. A comparison with Fig. 32 again reveals 
that mainly counterparts of old globular clusters (age ^2 Gyr) 
are detected in X-rays. 

10.4.2. Comparing GIG and candidates in XMM-Newton, 
Chandra and ROSAT catalogues 

The combined ROSAT PSPC catalogue (SHP97 and SHL2001) 
contains 33 sources classified as globular cluster counterparts. 
Of these sources one is located outside the field observed with 
XMM-Newton. Another two sources do not have counterparts in 
the XMM LP-total catalogue. The first one is [SHL2001] 232, 
which is not visible in any XMM-Newton observation taken be- 
fore December 2006 as was already reported in Trudolyubov & 



Table 22. Integrated V-I colours, dereddened V-I integrated colours, and age estimates of 
the globular clusters in which the X ray sources are located. 





Name 


Class 


Age 


mag 


VI 


VIO 


t 

Age' 


BOOS 


confirmed 


old 


15.69 


1.29 


1.15 


old 


SK055B 


candidate 




18.991 


0.388 


0.248 




B024 


confirmed 


old 


16.8 


1.15 


1.01 


old 


SKIOOC 


candidate 


na 


18.218 


1.181 


1.041 


old 


B045 


confirmed 


old 


15.78 


1.27 


1.13 


old 


B050 


confirmed 


old 


16.84 


1.18 


1.04 


old 


B055 


confirmed 


old 


16.67 


1.68 


1.54 


old 


B058 


confirmed 


old:: 


14.97 


1.1 


0.96 


old-inter 


MITA140 


confirmed 


old 


17 


9999 






B078 


confirmed 


old 


17.42 


1.62 


1.48 


old 


B082 


confirmed 


old 


15.54 


1.91 


1.77 


old 


B086 


confirmed 


old 


15.18 


1.26 


1.12 


old 


SK050A 


confirmed 




18.04 


1.079 


0.939 


old-inter 


B094 


confirmed 


old 


15.55 


1.26 


1.12 


old 


B096 


confirmed 


old 


16.61 


1.48 


1.34 


old 


B098 


confirmed 


old 


16.21 


1.13 


0.99 


old-inter 


B107 


confirmed 


old 


15.94 


1.28 


1.14 


old 


BllO 


confirmed 


old 


15.28 


1.28 


1.14 


old 


B117 


confirmed 


old:: 


16.34 


1 


0.86 


inter 


B116 


confirmed 


old 


16.79 


1.86 


1.72 


old 


B123 


confirmed 


old 


17.416 


1.29 


1.15 


old 


B124 


confirmed 


old 


14.777 


1.147 


1.007 


old 


B128 


confirmed 


old:: 


16.88 


1.12 


0.98 


old-inter 


B135 


confirmed 


old 


16.04 


1.22 


1.08 


old 


B143 


confirmed 


old 


16 


1.22 


1.08 


old 


B144 


confirmed 


old:: 


15.88 


0.59 


0.45 


young 


B091D 


confirmed 


old 


15.44 


9999 






B146 


confirmed 


old:: 


16.95 


1.09 


0.95 


interm 


B147 


confirmed 


old 


15.8 


1.27 


1.13 


old 


B148 


confirmed 


old 


16.05 


1.17 


1.03 


old 


B150 


confirmed 


old 


16.8 


1.28 


1.14 


old 


B153 


confirmed 


old 


16.24 


1.3 


1.16 


old 


B158 


confirmed 


old 


14.7 


1.15 


1.01 


old 


B159 


confirmed 


old 


17.2 


1.41 


1.27 


old 


B161 


confirmed 


old 


16.33 


1.1 


0.96 


old-inter 


B182 


confirmed 


old 


15.43 


1.29 


1.15 


old 


B185 


confirmed 


old 


15.54 


1.18 


1.04 


old 


B193 


confirmed 


old 


15.33 


1.28 


1.14 


old 


SK132C 


candidate 




18.342 


1.84 


1.7 


old 


B204 


confirmed 


old 


15.75 


1.17 


1.03 


old 


B225 


confirmed 


old 


14.15 


1.39 


1.25 


old 


B375 


confirmed 


old 


17.61:: 


1.02 


0.88 


interm 


B386 


confirmed 


old 


15.547 


1.154 


1.014 


old 



Notes: 

* : classification as confirmed or otherwise comes from the revised Bologna catalogue 
(December 2009, Version 4) http : //www . bo . astro . it/M3 l/RBC_Phot07_V4 . tab 

: age comes from Caldwell et al. (2009) 
V, and V-I are integrated colours that come from the revised Bologna catalogue (December 
2009, Version 4) http : //www . bo . astro . i t/M3 l/RBC_Phot07_V4 . tab 
(V-I)o is the dereddened V-I integrated colour, assuming E(B-V)=0. 10-1-0.03, which is the 
average of the reddenings of all M 31 clusters in Rich et al. (2005). (this E(B-V) corre- 
sponds to E(V-I)=0.14) 

^ : This dereddened colour (V-I)o is used to estimate the age on the basis of the plots (V-I)o 
versus log Age from Sarajedini et al. (2007). 

Priedhorsky (2004). The second source ([SHL2001] 231) cor- 
relates with B 164 which is identified as a globular cluster in 
RBC V3.5. In addition [SHL2001] 231 is fisted in PFH2005 as 
the counterpart of the source [PFH2005] 423. Due to the im- 
proved positional accuracy of the X-ray source in the XMM- 
Newton observations, PFH2005 rejected the correlation with 
B 164 and instead classified [PFH2005] 423 as a foreground star 
candidate. 

Three ROSAT GIC candidates have more than one counter- 
part in the XMM LP-total catalogue. [SHL2001] 249 correlates 
with sources N- 1 262 and N- 1 267, where the latter is the X-ray 
counterpart of the globular cluster B 185. [SHL2001] 254 corre- 
lates with sources N- 1 289 and N- 1 293, where the former is 
the X-ray counterpart of the globular cluster candidate mita 311 
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Figure 32. The distribution of (V-I)o for globular clusters (and can- 
didates), with the approximate age-ranges marked, hosting XMM LP- 
total X-ray sources. 




Figure 33. The distribution of (V-I)o for globular clusters (and candi- 
dates) from the RBC V4, located in the XMM LP-total field. 



(Magnier 1993). [SHL2001] 258 has a Icr positional error of 
48'' and thus correlates with sources 1 297, 1 305 and 
1 357.^"^ The brightest of these three sources (N- 1 305), which 
is actually located closest to the RO SAT position, correlates with 
the globular cluster candidate SK 132C (RBC V3.5). 

Table 23 gives the variability factors (Cols. 6, 8) and sig- 
nificance of variability (7, 9) for sources classified as GIC can- 
didates in the ROSAT FSFC surveys. For most sources only low 
variability is detected. The two sources with the highest variabil- 



24 In addition [SHL2001] 258 correlates with 1 275, 1 289 and 
N- 1 293. However these sources have each an additional /?0*SAr coun- 
terpart. 



ity factors found (N^ 1262, 1293) belong to ROSAT sources 
with more than one XMM-Newton counterpart. In these cases 
the XMM-Newton sources that correlate with the same ROSAT 
source and the optical globular cluster source show much weaker 
variability. Interestingly, a few sources show low, but very signif- 
icant variability. Among these sources is the Z- source identified 
in Barnard et al. (2003, N- 966) and two of the sources discussed 
in Barnard et al. (2008, N^ 1 057, N^ 1 692). 

The 18 X-ray sources correlating with globular clusters 
which were found in the ROSAT HRl observations (PFJ93) were 
all re-detected in the XMM-Newton data. 

From the numerous studies of X-ray globular cluster coun- 
terparts in M 31 based on Chandra observations (Kong et al. 
2002b; Di Stefano et al. 2002; Williams et al. 2004; Trudolyubov 
& Priedhorsky 2004; Voss & Gilfanov 2007), only eight sources 
are undetected in the present study. One of them ([TP2004] 1) 
is located far outside the field of M 31 covered by the Deep 
XMM-Newton Survey. The transient nature of [TP2004] 35, and 
the fact that it is not observed in any XMM-Newton observa- 
tion taken before December 2006 was mentioned in Sect. 7.2. 
The six remaining sources (r2-15, r3-51, r3-71, [VG2007] 58, 
[VG2007] 65, [VG2007] 82) are located in the central area of 
M 31 and are also not reported in PFH2005. Figure 34 shows 
the position of these six sources (in red) and the sources of 
the XMM LP-total catalogue (in yellow). If the brightness of 
the six sources had not changed between the Chandra and 
XMM-Newton observations, they would be in principle bright 
enough to be detected by XMM-Newton in the merged ob- 
servations of the central field, which have in total an expo- 
sure > lOOks. Two sources (r2-15 and [VG2007] 65) are lo- 
cated next to sources detected by XMM-Newton. Source r2- 
15 is located within 13'' of N^ 1012 and within 17'' of N^ 
1017 and has - in the Chandra observation - a similar lumi- 
nosity to both XMM-Newton sources. The distance between 
N^ 1012 and N^ 1017 is 17" and within 20" of N^ 1012, 
XMM-Newton detected source N- 1 006, which is about a fac- 
tor 4.6 fainter than N- 1 012. Therefore, when in a bright state, 
source r2-15 should be detectable with XMM-Newton. Source 
[VG2007] 65 is located within 17"of N^ 1 100, which is at least 
3.5 times brighter than [VG2007] 65. This may complicate the 
detection of [VG2007] 65 with XMM-Newton. The variability 
of [VG2007] 58, [VG2007] 65, and [VG2007] 82 is supported 
by the fact that these three sources were not detected in any 
Chandra study, prior to Voss & Gilfanov (2007). Hence, these 
six sources are likely to be at least highly variable or even tran- 
sient. 

Several sources identified with globular clusters in previ- 
ous studies have counterparts in the XMM LP-total catalogue 
but are not classified as GIC sources by us. Source N- 403 
([SHL2001] 74) correlates with B 007, which is now identified 
as a background galaxy (Caldwell et al. 2009; Kim et al. 2007, 
RBC V3.5). Sources N^ 793 ([SHL2001] 136, sl-12) and N^ 
796 (si- 11) are the X-ray counterparts of B 042D and B 044D, 
respectively, which are also suggested as background objects 
by Caldwell et al. (2009). Source N^ 948 (sl-83) correlates 
with B 063D, which is listed as a globular cluster candidate in 
RBC V3.5, but might be a foreground star (Caldwell et al. 2009). 
Due to this ambiguity in classification we classified the source as 
<hard>. Source N^ 966 correlates with [SHL2001] 184, which 
was classified as the counterpart of the globular cluster NB 21 
(RBC V3.5) in the ROSATFSFC survey (SHL2001). In addition, 
source N- 966 also correlates with the Chandra source r2-26 
(Kong et al. 2002b). Due to the much better spatial resolution of 
Chandra compared to ROSAT, Kong et al. (2002b) showed that 
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Figure 34. Image of the central field of M 3 1 
over-plotted with the positions of six possible 
transient sources (red) and the sources of the 
XMM LP-total catalogue. Sources r2-15 and 
r3-71 are listed as sources #17 and #28, respec- 
tively, in Di Stefano et al. (2002). The three 
"red" sources that are only marked with a num- 
ber (#58, #65, #82) are taken from Voss & 
Gilfanov (2007). 



source r2-26 does not correlate with the globular cluster NB 21. 
Barnard et al. (2003) identified this source as the first Z- source 
in M 31. The nature of source 1 078 is unclear as RBC V3.5 
reported that source to be a foreground star, while Caldwell et al. 
(2009) classified it as an old globular cluster. Due to this ambigu- 
ity in the classification and due to the fact that source N- 1 078 is 
resolved into two Chandra sources (r2-9, r2-10), we decided to 
classify the source as <hard>. Due to the transient nature (Kong 
et al. 2002b; Williams et al. 2006b) and the ambiguous classifi- 
cations reported by RBC V3.5 (GIC) and Caldwell et al. (2009, 
HII region), we adopt the classification of PFH2005 (<XRB>) 
for source N- 1 152. SBK2009 classified the source correlating 
with source N- 1 293 as a globular cluster candidate. We are not 
able to confirm this classification, as none of the globular clus- 
ter catalogues used, contains an entry at the position of source 
N- 1 293. Instead we found a radio counterpart in the catalogues 
of Gelfand et al. (2005), Braun (1990) and NVSS. We therefore 
classified the source as an AGN candidate, as was also done in 
PFH2005. 

For source 1 449 ([SHL2001] 289) the situation is more 
complicated. SHL2001 report [MA94a] 380 as the globular clus- 
ter correlating with this X-ray source. Based on the same refer- 
ence. Fan et al. (2005) included the optical source in their sta- 
tistical study of globular cluster candidates. However, the paper 
with the acronym [MA94a] is not available. An intensive liter- 
ature search of the papers by Magnier did not reveal any work 
relating to globular clusters in M 31, apart from Magnier (1993) 
which is cited in Fan et al. (2005) as "MIT". In addition the 
source is not included in any other globular cluster catalogues 
listed in Sect. 3.8. In the X-ray studies of Williams et al. (2004) 
and PFH2005 and in Magnier et al. (1992) the source is classi- 
fied as a foreground star (candidate). Hence, we also classified 



source N- 1 449 as a foreground star candidate, but suggest opti- 
cal follow-up observations of the source to clarify its true nature. 

A similar case is source 422 ([SHL2001] 76), which is 
classified as a globular cluster by SHL2001, based on a corre- 
lation with [MA94a] 16. Here again the source is not listed in 
any of the globular cluster catalogues used. We found one cor- 
relation of source N- 422 with an object in the USNO-Bl cat- 
alogue, which has no B2 and R2 magnitude. Two faint sources 
(V> 22.5 mag) of the LGGS catalogue are located within the X- 
ray positional error circle. Thus source N- 422 is classified as 
<hard>. While RBC V3.5 classified the optical counterpart of 
source 1 634 ([SHL2001] 316) as a globular cluster candi- 
date, Caldwell et al. (2009) regard SK 182C as being a source 
of unknown nature. Therefore we decided to classify source N- 
1634 as <hard>. 



11. Conclusions 

This paper presents the analysis of a large and deep XMM- 
Newton survey of the bright Local Group S A(s)b galaxy M 3 1 . 
The survey observations were taken between June 2006 and 
February 2008. Together with re-analysed archival observations, 
they provide for the first time full coverage of the M3 1 D25 el- 
lipse down to a 0.2-4.5 keV luminosity of ~10^^ erg s"^ 

The analysis of combined and individual observations al- 
lowed the study of faint persistent sources as well as brighter 
variable sources. 

The source catalogue of the Large XMM-Newton Survey 
of M 31 contains 1 897 sources in total, of which 914 sources 
were detected for the first time in X-rays. The XID source lumi- 
nosities range from -4.4x10^"^ erg s"^ to 2.7x10^^ erg s"^ The 
previously found diflTerences in the spatial distribution of bright 
(^10^^ erg s"^) sources between the northern and southern disc 
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Table 23. Variability between XMM LP-total and /^OS'Ar observations for sources classified as GIC candidates in the ROSAT FSFC surveys 
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383 


73 


68 


1.45E-12 


1.05E-14 


1.26 


11.58 


1.23 


9.04 


GIC 


* 


* 


403 


79 


74 


2.55E-14 


2.36E-15 


2.61 


2.09 


7.57 


4.92 


Gal 
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76 


2.01E-14 


1.43E-15 










<hard> 






694 


122 


113 


1.52E-12 


1.04E-14 


1.38 


10.56 


1.26 


7.46 


GIC 


* 


* 


793 


138 


136 


4.76E-14 


2.33E-15 


1.10 


0.55 


1.02 


0.14 


<Gal> 


* 




841 


150 


147 


1.40E-12 


1.61E-14 


1.77 


19.81 


2.11 


26.37 


GIC 


* 
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4.21E-13 


3.70E-15 


1.01 


0.14 


3.54 


48.75 


GIC 


* 
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1.56E-14 


1.60E-15 


1.70 


1.54 
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1.47E-13 


2.07E-15 
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3.44 


GIC 
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3.67E-14 


1.64E-15 


3.03 
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3.24E-13 


7.46E-15 


2.43 


12.78 


2.29 


12.16 
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* 




966 


184 
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1.00 


0.20 


2.31 
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2.67E-12 


2.05E-14 
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26.26 


1.79 


28.48 
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* 


* 
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3.23E-13 


2.93E-15 


1.06 


1.04 


5.51 


60.80 
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1 109 
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3.25E-13 


9.08E-15 


1.91 


9.70 


2.10 


10.72 
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* 
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1 118 


222 


216 


1.16E-13 


2.03E-15 


1.46 


3.63 


1.89 


7.46 


GIC 
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2.48E-13 


2.72E-15 


2.08 


12.02 


7.03 


73.27 
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7.59E-13 


4.48E-15 


1.06 


1.68 


1.08 


2.75 


GIC 
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4.68E-13 


4.93E-15 
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247 
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2.94E-14 


3.10E-15 


14.11 


18.40 


14.64 


20.24 
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247 
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4.57E-15 


1.16 


3.04 


1.11 


2.44 
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250 
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1.91E-15 
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8.53 
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<AGN> 






1296 
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1.58E-15 
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9.38 


1.25 


1.17 
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252 
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1.69E-14 
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1634 
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316 


7.70E-14 


2.91E-15 


3.14 


5.60 


1.89 


4.33 


<hard> 


* 
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1.15E-12 


2.00E-14 
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45.59 
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38.63 
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8.72E-13 


9.17E-15 


1.32 
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1.03 
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Notes: 

* : SI: SHP97, SII: SHL2001 

: XID Flux and error in erg cm"^ s"^ 

^ : Variability factor and significance of variability, respectively, for comparisons of XMM-Newton XID fluxes to ROSAT ^uxes listed in SPH97 and SHL2001, respectively. 

* : An asterisk indicates that the XID flux is larger than the corresponding ROSAT fiux. ROSAT count rates are converted to 0.2-4.5 keV fluxes, using WebPIMMS and assuming a foreground 
absorption of A^h =6.6x10^^ cm"^ and a photon index ofr=1.7: ECFshp97 =2.229x10"^^ erg cm"^ cts"^ and ECFshl2001 =2.249x10"^^ erg cm"^ cts"^ 



could not be confirmed. The identification and classification of 
the sources was based on properties in the X-ray wavelength 
regime: hardness ratios, extent and temporal variability. In ad- 
dition, information obtained from cross correlations with M 3 1 
catalogues in the radio, infra-red, optical and X-ray wavelength 
regimes were used. 

The source catalogue contains 12 sources with spatial extent 
between 6V2 and 23 VO. From spectral investigation and compar- 
ison with optical images, five sources were classified as galaxy 
cluster candidates. 

317 out of 1 407 examined sources showed long term vari- 
ability with a significance >3(T between the XMM-Newton ob- 
servations. These include 173 sources in the disc that were 
not covered in the study of the central field (SPH2008). Three 
sources located in the outskirts of the central field could not have 
been detected as variable in the study presented in SPH2008, as 
they only showed variability with a significance >3(T between 
the archival and the "Large Project" observations. For 69 sources 
the flux varied by more than a factor of five between XMM- 
Newton observations; ten of these varied by a factor >100. 

Discrepancies in source detection between the Large XMM- 
Newton Survey catalogue and previous XMM-Newton cata- 
logues could be explained by difl'erent search strategies, and 
difl'erences in the processing of the data, in the parameter set- 
tings of the detection runs and in the software versions used. 
Correlations with previous Chandra studies showed that those 



sources not detected in this study are strongly time variable, tran- 
sient, or unresolved. This is particularly true for sources located 
close to the centre of M 31, where Chandra's higher spatial res- 
olution resolves more sources. Some of the undetected sources 
from previous ROSAT studies were located outside the field cov- 
ered with XMM-Newton. However, there were several sources 
detected by ROSAT thai had a ROSAT detection likelihood larger 
than 15. If these sources were still in a bright state they should 
have been detected with XMM-Newton. Thus, the fact that these 
sources are not detected with XMM-Newton implies that they 
are transient or at least highly variable sources. On the other 
hand 242 <hard> XMM-Newton sources were found with XID 
fluxes larger than 10"^^ erg cm"^ s"\ which were not detected 
with ROSAT 

To study the properties of the difl'erent source populations 
of M 3 1 , it was necessary to separate foreground stars (40 plus 
223 candidates) and background sources (11 AGN and 49 candi- 
dates, 4 galaxies and 19 candidates, 1 galaxy cluster and 5 can- 
didates) from the sources of M 3 1 . 1 247 sources could only be 
classified as <hard>, while 123 sources remained without iden- 
tification or classification. The majority (about two-thirds, see 
Stiele et al. 2011 in preparation) of sources classified as <hard> 
are expected to be background objects, especially AGN. 

The catalogue of the Large XMM-Newton survey of M 3 1 
contains 30 SSS candidates, with unabsorbed 0.2-l.OkeV lumi- 
nosities between 2.4x10^^ erg s"^ and 2.8x10^^ erg s"^ SSSs 
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are concentrated to the centre of M 3 1 , which can be explained 
by their correlation with optical novae, and by the overall spa- 
tial distribution of M 3 1 late type stars (/. e. enhanced density 
towards the centre). Of the 14 identifications made of optical no- 
vae, four were presented in more detail. 

The 25 identified and 3 1 classified SNRs had XID luminosi- 
ties between 1.1x10^^ erg s"^ and 4.3x10^^ erg s"^. Three of the 
25 identified SNRs were detected for the first time in X-rays. 
For one SNR the ROSAT classification can be confirmed. Six of 
the SNR candidates were selected from correlations with sources 
in SNR catalogues from the literature. As these six sources had 
rather "hard" hardness ratios they are good candidates for "ple- 
rions". An investigation of the spatial distribution showed that 
most SNRs and candidates are located in regions of enhanced 
star formation, especially along the lOkpc dust ring in M 31. 
This connection between SNRs and star forming regions, im- 
plies that most of the remnants are from type II supernovae. Most 
of the SNR classifications from previous studies have been con- 
firmed. However, in five cases these classifications are doubtful. 

The population of "hard" M 31 sources mainly consists 
of XRBs. These rather bright sources (XID luminosity range: 
1.0x10^^ erg s"^ to 2.7x10^^ erg s"^) were selected from their 
transient nature or strong long term variability (variability factor 
>10; 10 identified, 26 classified sources). The spectral properties 
of three transient sources were presented in more detail. 

A sub-class of LMXBs is located in globular clusters. They 
were selected from correlations with optical sources included in 
globular cluster catalogues (36 identified, 16 classified sources). 
The XID luminosity of GlCs ranges from 2.3x10^^ erg s"^ to 
1.0x10^^ erg s"^ . The spatial distribution of this source class also 
showed an enhanced concentration to the centre of M 3 1 . 

From optical and X-ray colour-colour diagrams possible 
HMXB candidates were selected. If the sources were bright 
enough, an absorbed power-law model was fitted to the source 
spectra. Two of the candidates had a photon index consistent 
with the photon index range of NS HMXBs. Hence these two 
sources were suggested as new HMXB candidates. 

Follow-up studies in the optical as well as in radio are in 
progress or are planned. They will allow us to increase the 
number of identified sources and help us to classify or identify 
sources which can up to now only be classified as <hard> or are 
without any classification. 

This work focused on the overall properties of the source 
population of individual classes and gave us deeper insights 
into the long-term variability, spatial and flux distribution of the 
sources in the field of M 3 1 and thus helped us to improve our 
understanding of the X-ray source population of M 3 1 . 
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Source 


fvar 


svar 


fmax* 


efmax* 


class 


Comment^ 


523 


692.64 


63.33 


106.61 


1.68 


XRB 




1 032 


660.61 


54.21 


33.89 


0.62 


<G1C> 


l(t, 92.2), 11, 16, 27(831.10) 


57 


644.03 


96.796 


147.76 


1.52 


XRB 




1 131 


558.48 


80.18 


38.93 


0.48 


<XRB> 


l(t, 954.2), 2(t, 2163), 19, 20(t), 22, 26, 27(624.05) 


944 


236.47 


44.74 


24.24 


0.54 


XRB 


l(t, 370.5), 14(t, BH-XRN), 20(t), 21(v,t), 26, 27(353.67) 


705 


185.02 


39.85 


32.52 


0.81 


<XRB> 


3(t), 27(178.79) 


1007 


161.84 


42.33 


44.30 


1.04 


<XRB> 


1, 12, 13, 26, 27(12.96) 


1 195 


144.47 


36.95 


28.51 


0.76 


<SSS> 


2(t, 694), 14(v,t), 20(t), 26, 27(82.68) 


1 152 


123.95 


55.31 


31.86 


0.56 


<XRB> 


l(t, 96.3), 2(t, 107), 9, ll(v), 12, 13(v), 19, 21(v,t), 26, 
27(131.73) 


788 


117.83 


20.99 


12.27 


0.58 


<XRB> 


l(t, 20.8), 2(t, 93), 11, 14(t), 20(t), 26, 27(97.89) 


1 136 


87.84 


25.00 


9.79 


0.39 


<XRB> 


l(t, 46.1), 2(t, 155), 19, 20(t), 23, 26, 27(97.65) 


934 


84.94 


28.35 


21.26 


0.90 


<SSS> 




949 


80.16 


25.55 


17.75 


0.68 


<XRB> 


2(t, 13), 20, 25, 26 


1 194 


75.22 


39.69 


12.79 


0.31 


<SSS> 


1, 2(t, 96), 9, 12, 13(v), 14(v), 19(v), 21(v), 26, 27(85.93) 


1 416 


72.25 


21.00 


10.51 


0.49 


<SSS> 


26, 28 


1 084 


59.11 


8.70 


6.64 


0.75 


<XRB> 


l(t, 79.0), 2(t, 260), 14(t), 18, 20(t), 21(v,t), 26, 27(57.41) 


985 


51.66 


80.00 


80.46 


0.97 


XRB 


3(t), 5(t, LMXB), 9, 21(v), 27(76.27) 


1 017 


49.85 


20.73 


9.64 


0.45 


<XRB> 


l(t, 99.5), 2(t, 158), 21(v,t), 26, 27(38.07) 


990 


48.90 


61.77 


82.16 


1.28 


XRB 


l(t, 468.8), 2(t, 285), 14(v,t, BH-XRT), 18(t), 21(v,t), 26, 27(201.71) 


1 059 


48.10 


23.64 


9.16 


0.37 


<XRB> 


l(t, 64.6), 27(39.69) 


698 


41.20 


18.73 


9.33 


1.37 


<XRB> 


13, 23, 26, 27(44.67) 


884 


32.88 


9.59 


2.84 


0.61 


<XRB> 


1, 27(24.95) 


1 153 


31.06 


17.96 


4.22 


0.22 


XRB 


l(t), 3(t), 6(t, LMXB), 27(30.50) 


1000 


26.86 


20.17 


5.75 


0.48 




l(t), 19, 21, 23, 26, 27(34.47) 


1 024 


25.66 


78.46 


144.02 


1.89 


<XRB> 


1, 9, 11, 12, 13, 19, 21(v), 26, 27(34.02) 


1 177 


24.00 


19.31 


4.59 


0.22 


XRB 


3(t), 27(27.06) 


939 


22.48 


10.13 


2.20 


0.20 


<XRB> 


l(t, 65.2), 27(28.02) 


960 


21.56 


9.497 


5.51 


0.55 


<fg Star> 


1, 12, 13, 19, 21, 23, 26, 27(21.97) 


1 180 


16.48 


17.16 


7.53 


0.40 


<XRB> 


1, 12, 19, 21(v), 26, 27(13.01) 


378 


16.21 


9.39 


6.83 


0.98 


<XRB> 


25, 26 


714 


14.89 


13.77 


10.35 


0.68 


<fg Star> 




748 


14.45 


13.86 


3.62 


0.23 


<sss> 


25, 26, 27(14.18) 


92 


14.23 


5.73 


1.86 


0.30 


<sss> 




814 


13.66 


7.98 


4.94 


0.56 


<XRB> 


1, 21(v), 25, 26, 27(15.76) 


1 006 


11.30 


14.34 


1.44 


0.08 


<SSS> 


2(t, 51), 19, 23, 26, 27(10.70) 


542 


10.90 


3.26 


2.50 


0.69 


<XRB> 




1 016 


10.46 


20.66 


6.37 


0.26 


<XRB> 


1, 11 ,12, 19, 21, 26, 27(11.69) 


1 144 


10.23 


10.53 


1.54 


0.12 


<SSS> 


2(t, 38), 26, 27(10.05) 


1034 


10.11 


28.93 


9.81 


0.28 


<XRB> 


1, 19, 21, 23, 27(10.24) 


1099 


9.95 


4.50 


5.12 


1.01 


<hard> 


25, 26 


904 


9.16 


8.65 


8.55 


0.84 


<hard> 


12, 26 


872 


9.13 


17.16 


7.31 


0.38 


<G1C> 


1, 7, 19 , 21(v), 26, 27(7.48) 


1422 


9.06 


2.99 


1.87 


0.55 


<hard> 




1 183 


7.66 


11.93 


8.43 


0.57 


<hard> 


12, 13, 26 


422 


7.23 


9.65 


2.10 


0.16 


<hard> 


13 


823 


7.16 


28.29 


20.25 


1.30 


<G1C> 


1, 4, 11, 14(v), 15, 19, 21(v), 26, 27(6.45) 


1 250 


6.94 


2.80 


0.89 


0.27 


<SSS> 




398 


6.79 


2.93 


1.79 


0.52 


<hard> 


13, 26 


1 266 


6.67 


7.54 


2.37 


0.24 


QIC 


4, 26, 27(8.27) 


975 


6.58 


7.98 


1.76 


0.35 


QIC 


1, 4, 13, 15, 19, 20, 21(v), 26 


522 


6.47 


6.62 


2.21 


0.26 


<hard> 




805 


6.22 


3.66 


0.96 


0.21 


<hard> 




916 


6.21 


3.73 


2.00 


0.44 


<hard> 




862 


6.17 


5.33 


1.38 


0.20 


QIC 


1, 4, 15, 21, 26, 27(5.35) 


1 124 


6.04 


12.93 


5.73 


0.23 


<G1C> 


1, 19, 20, 21(v), 26 


430 


5.98 


4.95 


1.01 


0.16 


<hard> 




1 825 


5.57 


6.11 


2.97 


0.39 


fg Star 


13, 26 


1494 


5.55 


4.06 


0.95 


0.18 


<hard> 


26 


624 


5.54 


7.58 


1.49 


0.13 


<fg Star> 




1 167 


5.44 


7.69 


3.22 


0.57 


<hard> 


1, 19, 20, 21(v), 26 


1450 


5.44 


6.51 


2.24 


0.24 


<hard> 




66 


5.39 


3.80 


2.31 


0.49 


<hard> 




1655 


5.34 


4.72 


1.50 


0.24 


<hard> 




964 


5.30 


6.36 


1.09 


0.12 


<hard> 




226 


5.28 


3.29 


0.56 


0.13 


<hard> 




244 


5.27 


9.08 


3.33 


0.28 


<hard> 


13, 25, 26 


1361 


5.21 


3.12 


1.36 


0.34 


<hard> 


27 


933 


5.19 


4.99 


8.73 


1.30 


QIC 


4, 11, 12, 25, 26 


1366 


5.13 


3.02 


1.05 


0.27 


<hard> 





Notes: 

* : maximum flux and error in units of 1x10"^"^ erg cm"^ s"^ or maximum luminosity and error in units of 7.3x10^^ erg s"^ 
: class according to Table 6 

^ : 1: Voss & Gilfanov (2007), 2: Williams et al. (2006b), 3: Trudolyubov et al. (2006), 4: Trudolyubov & Priedhorsky (2004), 5: Williams et al. (2006a), 6: Williams et al. (2005b), 7: 
Pietsch & Haberl (2005), 8: Kong et al. (2003), 9: Trinchieri & Fabbiano (1991), 10: Collura et al. (1990), 11: Primini et al. (1993), 12: Supper et al. (1997), 13: Supper et al. (2001), 
14: Osborne et al. (2001), 15: Di Stefano et al. (2002), 16: Fan et al. (2005), 17: Pietsch et al. (2007), 18: Garcia et al. (2000), 19: Kaaret (2002), 20: Williams et al. (2004), 21: Kong 
et al. (2002b), 22: Williams et al. (2005a), 23: Di Stefano et al. (2004), 24: Barnard et al. (2003), 25: Shaw Greening et al. (2009), 26: Pietsch et al. (2005b), 27: Stiele et al. (2008), 28: 
Trudolyubov et al. (2002) ; t: transient, v: variable, sv: spectrally variable, r: recurrent, d: dipping, z: Z-source candidate; BH: black hole, XRN: X-ray nova, XRT: X-ray transient, LMXB: 
low mass X-ray binary, NS: neutron star; numbers indicate the variability given by the corresponding paper 
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Figure .1. Logarithmically scaled XMM-Newton EPIC low background images integrated in 2"y!1" pixels of the M 31 observations combining PN 
and MOS 1 and MOS 2 cameras in the (0.2-4.5) keV XID band. The data are smoothed with a 2D-Gaussian of FWHM IC, which corresponds to 
the point spread function in the centre area. The images are corrected for unvignetted exposures. Contours in units of 10~^ ct s~^ pix~^ including a 
factor of two smoothing are at (8, 16, 32, 64, 128) in the upper left panel, at (6, 8, 16, 32, 64, 128) in the upper right panel, and at (4, 8, 16, 32, 64, 128) 
in both lower panels. Sources from the XMM LP-total catalogue are marked as 60''x60'' squares. 
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Figure .1. (continued) Contours are at (4, 8, 16, 32, 64, 128) in the upper left panel and lower right panel, at (6, 8, 16, 32, 64, 128) in the upper right 
panel, and at (8, 16, 32, 64, 128) in the lower left panel. 
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Figure .1. (continued) Contours are at (8, 16, 32, 64, 128) in the upper left panel and lower right panel, at (4, 8, 16, 32, 64, 128) in the upper right 
panel, and at (16, 32, 64, 128) in the lower left panel. 
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Figure .1. (continued) Contours are at (4, 8, 16, 32, 64, 128) in both upper panels and the lower left panel, and at (6, 8, 16, 32, 64, 128) in the lower 
right panel. 



H. Stiele et al: The deep XMM-Newton Survey of M 31, Online Material p 5 




ffl M 31 SN2 

1M July 2006 




432 404' 

XMM-Newton^ 5 □ CtJ^S^^' ^-^ 
EPIC PN + M6d ■ r^^jiT ^ 



"""" ° 



J13 ^ J 



-€WC PN + MOS .i?" "R.J ^'•>^J^3?,1s^ ,^ 



„jm □ -'^Jh • i'", 3 a 



u / 



rn H77H ittiSil:' 

OH .-^T^ 



July 2006 



□ M31 SN3 



-40:40:00 fTl ^""^ 

Pi. ^ . gn ^ 

Pig-' □ siSpl-'^" 



December 2007 



W ST;.-,' 




™ rn '— — 

R ° 2. □ H 

_0.2yr.^keV ^„ 
. XMM-Newton |J . ■ ' / 

[EPICPN + MOS 



■ -f^ # In 



Figure .1. (continued) Contours are at: (8, 16, 32, 64, 128) in the upper left and lower left panels, (6, 8, 16, 32, 64, 128) in the upper right panel, and 
at (4, 8, 16, 32, 64, 128) in the lower right panel. 
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Figure .1. (continued) Contours are at (4, 8, 16, 32, 64, 128) in all panels. 
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Figure .1. (continued) Contours are at (4, 8, 16, 32, 64, 128) in all panels. 
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Figure .1. (continued) Contours are at (4, 8, 16, 32, 64, 128) in both upper panels and the lower left panel, and at (4, 8, 16, 32) in the lower right 
panel. The inner area of the image shown in the lower right panel is shown in detail in Fig. .2. 
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Figure .1. (continued) Contours are at (6, 8, 16, 32) in both upper panels, at (8, 16, 32) in the lower left panel, and at (4, 8, 16, 32) in the lower right 
panel. The inner area is shown in detail in Fig. .2. 
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Figure .2. Inner area of M 31 enlarged from Fig. .1. Contours are at (4,8, 16,32,64, 128,256) x 10~^ ct pix~^ including a factor of one 
smoothing. Sources from the large catalogue are marked as 30''x3C squares. The images are ordered as follows: Centre 1 (upper left), Centre 2 
(upper right). Centre 3 (middle left). Centre 4 (middle right) and Centre B (lower left). 



